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SYNOPSIS 


LAGRANGIAN SIMILARITY TOEORY APPLIED TO DIFFUSION IN 
SURFACE LAYER OF THE ATMOSPHERE 

P» Swadas 

Depaxtaent of Civil Enginaaxing 
Indian Institute of TechnologY, Kanpur 
June, 1974 

Diffusion models based on gradient-transfer theory or qn 
statistical theory (such as 6f Bosanquet and Pearson, Sutton, 
Pasquill, etc#) are inawqiate to rationally and corapletely ■ 
describe the diffusion of pollutants released front a source in 
the surface layer of the atmosphere# Monin applied Lagrangian 
similarity ttoory to solve this problem. His work, however, did 
not cover practical ranges of investigation for the three - 
unstable, neutral and stable - condition^ of the atmosphere, nor 
was tested on field data# Monin's approach has been extended in 
this work to cover all these ranges# Besides, unlike Honin' s 
approach, where he assumed that the equations of moticn repress* 
nted boundaries of the plume, the author has proposed that the 
equations of motion be tidcen to reinresent the mean height of 
ensemble of particles# Equations have been derived mathematically 
relating mean vertical spread z, mean downwind distance travelled 
X , surface roughness laagth and Mcnin-Obudihov stability length 
L« Solving these equations with the help of a digital oeiifatter, 
tables and graphs are presented to facilitate estimation of 
vertical j^reads under given meteorological conditicms# 


The mathematical eqtkitions proposed were tested cai field 
data from two comprehensive diffusion studies : one carried out 
in open unobstructed terrain, viz,. Project Prairie Grass, O’Neill, 
Nebraska, U.S.A,, and the other conducted in urban area ot 
metropolitan St* Louis, Missouri, U.S,A, It was observed that 
there was a very close agreemont in the observed vertical spreads 
and that computed by the proposed equations for both mstable as , 
well as neutral conditions of the atnwsphere for both the cases 
of low and high surface rougjhnesses. In comparison the extended 
Nbnin's approach gave smaller values of wrtical spreads iar 
both these stability conditions for both surface roughnesses, 
hence it may overestimate the concentration distribution* During 
stable atmospheric conditions, estimates of vertical spreads 
by both the approaches (Honin’ s as well as the proposed by author) 
were nearly equal, and generally agreed well with the observed 
spreads for z/L <,0.1, beyond which the agreement was poor, Once 
the vertical spreads were estimated for any diffuslcn situation, 
concentration distribution could be easily canputed assuming 

e3q)onQntial distribution of particles, 

% 

The mathematical models were also tested on lindttd 
data obtained from the city of Kanpur and they gmaerally agreed 
well with the observed concentrations. 
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Chapter One 


1 . IWiaODUCTION 

This God-iiiade space ship Earth is enveloped with a thin and frail 
veil of ain that is pre-requisite to the life, For too long men have 
treated the atmosphere as a sink, dumping vast quantities of efflaonts * 
into it, oblivious to the fact that the air is not an unlimited resource 
and that the ability of air to cleans© itself of poison is not infanite. ’ 

Of all the threats to our onviironment, air pollution is potcnticuLly 
the most dangerous, because whereas wo may replenish soil or purify water 

before use, vvo breathe the air as it comes to us - and each day it oomes : 

i 

to us more heavily l.*adcn with the by-product wastes of our industrialized 
society. 

Abundant evidence has establi^od that air pollution adversely 

affects man and his environment. It is a definite factor in human and i 
animal morbidity and mortality. It damages 'vegetation, accelerates 

I 

corrosion and deterioration of materials, soils property, and interferes ; 
with the oomfortable enjoyment of life and property. It affects climate, 
reduces visibility and solar radiation and contributes to safety hazards, ; 
In spite of the overwhelming evidence of doloteilous effects of air 
pollution, complete elimination of the sources of air pollution is not 
possible at iiio present or in the foreseeable future. It is woven 
throughout the fabric of our daily lives. 

Major proportion of our present-day air pollution is contributed 
by urbffl.ization and industrialization. 'Bie pace of urbanization is very 
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great among 'both the developed as well as the developing countries. This 
has increased the levels of air pollutants, and it will continue to do so 
unless something is done about it. The existing sources of air pollution 
should be judiciously controlled and the future development should be 
rationally planned. The effects of air pollution are obviously proportioi 
to the concentrations of pollutants prevailing in the ambient air* The 
concentrations, among oiiier things, will largely depend on the ability of 
the atmosphere to disperse the xiollutants. Thus, to prescribe a degree of 
control in any nir pollution control programme one must take into 
consideration the dispersion characteristics of the atmosphere. This has ; 

I 

led to the search for on adequate mathematical model of diffusion pattem| 
of air pollutonta over the urban areas. Quite a few diffusion theories | 
have been proicsed for individual sources, either instantaneous or I 

continuous. ; 

The study reported herein is concerned mth the application of 
the L'lgrangian similarity theory to predict the spreads of the pollutants 
in the field. Both the adiabatic and diabatic stability conditions of 
the atmosphere axe studied. 

The effects of the variation of aerodynamic roughness on the 
spread, and hence on the concentration, of pollutants have also been 


studi ed 



Chapter Two 


2. LITERiiTaBE BBVIEf 

2*1 Transport of Pollutants in Air 

Just as in the oase of a polluted streaa or river, a polluted 

air environnent undergoes self -purification. The self -purification of 

a strean is priaarily the result of biological action and dilution. In 

the atmosphere, several self -purification mechanisms, such as sedimentatio| 

of particulate matter, washing action of rain, photochemical reactions, 

absorption by vegetation, soil, etc., are at work along with the dilution ;! 

of pollutants. Unlike stream pollution, where the water flow confines : 

within certain finite boundaries, the spread of air pollutants is tliree | 

dimensional and may not confane within finite boundaries. The general 

area into which pollutants released from a source move is specified by E 

the wind direction, which itself changes from time to time, 

■ 

The dilution of air pollutants in the atmosphere is of prime ! 
importance in the prevention of undesirable levels of pollution, and it j 
is due to two physically separable processes of advection and diffusion i 
(Monin, 1959S')* iidvection brings about the transport of pollutants due 'to I 
mean wind. Other things being equal, the concentration of a pollutant | 
downwind from a source varies inversely with the wind. speed (Howson, 1955) 4 

i 

Diffusion transfers the pollutant from regions of higher concentration to 
regions of lower concentration, and proceeds by two different processes ’ 
of molecular diffusion and eddy dispersion. Molecular diffusion, caused > 

• ' i 

by Brownian motion of molecules, is several thousands of times slower than ! 




eddy diffasion or dispersion, which is due to the transport hy fluid 
eddies. Hence, in atmospheric diffusion problems molecular diffusion 
is often neglected (Sutton, 1953). 

Atmospheric turbulence produces eddy dispersion, and the study of 
atmospheric diffusion problems has closely followed developments in the 
study of turbulence in general. i 

Atmospheric turbulence is classified as mechanical or thermal 
as they aro produced primarily by shearing stress or by convection ( 

(Wanta, 1968). While both types are often present in combination, one ! 
or the other may dominate under particular meteorological conditions* | 

Mechanical turbulence prevails in stable air, whereas in unstable air ! 

|. 

thermal turbulence generally predominates, 'j 

Turbulence increases with increase of wind shear, surface roughnesi 
and instability (Gifford, 1968). High winds and unstable thermal | 

stratification favour rapid diffusion, whereas inversion conditions of 

. i: 

the atmosphere along with li^t winds impede the diffusion of pollutants, i 

[ 

It is not particularly difficult to describe qualitatively the 
concept of the diffusion of air pollutants in the atmosphere, but it is i 

more difficult to express these facts quantitatively in a form so that i 

predictions Cfm be made* Many theories have been proposed to predict 
the downwind concentration of pollutant in terms of source strengiii 
and meteorological parameters. Meteorological conditions are obviously 
of prime importance in such predictions. The existing concepts regarding 
such parameters and the major theories for interpretation of trensport; of 
pollutants are briefly reviewed in the followa.ng sootions. 
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Me taoro logical Elements Useful for Jiir Pollu.tion Studj 


ito eDdiaustive discussion of the meteorological elements useful 
for the nir i»llution study has been given by Sutton (l 953 )j Priestley 
( 1959 )» and Lumley and Panofsky (1964)* Matveev (1967) has also 
discussed the topics in great details, with special reference to 
contributions of Rissian rese,3rch workers on the topic. 


2 . 2 . 1 . Basic Representation of Turbulent Motion 

Sutton (1953) describe 3 the components of velocity measured at 
a point (x,y,z) by u, v, w. In turbulent flow all three oomponente are 
functions of time as well as of position. The mean velocity, with 
comj.ononts u, v, w, is defined at a fixed ijoint and at a time by 


u w 7= / u dt 

t - 4 t 

0 

1 t +iT 

V « -5 / ° V at 
t^-iT 

- 1 A'"*" 

w ~ r j w dt 


(2.1) 


^ t -^T 

0 

where T is an arbitrary interval called the i^oriod of sampling. 


A fluctuation, or eddy velocity, u’, v' , w* , is defined as tho 
difference between tho total velocity at any instant and the mean veloeity, 
so that 


u' 

a U 

^ u. 



» T 


(2,2) 

w’ 

m 'yy 

»',W 
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2.2,2* Boundary layer 

Following Prandtl’ s boundary layer theory (Goldsteini 1950? 

Rouse, 1959 )? ■fcho atmosphero near the earth's surface is divided into 
the planetary boundary layer close to the surface, and above this, the 
free atmosphere. In the planetary bound jury layer (upto 1 to 1.5 Inn) 
the character of motion is considerably influenced by 'tiie underlying 
surface and turlulont friction (Matveev, I 967 ). Ihe diurnal variations 
in the meteorological elements are pronounced in this layer. 

Within tho boundary layer, a surface layer of the atmosphor© 

(50 - 200 ra) is distinguiehod (Sutton, I955)i iii vfhich tho eddy transport 
of heat and tho eddy viscosity stress can bo oonsidoi^ed constant (i,Q* ! 

independent of height) to a aufficiont accuracy (I’rnndtl and Odotjons, i 

1954-). : 

The laws governing processes in the surface layer are very differd 
from those ojjorating elsewhere in the atmosphero (latveov, 196? )• 1 

state of tho surface layer is very olosoly related to the state of the 
active layer of tho earth's surface. Meteorological elements undergo 

i 

sharp variations in height and time. The vertical ®csdicnts of the 
meteorological elements are one or two orders larger than aloft | 

(Trie at ley, 1959). ' 

■ , " ■ , ■ i; 

2 . 2 , 5 . Vertioal Wind Profile | 

2 *2.5.1 Friction Velocity 

The assumptions that the lowest layer of the atmosphere can be 
regcaded as part of a fully developed turbulent boundary layer and that 
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the eddy shearing stress is invariant with height in the surface layer 
of the atmosphere, provided the basis to define the friction velocity, 

, (Goldstein, 1958) as 

% “ (2.5) 

where t is the constant eddy shearing stress and p is the density of air 
The friction velocity, (constant wilii height), depends on the nature of 
the surface and magnitude of the moan velocity (lumley and PsnofslQr, I964, 

2. 2. 5. 2 Velocity Profile under Miabatio Oonditiona 


For adiabatic conditions of the atmosphere, following the class- ; 
ioal experinents of Nikuradse (Priestley, 1959 )» variation of mean ; 
wind velocity u vdth height z from the earth's surface, is given by ’ 
„ u^ 

u a In z + constant (2.4) 


where k is Von Karmen* s constant and is approximately equal to 
The eepation is also written in the forms (Sutton, 1955) * 

- , z 

0 

- ^ + ^0 , 

or ^ 


0*4. 


i 

(24a) i 

■■ I 
I 

P 

(241) : 


0 

Where is the rnu^neas length (the height above surface at which 
tI « 0) which is relented to the size of the proturbenances, and varies 
from a fraction of a centimetre over smooth surfaces like mudflats, 
ice, etc. to several metres over urban aroas. In Table 2.1 , some 
values of are listed. 



TilBIE 2.1 


ROUGHNESS LENGTHS OP V/iEEOUS SUEPACES 


Type of Surface 

z , cm 

Reference 

Smooth mud flats, ice 

0.001 

Sutton (1953) 

Mown grass, upto 4*5 om 

0.2 - 2.4 

Priestley (l959) 

Long grass, upto 60 - 70 cm 

5.7 -9.0 

Priestley (1959) 

Air field 

2.5 

Matveev ( 1967 ) 

Urban areas 



Copenhagen 

0 

0 

0 

Davenport (1961) 

Tokyo suburbs 

5000 

Davenport (1961) 

Reading 

700 

Pasquill ( 197 On) 

London 

2300 

Pasquill (l97l) 


The roughnosB length, ss^, cannot bo intorproted as a function 
only of the surface elements directly below tho ahoraomoters. In fact 
the roughness parameters, as suggested by O'Brien (l965)» niay vary 
considorably according to the influence of upstream topography* He 
obtained by least -square analysis of tho wind-profile data near tho 
surface of the earth, T© get roughness length, z^, Anisimova et al. 


(1965) use Eq, (2,4a) at two heights with half-hourly avera^ valuos 
of wind velocities, thus give 


In z__ 


Ug In z^ - u^ In Zg 


(2.5) 


For very rough surfaces, as found in urban localities, tho 
height z in Eq. (24) is taken as (z-d) whore d is an arbitrary zero-plane 
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displacement length (Sutton, 1955)» A simple method of approximating 
and d has been given “by Iditzhach (196I, q.uoted by Hanna, 1969) « If A is 
the ratio of the total area of a given region to the area of the buildings 
in that region, and h is the average height of ihe buildings, then 
Kitzbach suggests the laws * 


d 

h 



(2 .6a) 


h 


0 ,-1 .1 

- K A 


(2.6b) 


Ihe parameter A varies from 2 to 10 for cities, consequently, if h is 
about 10 m, the roughness length is betiveon 0,5 and 2m and ihe displace- 
ment length d is between 6m and 9iii‘ 


2.2, 3,3 Wind Profile under Piabatio Conditions 


When the air is being heated or cooled from below, the effects of 
thermal stratification and buoyancy have also to bo considered, so that 
at least one more basic parameter is required before the wind profile can 
be specified (Priestley, 1959)* Eiohardson (1926) number has been used 
for this purpose (Batchelor, 1953) • Ihe Richardson number. Si , is a 
dimensionless ratio of buoyancy forces to inertia forces in the atmoaihorre. 
Thus 


ii 


(31/32 + r) 

(au/3z) 


(2.7) 


where g is tho acceleration due to gravity, I is the absolute tamporaturfe, 
and r is the adiabatic lapse rate (« 1® C/lOO m), Following Monin and 
Obukhov (I954)j Sussian workers have used z/l as the basic parameter in 
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place of the Richardson numhor, where L is the scaling length defined on 

% 

dimonsional grounds as , 

5 

u^ 

Xl a - (2 .e) 

k(g/T)(H/Op p) 

where H is the heat flux in the vertical direction, and c^ the specific 
heat of sir at constant pressure (« 0,24 cal/° 0 - gram). When the 
atmosphere is unstable, both Ri and z/L are negative | and for stable ' 
atmosphere, both are positive (liumley and Panofaky, I 964 )* 


Of the purely empirical extension of the wind profile law to the 
diabatic atmospheri c condition, the most applicable ids that due independen 
tly to Laikhtmm (1944, cited by Priestley, 1959) anil Deacon (1949) Who 
proposed ttie form ' 

' (2.9) 


du 

dz 


A z 


whoro A is independent of z, and 6 is greator or losp than unity in 

\ ^ 

unstable or stable conditions respectively. Deacon’s observations showed 

' . i 

that the magnitude of 6 determined from mean profiles was related to the 
Richardson number. Rider ( 19 , 54 ) found satisfactory agreement of this 
relationship in all but stable conditions. Deacon (l9b7) pointed out that 

I 

6 is not truly constant wit|i height. Pasquill ( 1962a)', observed that the 

2/2 

determination of 6 implies/a knowledge of d u/dz which is difficult to 
achieve with sufficient accuracy. 

itoother empirical relation used to represent the wiM profile is 
by the power law (Frost, 1947; as reported by Sutton, 1955) 


u 


(' 




( 2 , 10 ) 
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where is the mean velocity at a fixed reference height , and p is ■ 
an index which depends on the stability of the atmosphere, frost reported 
different values of p related to the temperatuTO difference between 5 
and 400 ft. Although the power law proposed by Deacon (Eq.(2.9)) is 
different from that given by Bq,. (2,10) (Sutton, 1955)f the difficulties 
in determining p are similar to those in determining 3 , 


Most of the recent progress in the understanding of the wind 
profile has been brou^t about by the Monin - Obukhov (1954) 'similarity 
theory' . According to this theory, there exist, near the ground, a 
friction velocity u* and a scaling length 1 that are essentially invariant 
with height, and the flow characteristics in the surface layer of the 
atmosphere oaa entirely be dotermiined by tliem. On dimensional grounds, 
Monin and Obukhov (1954) proposed the generalized wind profile relation- 
ship as 


du 

dz 


u 


* 



k 


. <i>(zA) 


(2.11) 


where is tho eddy viscosity and i}>(z/l) is the universal function of 
the dimensionless ratio z/l, which, for small values of z/i , they gave 

as 


^(z/L) « 1 + o: z/L 

where « is a universal constant. Thus, for adiabatic conditions, 
tends to infinity and 4 )(z/l) tends to unity, giving 

du 1 

dz " F * z 


(2.12) 

L 


(2*13) 


which, on integration, gives Bq,. (2,5). for diabatio conditions and 
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small values of z/L, Monin - ObuMhov relation gives 

11 „ z “ z 

n - ^ (Iji « -— 2 .) (2*14) 

0 

The constant « v^ras originally assigned the value of 0.6 by 
Monin and ObuKhov (1954) on the basis of experimentally determined wind 
profiles, but there is a controversy regarding ih.e numerical value of «, 
and the question of the numerical value of the constant « and its 
universality can not bo considered to have been settled (Matveev, 1967 ) * 
The value of « depends on the stability conditions ( Talceuclii, 1961). 

For unstable conditiono, the values of suggested are 0,7 (Ellison, 
1957), 4.5 (Panofsky at al., 196O), 5.0 (Taylor, I96O), 0.6 (O'Brien, 
1965), ©to.' ; whereas for stable conditions 6.0 (Taylor, 19^^), 7»0 
(MoVehil, I9625 quoted by lumley and Panofsky, I964), 5*0 (O'Brien, 1965)» 
4.7 (Basinger ot al., 197l)j ©tc. 


Avoiding the evaluation of « , Swinbank ( 19^2 5 19^4 ) suggested 
an cijqjonential w3.nd profile as 
u. 


^ II 1 

dz " kL (1 - 0 xp(-z/'l") ) 


and 


u 


CQS J,\ 
k ^exp z^7^ 


(2.13a) 

(2,14a) 


The wind profile for hi ^ly unstable stratification (free 
oonveotion condition) has been given by Kazsnskii and Monin (Kazanskli, 
1965? Gurvioh, 1965) as 


u <■ "o(*-z/L) + constant 


(2.15) 


whore 0 is a constant, which is evaluated as very near unity (Gifford, 
1962). Qurvioh (1965) t^eoommends 0 = I.3 - 1*4* Pidestley (1959) I'QPo^Jted 
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that the free convection begins at about Hi » -O.OJ. Gfuivich (l965)» 
and Kazanskii (1965) take free convection conditions at Bi » -0.05 to 
-0.05, whereas Taylor (196O) takes at = -0,Q8. 

In general, according to Monin - Obukhov similarity theory, 

u= ~ f(z/L) - f( (2.16) 

where f(z/L) is a universal function of z/L and f(z^/L) is the 
similar function of z^/l. 

O' 

The relation between the Eicharison number K. and the Monin- 
Obukhov stability parameter z/L can bo theoretical establidaed (for 
oxamplo, Kazanskii, 1965) as 

z/L«(yi^).Ei (z/L) (2.17) 

where is the eddy conductivity, the eddy viscosity, and (z/L) 
is same universal function defined in Eq., (2.II), For the unstable 
conditions of the atmosphere, Pandolfo (1966) has suggested the relation 

Kjj - Kg . 4 >(z/l) (2.18) 

Hence, from B<is. (2.1?) and (2.18), for unstable conditions, the stability 
parameters z/L and Hi are equal. This has also been tested experimentally 
by Busingor ot al. (l97l)» 

2, 2 ,4. Turbulent Fluotuations 

Lumley and Panofsky (1964) reported that turbulent fluotuations of 
velocity, which are responsible for the transfer processes in the 
surface layer of the atmosphere, are relatively stable when averaged 
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over a time of ■&© order of 10 to 20 minutes, 

jkjcording to Monin (1962), the joint probability distribution 
for the values of dimensionless fluctuations at a fixed point in space- 
time depends only upon the parameter z/L. He reports that the first 
moments of this distribution are equal to zero, and all the second 
moments can be expressed by means of standard deviations o^/u^ » 

0^/u^ and 0^11* considered as functions of z/L, whore a^, are 
the standard devia.tions of axial, lateral and vertical velocities 
respectively. 

2 . 2 . 4.1 Vortical Fluctuations 

On the basis of the Monin-Obukhov similarity theory, the 
stendard deviation of vertical velocity, o^, is given by on equation of 
the form (Monin, 1959 b) 

» b u^ h(z/l) (2.19) 

where b is a constant and h(z/L) is a universal function of the 
dimonsionloss parameter z/l, Monin ( 1959 b) further showed from the 
equation of energy balance that the function h(z/l) was of the form 

h(z/L) -[1 - J (2,20) 

where <(i (z/l) is as defined earlier in Eq» (2,1l). 

Panofsky md Mo Coxmick (196O) suggest the relation 

0^ « b u^ [<1 >(z/l) - 2.4 (z/L)]'’/^ (2,21) 


on dimensional arguments. 
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Hanna (I968) reports that both the above formulae correspond 
to the ©xperimontal data for z/l ^-5j "bu-li for z/L < -5» the a^eement 
is poor, and hence he suggests the form 

= b u^(l +7 (2,22) 

for unstable atmospheric conditions. 

For adiabatic conditions of the atmosphere, L tends to 
inf inity and 

(2.23) 

thus giving that b is the ratio of the standard deviation of vortical 
velocity under adiabatic conditions to the friction velocity u^ 

(lumloy and Panofsky, 1964). IliQ value of b is in dispute between 
Western and Russiai workers. Panofsky and Mo Coaamick (196O) obtained 
b as 1,25» Pasquill (l962b) as 1.33, Pasquill (1966) as 1.25* 

Monin (1962) reported b as 0.7, Qurvioh as 0,7, Porepelkina as 0*87 
(reported by lumloy and Panofsky, I964). Ibe ostimate of Busch and 
Panofsky (1968) leads to more generally accepted value of b as 1*3* 

2* 24*2 Latoral Fluctuations 

Unlike sr , the standard deviation of lateral velocity, a,., is 
insensitive to ohangos in height (Lumley and Panofsky, 1964)* and 
is given by 

“ Oq ^ (2.24) 

whoro 0g is the standard deviation of azimuth, which is usually 
meaw-TOd, and Og and u are measured at the saa© height* 
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Peoiof slcy and Prasad (1965) related with a non-dimensional 

stability ratio B which they define as 


T -2 ^3z / 


( 2 . 25 ) 


They hasro also shown a graphical relationship between B and z/li. 

Monin (l965)ji’^ 'his, .theoretical analysis of turbulence, gives 
the following relations 


0^ V ^(z/L) 

ol oj <|>Cz/L) - z/L 

2 0^ 0 ^ 

"w % . "0 , */L 

and -:r * -s C—s IJ C j 

al 0; 0; <t>Cz/L) - z/L 

0 0 


Kz/L) - z/L 


( 2 . 26 ) 


( 2 . 27 ) 


whore 0 ^ is tho standard deviation of longitudinal velocity, (z/l) is 

as defined earlier in Eq., ( 2 .II), and 0 ^ , 0 ^ , a are the values of 

o 0 ^0 

0.5 0 5 0 under neutral conditions of the atmosphere, 

u V w 


Lateral Standard Deviation, 0 ^ is insensitive to changes in 
hoi^t, but is very sensitive to changes in stability. Hence 
Cramer (1959) has used 0g(“ Cy/u) as predictor for other charactori sties 
of turbulence. Inoue (1959) has also shorn a relation to relate 
with stability parameters, as 

g o, — (2»24a) 

® u 


where B is a constant, and Og and u are reckoned at same height. 
Prom Eq. (2.l6) 
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2 , 3 - Diffusion Theories 

The problem of turbulent diffusion in the atmosphere has not 
yet been uniquely formulated in the sense that a single basic physical 
model capable of explaining all the significant aspects of the problem 
has not yot been proposed. There axe three altoinative approaches 
available, none of which can bo categorically eliminated from 
consideration since each has axeas of utility (Giffoid, 1968)# 

The theoretical treatment of diffusion in a turbulent field 
has developed along the following main lines * 

1, Gradient -transfer theory 

2. Statistical theory 

5. Lagrongimi Similarity theory. 

In the first, the eddy transfer of material across a plane is 
represented as a product of the gradient of material (normal to tibe 
plane) and an eddy diffusivity. This theory is Eulerian in nature 
in that it considers properties of the fluid motion relative to a 
spatially fixed coordinate system. The second and third theories 
involve statistical descriptions and laws for velocities and 
trajectories of typical particles in the fluid, and thus can be 
dewribad as Lagrangian. 


' (2,24b) 
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2 , 3 . 1 . The Gradient - Transfer Theoay 

In his studies of molecular diffusion, Fiok showed that 
(Goldstein, 1958 ) the diffusive flux of a substance is proportional 
to the gradient of its concentration. Thus, mass flux of a 
material is given "by 


Ji 


3 c 
% 3 i 


( 2 . 28 ) 


where is the coefficient of molecular diffusion, o the concentration 
of material, and i the coordinate direction parallel to flow and 
perpendicular to reference plane. 


inalogous to the molecular diffusion given by Pick, Schmidt 
developed a theory in 1925 for the diffusion in turbulent flow 
(Sutton, 1953)* He assumed that a fluid in turbulent motion has a 
kind of granular structure, in which 'lumps’ or 'eddies’ of fluid' 
break away from the mean motion and lead on independent life for a 
short timej and absorb again into the mean flow at some different 
level. These eddies transfer the material or momentum in the 
body of fluid. 

Sohmidt showed that for a turbulent flow, the instantaneous 
flux of material through a unit area of x-y plane con be given by 

\ H + pc w’ (2.29) 

whore the first term on the right hand side is the contribution 
arising from the molecular agitation, and the second tem shows 
that due to turbulent fluctuations. Taking an instantaneous values 
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of c to Tdo mode up of mem component c and fluctuation o', i.e* 

0 » 0 + 0* , ho got the mem flux as - 

J2,(mem) “ -^ + o' ' (2.30) 

Using the exchange -coefficient hypothesis, the tem o’ w’ 
can he expressed as th© product of a virtual coefficient of diffusion 
and the gradient of the mean entity^ (Schmidt, 1925j reported hy 
Sutton, 1955)* If K is this coefficient, on this hypothesis 

J2,(mean) « (k^ + K) (2.31) 

Usually K » (Sutton, 1953) > hence 

Jjj(mem) “ K ~ (2.32) 

which is analogous to Eq,. (2,28). 

Richardson independently applied similar hypothesis to the 
atmospheric diffusion in I926 (Gifford, I968). 

Using Bq.. (2,52) with a continuity equation, Schmidt derived 
a general equation of diffusion in three diirensions, as 


a JL (K 1% + (K li.) + -i- (K li) 
9t 3X sx** 3y '■ y ay-* 3Z z 3z^ 


(2.33) 


where TS^, E.r& ere the ediy diffeeion eoeffieiehte ih the r>- 

y», and directions. 


The problem of the atmospheric diffusion is reduced to that 
of solving Bq. (2.33) linder appropriate boundary conditions, ®ais 
is called the gradient - transfer theory or tbs K-theory, If IT 
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and are independent of spatial coordinates, the diffusion 100011108 
Piokian (Pasquill, 1962a). 

2 , 5 . 1 . 1 . Limitations of Gradient "■Transfer ITheory 

The entire gradient-transfer theory of diffusion hinges over 
the selection of the eddy coefficients of diffusion, K , Z and K • 

Ju ^ 44 

There are two main points to he noted regarding this. 

The exchange-coefficient hypothesis, as stated earlier, is 
based on the postulate that the rate of transfer of a conservative 
entity is proportional to the gradient of the mean value of ihe 
entity? the exchange coefficient is simply the factor of proportiona- 
lity. Provided that this ho accepted as a definition, no ohjection 
can he raised to the gradient -transfer theory in a mathsmatioal 
sense (Sutton, 1955)» In the physical sense, however, there is no 
real justification of the exchange-coefficient hypothesis. There is 
no defendahle ground for the use of Eq.. (2,51 ) and its use, at the ■ 
present stage, is considered as empirical only (Pasquill, 1970a), 
Priestley (1959) points out that there is no precise physical hasis 
for the use of this assumption as the foundation for the description 
of turhulont diffusion in the atmosphere, and consequently the 
validity of the gradient-trsnsfe r theory "is normally judged from 
the degree of success achieved in ... predicting particular diffusion 
phonomona". C alder (quoted hy G-ifford, I968) studied the sppliegibility 
of the diffusion equation to the atmospheric case and concluded that 
the standard K-theory form, Eq. (2.33), can not he generally valid, 
Hussian workers, e.g. Monin (1959®’)» refer to the gradient-transfer 
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theory as a semi-empirioal theory of diffusion. Pasquill (19703-) 
noted that the fundamental validity of the gradient-transfer approach 
has always heen regarded as a matter of dehate. The hasic nature 
of the gradient-transfer theosy must he kept in mind as the chain of 
deductions from the original equation grows larger and more involved. 

Over and above the theoretical difficulty mentioned above, the 
determination of actual values of the coefficients of diffusion has 
always remained problematic. Many expressions for evaluation of K 
in the atmosphere, both, under adiabatic as well as thermally stratified 
conditions, have been suggested. There is abundant literature on this 
phase, and has been reviewed by Priestley (l959)> Lumley and Panofsky 
(1964) and Matveev (I967). Also reference should be made to the 
recent profiles suggested by Hanna (1968), Pasquill (l970s-)» O'Brien 
(1970), and Cohen et al. (1972). The expressions for K-profiles 
within the surface Isyer of the atmosphea^e have been based on Monin- 
Obukhov (1964) similarity theory, but the expressions for entire 
boundary layer remain empirical. 

Some plausible universal forms have been suggested (Hanna, 

1968} Pasquill, 1970a) which can be used irrespectivo of height and 
stability in the atmosphere, but they have yet not been put to tost 
for diffusion studies, mainly because the expressions require 
estimation of such quantities which require special turbulence 


ffl©asurem®its 
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2. 3.1. 2. Solution of Eg. (2.33) 

Once tho soleotion of K-profile is properly done, the Eq.» 
(2.35) is required to be solved to get mean concentration o at any 
point (x,y,z). Depending upon the form of K-profile, many types of 
solutions of Eq, (2.35) have been attempted, 

Sutton ( 1955 ) states the early attempts of Boherts to solve 
the Eq. (2. 35) t assuming that K^, K^, and axe constants (i.e, 
Piokian) and that the wind velocity u is constant at all points. The 
failure of these exact solutions to confom to the observations 
produced unassailable evidence that oddy diffusion in the atmosphere 
can not bo represented by the Piokian equation. Eoberts also gives 
solutions assuming power-law forms of u and and a line source, 
but it considers neither varying conditions of ground rou^mess nor 
varying thermal stability. Thus it could bo expected to apply only 
on a©3X)dynamically smooth terrain, and at times of neutral conditions 
of atmosphere. 

Using Deacon (l949)'s wind -and Kj^-profiles for both adiabatic 
as well as dlabatic conditions, Oalder (l949) gew'e solutions for 
two-dimensional diffusion. Rounds (l955) extended tho work and 
gave solutions of Eq. (2.35) for different types of sources* 
Limitations of Deacon’s profiles in determining the value of g 
(Bq, (2,9)) reflect in both the above applications. 

Smith (1957) treated the two-dimensional equation, for -tie 
case of an elevated line source, with the con^ugate-power-lawa for 
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■wind and diffacivity. The index of the power-law form depends on -the 
stability of the atmosphere, and has to he decided empirically. 


Ihile the two-dimensional case was being studied, some woAers 
also attempted to solve the case of continuous point source, Davies 
(1950) gave a solution where and u have conjugate-power-law forms, 
and K « z^, whore p is the index, same as that for u. Smith ( 195 ?) 
and Godson (l 958 ) also gave solutions for different sources and 
stability conditions using power-law forms for K and K . All those 

V ^ 

works have tho same limitation of empirical determinaiion of powor-law 
index. 


For a continuous elevated source, Bosanq.uot and Pearson (1956) 
solved the diffudon equation mdcing use of dimensional analysis. The 
average ground level concentration at any point (x,y, 2) according to 
them, is given by 


c Cx,y,z) * 


p q u x'* 


exp 


C-C-^ 


2q2x2 


He 

+ — .)] 
px-'J 


(2.34) 


where Q is the rate at which pollutant is released, H is the effective 

0 

height of the stack, pis the index of vertical diffusion, and q is the 
index of ho li sc ntal diffusion. They gave values of p and q for three 
different turbulence fields. Empirical selection of p and q, and 
constnntancy of p with height are the main limitations of this 
approach. Strictly speaking, this equation can 'be applied only in 
neutral conditions of stability (Pasquill, 1962a). 

Monin-Obukhov ( 1954) similarity theory gave gene rail 2 ed wind 
and profiles within the surface layer of the atmosphere. Using 
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these profiles, an exact solution of Eq.. (2»55) is not availahle, hut 
comjuter methods using finite-difference scheme have been usod to 
solve it (Jaffa, 1967 5 Hino, 1968a} Shir, 1970)# Their attempts 
show that this method is very promising# However, determination of 
-profile heyond the surface layer has remained ahhitrary (Cohen et 
al., 1972). 

Walters (1962), Peters and Klinzing (l97l)» Heines and Peters 
C 1975) gave solutions of Eq., (2.55) for different source positions, 
using the form of the eddy diffusion coefficient as 


xr A ^ 

\ . X 


(2.55) 


where x is the downwind distance from the source, q and n are indices 

which depend upon stability conditions, and are coefficients. 

All the previous workers had usod the forms of K and K which did not 

y 2 

vary with tho downwind distance. But it is always difficult to got the 
universal forms of and of the kind shovm. by Eq. (2,55). 


In general, it can bo noted that the basic theoretical limitation 
of the gradient-transfer theory is in the use of Eq, (2.51), whereas, 
on practical side, the dotemination of K profiles has been empirical 
at the present stage. When the diffusion in tho vertical direction 
is to bo estimated in the surface layer of the atmosphere, the 
mathematical analysis of the gradient-transfer theory bocomes vary 
oomplLoated# 
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Despite its limitations, the gradient-transfer theory has great 
appeal to researdi workers in atmospheric turbulent diffusion* If the 
soleotion of K profile is properly done, this theory does provide many 
useful, praotioal results, in approach to the difficult problem of 
the deposition of polydisperse aerosols has been made via K theory 
(Gif ford, 1968)* This theory has been applied to diffusion problems 
involving non-linear chemical reactions and thus simulating the 
concentrations of ohemioally reactive pollutants (Lamb, 1975). For 
fflultiplo-souroe urban diffusion problem, the gradient -transfer theory 
han boen used and both the analytical solutions (Gifford, 19705 
Gifford and Ilnnna, 1970 ) «id oomputor-oriented numerical solutions 
(Randoraon, I97O) ktwe boon offered. 

2.5*2. The S t at iatioal Theory 

The statistical theory of turbulent diffusion, based on G# I, 
Taylor’s pioneer woric (Taylor, I921), studios the histories of the 
motion of individual fluid particles and tries to determine from these 
the statistical properties ncjoesaary to represent diffusion* 

In a turbulcmt field, the diffusive spread of tagged particlos 
from a oouroo is a rmdom process. Taylor oonaiderod idio displacement 

X of a single tagged particle as it is carried away from the origin 

' ' ' ' 

by a turbulent wind u. The mean-square of the displacement, x , for 
n homogeneous and statiomry field of turbulence is given by 

3j2 , 2 u'2 RCC) d? 

0 


(2.56) 
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where u» is the moan-square of eddy velocity and r(^) is the 
Ln^r(jngian oorrolation ooefficiont betwoon the velocity of the single 
particlo at timo t and that at time (t + 5}, Ehe latter is given by 


Hi) * 


u’ct+s) 

..->- 1 .^ 1 . .... I I ■ ■ ■ I-Iii ' J ^ 

{u'^Ct+C) 


(2.37) 


The quantity (x^)'^is the one-dimmsional standard deviation, a , of 

♦JC 

the probability distribution of the particle positions after some 
given time t, Tlio Analysis is extended to three-dimensions by 
introducing Ho|im.*ato oxprossion for (y^)’^ and (z^)^, 

Tho basic assumption on which the statistical theory progresses 
is that tho turbulence is homogeneous and stationary. This idea of 
turbulence homogeneity is a simplification introduced into tho theory 
to avoid the mathomatioal complications (Rouse, 1959 )• Dae Lagrangian 
oorrolation cooffioient, R( C ) of Eq. (2.57 )» is baaed on this 
assumption. In the real atmosphere, this assumption of turbulence 
homogunoity in not valid. Turbulence in tho upper portion of the 
pltnetary boundary layer may approximate the homogeneous typo, bit 
mirfaco Ijyer turbulence is decidedly inhomogeneous (Gifford, 1968)* 
Tho use of this theory, therefore, wilhin tho surface layer of the 
fttmoophoro is alwnya questionable (Randerson, 1970 ). 

2.5»2#1. Sutton* 0 iPCTXiQch 

Tho fora of R(C) mustbe ostabli^ed before the x can be 
evaluated, Prenkiel (l953) has oeaxsidered various possible functional 
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forms of R(5)» All those foms retjuiro prior evaluation of time-scale 


* / RCO ds (2,58) 

which is rarely avoilahlo. in explicit fomulation of R(c) is proposed 
by Sutton (l934)» which later he generalizes to the case of aerodyna- 
mically rough surfaces (Sutton, I955), as 

R(5) « ( L-- .)” (2.59) 

N 4- • 5 

who, re ho defines N as maorovisooaity (» z^)i aiid n depends upon 

Btability conditions* Using this form of R(5,)» Sutton derives 

" 2" "”2 

exprosaions for x , y and z , 


Sutton’ 8 approach was questioned on two main points. One, the 
study is applied to the surface layer of the atmosphere where the 
assumption of turbulent homogeneity is definitely not applicable. 
Second, the integral of the explicit form of R(5 )j Eq. (2.59)» is 
not convergent (Pasquill, 1962a), Moreover, Sutton uses the index 
n no tho stability factor in describing R(s) (and hence x , etc). It 
wna originally assumed that n oould bo determined from the power 
relation 


' 111 *) 


z n/(2-n) 


( 240 )- 


lui 2 2 


whore tho subscripts refer to the different elevations. But l(5) 
assumes a homogenoous turbulent field that is true for Eq. (2,40) 


only if n » 0 (Gifford, 1968). 



Sutton’ s model has been widely used in spite of the theoretical 
diffioultios mentioned above. It should certainly be regarded as a, 
model which is theoretically oriented and physically motivated. But 
it should not bo accorded tli© unequivocal status of a law of nature. 
Paoquill (l9Cjl) observes that those formulae are reliable for 
specifying the average distribution, over a few hundi^ed metros down-* 
wind of a source operating for a few minutes on level unobstructed 
terrain, with a steady wind direction and neutral conditions of 
atmospheric stabilityi but extension to other oircumstanoes has 
depended on empirioal and often speculative adjustaents of Iho 
diffusion paromotors, 

Sutton’ s sohome has been applied empirically to greater 
distances by introducing the separate stability parametorsi n^ and 
rij, , for each direction (Leonard, 1959? Barad and Hsugon, 1959)* 

2,5*2,2, Pasquill’s Method 

"'"'2 """^2 

Alternative approach in evaluating x , y and z (or, 0 ^, 
and <^^)t avoiding the determination of the form of R(5), has been 
suggested. Hay and Pasquill (1957) have presented experimental 
evidence that the vertical distribution of spreading particles from 
an, olwatud point is related to the standard deviation of the wind 
elevation angle, 0 . , at the point of release, Cramer (1959) used 

spootral oonsiderations to link 0 ^ and to the standard deviations 
of Wind vano fluctuations in azimuth and elevation ( 0 ^ aid 0 ^ ), 
respootivoly, by 
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Oy » Oq and 0^ = x® (2,41) 

where x ia the distance downwind and f and g depend on statiility 
conditions. Hay and Pasquill (1959) showed that 

Oy ■ X (2,42) 

where is the standard deviation of azimth, averaged on 

moving intervals} it is the maxinum angle of lateral wind fluctuations 
averaged over period s - T/g , where T is the travel time and t is 
sflffiipllnf' time of pollutant* Hino (196811) has shown a method of 
evaluating tlie value of . He has given a relationship 

between the ratio of > sampling time x aid the averaging 

time s(»T/6)* 

On the basis of available data md guided by theoretical 
expectations, Pasquill (1961) proposed a method for estimating 
diffusion when detailed wind data are not available. !Ehe standard 
deviations of tiio lateral and vertical distributions are expressed 
in terms of an angular lateral spread 6 and a vertical spread h, 
defined by concentrations one -tenth of the axial or ground -level 
values rospeottvely. Pasquill tabulated a weather classification 
aid gave families of curves of h and 6 as a function of distance 
downwind x, for six stability categories. Gifford (1961) oonverted 
Pasquill* s values of angular spread and'vertical spi^ead into stmdaa;d 
deviations of plume concentration distribution, 0 and 0„. Gifford's 
work was based primaally on diffusion data cited by, Pasquill (1961). 



'Pu.m.sr (1964) has modified the Pasofuill’ s stability categoriesf so 
also has been done by Meteorological Agency of Japan (Hino, 1968b), 

Pasq,uill's stability categories are qualitative, and hence 
they can not directly be related to unique values of stability 
length L (Pmiofsky and Prasad, 1965). Hino (1968b) has shorn the 
ranges of stability length L for different Pasquill's stability 
categories, 

Pasquill-Grif ford -’Turner method of estimating a (x) and a (x) 

. y ' ' ■ 2i ' 

has been very widely used in the last decade (Clarke, 19645 Miller 
and Holzworth, 1967 l Koolger et al. I9675 Bowne, 19691 Roberts et 
al», 197 0; Fortak, 1970} Hanna, 1971? Martin, 197 1)» ®ae estimates 
of Oy(x) and ^^(x) wore given from diffusion studies carried out 
at Porton, England, Theso studies were conducted in the open country 
witti tiie avorqge surface rou^ness of 5 cm (Pasquill, 1966). Biese 
graphs, therefore, are applicable to diffusion in open countries only. 
In tho original paper, Pasquill (196 1) observes that tho treatment 
has only boon verified for grassland, but implication is that 
increases in roughness of a moderate order (for example, due to 
crops, Bm;ill budaos and hedges) are unlikely to make avast difference 
to the vortical diffusion; however, for large disturbances such as 
those cm sod by woods, buildings or sharp changes in contour, this 
method Q£in not bo representative. Despite this, urban pollution 
modulo derived in the past decade (Turner, I9643 Bowne 19695 Roberts 
et al. 1970, etc,) have been derived from Pasquill's original estimates 
of spread, without any serious attempts to modify his estimate for the 
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effects of increased urban rou^ess. This mad© Pas(iuill ( 1970b) to 
observe s 

"Some modols appear to generalize about the standard dsfViation of 

particle spread in •vertical and hori®ntal directions, a and by 

y z' 

using pranticnl observations of the distribution of urban air 
pollutants. Although not unreasonable, it should be emphasised 
■that aich a procedure cannot be accepted ecu tomo tic ally as justifioa'bion 
for applying the model to flow and terrain conditions other than 
those obtaining in the pollutant survey. For the most useful 
generalization we iruat seek independent specifications of and 
in terms of me tooro logical parameters that con either be measured at 
the tlB» in question or estimated from routine meteorological data" 

The s'tudioa carried out at St. Louis is perhaps very indicative 
Mo Elroy and Pooler (1966) report their analysis of diffusion 
es^oiiments carried out nt the urban .area of St. Louis. Their 
Pig. 16 (a and b) (reproduced here as Pig, 2,l) compares their 
results according to Pasqui 11-Gif f ord stabili-ty classes wi-th those 
given by Pasqui 11-Gif ford (196I) for the open country. The effective 
surface rougtoess for the St. Louis city is reported to be between 
1 and 5 m (Pasquill, 1970a), which is of the same order of rou^mess 
observed in other cities. It is observed that in the range of 
700 - 100CX) m distance from the source, the vertical spreads in 
the St. Louis city were 5-4 tiiois more than those in the open country, 
except in ■the oaso of B category of stability whore also the vertioal 
spread was about 2 times that for open country at 500 m distance* 
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The lateral spreads, a , also showed tho effect of increased urban 

roughness. However, in this case, the effect of increase in surface 

roughness was not as proncunoed as in tho case of vertical spread. 

Iround TOO m distance the spread was about 2 times which reduced to 

# 

about 1.2 times around. 5OOO m. 

Althou^ the concept that the enhanced roughness in an urbsn 
area would emse an increase in vertical spread, as compared with 
open nnd relatively smooth country, seems to have been accepted 
(Pasquill, 1966), very little thooretioal estimate of the magnitude 
of the ©ffoot has yet boon offeMd. 0?his point is referred to 
iigain in tho following section, but point in order here is that 
tho ostlmates of and do vary with the surface rou^neas. 

Tho estimate of lateral spread is based on Eq.« (2,42), 

but o„ is not constant with tho hoi^t, hence a should vary 

V y 

with tho height, Boemse of those reasons, PasquiU-^ifford 
ostimatos can. not bo used indiscriminately, 

Singor and Smith (1966), Singer, Erizzola and Smith (1966) 
gme the power law relation for a and o„ in the form of 

o-bx*^ (2,45) 

whoro tho coofficient b and power q. vary according to tho meteoro- 
logical conditions and the scale of the problem, 

2, 3. 2. 3. Oonoentration Distribution 

Once tho standard deviations of oonoentration distribution 

or Oy, Og) are determined, the average concentration of 
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pollutant at any distance from the centre of the cloud or aDds of 
tho plumo may ho assumed to he given, hy an exponential function 
(PasqpiH, 1962a), The Gaussian distribution has been assumed 4on 
a continuous-source diffusion model hy Sutton (1954), Prenkiel (1953), 
A.S.M.E. (1968), Turner (1970) and many other iirorkers. Thus, for 
a continuous point source, tho average concentration 0 at any point 
(x,y, z) from the source is given hy 

2 2 

cCx,y,z3 " — ~~ — exp [- j (^ + •—)] (244) 

2ir u c_ ^ a_ 

y z y 2 

whore Q, is tho quantity of material roloased per unit time, ond u 
is tho effective moan wind speed invatisnt with tho height. 

Strictly speaking, the Gaussian diffusion model applied only 
in the limit of large diffusion timo and for homogeneous, stationary 
conditions (Gifford, 1968). Tho empirical support for the Gaussian 
shape oomos only from lateral distributions from continuous sources, 
or, of vertical distributions from effectively elevated sources 
(Pasquill, 197 0a), On the other hand, Pasquill (197^^) observes, 
there is no real basis for assuming vertical Gaussian distribution 
from a ground-lovol source. The theoretical and empirioal indications 
point to a ^apo between Gaussian and simple ej^onential, and the 
experimentally determined values of the exponent s in the exponential 
form, ©xp(-bz®), lie in the range of 1 ,1 to 1,5 (Pasquill, 1962a), 

The workers like Bosanquet-Pearson (1936), Calder (quoted by Pasquill, 
1962), etc, have used s <» 1. 
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In general, witiiin the limitations cited above, hoth Satton's 
as woll ns Pasquill's methods are very useful in the diffusion 
studies, "but their use beyond the specified circumstances becomes 
more of empirical nature. 

2. 5*5* The Lagrangian Simlarity Theory 

In the last few years, attempts have been mode to apply the 
so*»Cftllod Logrfmgian similarity principles to the problor of diffusion 
from a source of passive material on the ground. 

Ellison (l959)j prompted by a remark by Batchelor (l959) on 
Monin ’ 0 (1959^) paper, applied a dimensional method to the dotormi* 
nation of the diffusion downwind from a continuous point source in 
an adiabatic surface layer. Ho uses Monin’ s (19591>) (i) similarity 
hypothesis that the average vertical velocity of a particle released 
from a source is proportional to u^ , times some universal dimensionless 
function of Monin-Obukhov stability length L given by Eq. (2.19); 

(ii) Bimilsrity assumption that the average horizontal velocity of 
a particle is given by Eq. (2.16). Ellison, applying the lagrangian 
similarity principles, shows that the probability iji that a particle 
diffusing in the surface layer will reach a distance r » (x +y +Z5 )® 
from its average position characterized by coordinates x, y, z must 
bo a universal function of (x-x)/z , (y-y)/z and (z-z)/z. On 
dimonsionlal reasoning ho ^ows that for a odntinuous point source 
located in a neutral atmosphere, the average axial ground -level 
concentration c is given by 
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0,0) a(JS) 

0=4/ = — (2.45) 

z dz 


whore b is a ■universal constant, and Q, is the rate of emission. On 
simplification, he shows that 


c 

Q ‘ 

where from Eqs. 


bu^« 


-X 1 T 


(2.16) .and (2.19), 


(2.46) 


[z (]n-|- - 1)] (2.47) 

0 

By evaluating Eq,. (2.46), Ellison (1959) showed that in 'the 

ndinbatio surface layer the downwind concentration from a continuous 

point sou roo varies as jc^, where p varies approximately in the range 

-I.S to -1*9. For a continuous, infinite crosswind line source, the 

-1 

downwind concentration was found to vary as x , These results are 
in excellent agreonont' with data, and hence provide an alternative to 
Sutton’s solution that leads to essentially the same results. 

For the diabatic case, Gifford (1962) proposes 


^ j- (248) 

z‘^ u(z) 

which is 'the countej^art of Eq. (2.46). For a relation between axial 
ground -level concentration and mean downwind distance x to be obtained, 
a relation between x and z must be established, which, according to 
Monin (1959b), is given from Eqs. (2,16) and (2.19) as 
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bk - [fCz/L)-f(Zo/L3] 

— X a / dCz/L) 

L zjh h(z/L) 


(2.49) 


The results obtained by Gifford are in reasonably close 
agroesnent with detailed experimental atmospheric concentration 
meamrements. Wind-tunnel diffusion studies by Cermak (1965) have 
provided additional verification. 


The most importmt feature of the Lagrongion similarity theory 
is that it treats surface layer diffusion without postulating a 
diffusivlty. Thus it takes into account the inhomogenaity of 
turbulence in the surface layer of the atmosphere .(unlike the 
statistical theory), but does not assume a diffusivity coefficient 
(as done in -Kho gradient-transfer theory). 


On the basis of the good agreement between the theoretical and 
observed axial concentration values, ttie hypothesis of Lagrangian 
similarity in both the adiabatic and diabatio surface layer conditions 
seems to be strongly supixjrted. However, it has been tested only 
for axial ground-level concentrations. Moreover, in the application 
of this theory to a point source, the rate of increase of the 
lateral spread of Iho particles about their mean position is not 
explicitly stated. GaMer (reported by Pasqui 11, 1966) observes 
that the lateral spread is implicitly proportional to the vertical 
spread. The physical acooptibility of this implication of the 
similerity theory is questioned by Pasq.uill (1966). 
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^ Lagrangian Similarity Theory to Estimate the Vertical Spread 

Although tho Lagrangim similarity theory was put forward to 
predict the variation of a continuous point-source concentration with 
distance downwind, it can also te used to estimate the vertical , 
spread. 

In fact, the sirailarity theory approach to problem of diffusion, 
started by Ellison (l959) ai^d Batchelor (1959)» '’''''s.s stimalated by 
Monin’s ( 195913) Paper on smoke propagation in tho surface layer of 
tbo atmosphere, in which he gave equation of the shape of the 
boundary of the smoko plume, which took into account tho variation 
of wind velocity with altitude. He shows that the height z and 
distance x reached by a particle i?eleased from a source are given by 

||-bu*h(z/L) (2.50) 

||-~tf(z/L) - f(zyL)] (2.51) 


from which 


dx 1 

dz “ bk h( z/L) 


Integrating, 


(2.52) 


bk 


K * 


/ 


z/L If(z/L) 


fCyL)] 


dCz/L) 


(2.53) 


^ \/l hCz/L) 

Monin integrates Iq. (2,53) numerically and the results are expressed 


in a aeries of curves of z/L (upto z/L » 5) against bk three 

values of Gifford (1962) re-evaluated the curves for the 




stable conditions. Monin gives b to be close to or a little 
excooding unity. Gifford does not give any value of b. 


59 


The main point to be noted ■with Monin* s approach is that he 
takes the Eg.* .(2.54) to represent the boundaries of the plaume. 
Unlike infinite Gaussian spread assumed by all Western rese aj^chers , 
he takes the spread to be finite and assumes that the Bq.. (2.55) 
represents the boundary of the spread. 

A slightly modified approach suggested by Panofsky and 
Prasad (1965) makes use of the root-mean-sguare vertical component 
of eddy velocity, a , and thus they take Eg. (2.52) as u/o which 
is proportional to l/o^ , whore is the standard deviation of 

wind elevation angle, henco 


~ fZ 1 . 

X « / ^ dz 

0 

Kg 

Empirical data on as a function of z, z^ and L' ™'hero 

is eddy conductivity and eddy viscosity) have been summarized 
by them, and numerical integration gives x - z curves. 


Pasgulll (1966), trying to examine the extent to which a 
satisfactory representation of vertical spread was provided by 
similarity ■Uaeory, emphasized that z in Egs. (2.50), (2.52) and 
{2.53) is the moan hoi ^t of an ensemble of particles. Hence, 
these eguations become 

||-bu^h(w/L) 

— tfCi/L) - fCz^/L)] 


(2.50a) 

(2.51a) 
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and 


ii . 1_ 

dz h(z/L) 


iff - 

z^L h(z/L) 


He further showed that 


dz 

dt 



dz 



CO 

0 dz 


(2.52a) 


(2.55a) 


(2.54) 


where 0 is the ooncontration of pollutant at any hei^t z, is the 

eddy diffusivity. He proved, on this "basis, b to be eq.ual to k 

(■ 0,4) and got good agreoieat at IOO m distance with observed 

vertical sproai in odiabatio conditions. However, using z as the 

mem height of m ensemble of particles and b = k for the theunallj” 

stratified ool^ditions, he observed theoretical as well as practical 

inaloquacies . For such cases, he demonstrated that vertical spread 

at a distance of lOO m from a source for z = 0,03 cm, was in much 

0 

better agireement of observed values if vertical spread was 3 ^epresented 
by the hei^t of cloud z^^^^ (concentration at is lOfo of axial 
concentration) and b » 1,2, substantially similan to the treatment 
of Monin, 


FaBluill (1970a) gave additional verification for neutral 
conditions of the atmosphere, and showed that the similaxity theory 
could predict the vortical spread, in tae neutral conditions, to a 
oonsiderablo rang® of distance upto 5 Km. 
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This method of piredicting vertical spread of pollutants iwith 
the help of the lagrangian similarity theory has a great potential# 

It makes use of both the important points of the similarity theory, 
viz#, taking the atmo spheric turbulence as inhomogeneous and at the 
same tiHB avoiding use of diffusivity concept, and steers clear of 
the contTOversy of the use of implicit crosswind spread of the 
similari.ty theory in predicting the concentratton downwind. But 
it should be noted that althou^ the treatment has a fairly 
convincing theoretical background, it has so far been tested in a 
very limited way. 

Thus, the prosent-day status of mathematical modelling of 
atfflosihorio diffusion is such that even when both the gradient- 
tronsfor theory and statistical theory are in use, the Pasquill- 
Gifford method has become more popular in the last few years, and 
is used even for urban diffusion studies, despite its limited 
applicability to open, unobstructed terrain only. It is interesting 
to investigate vihother the Lagrangi an similarity theory can be used 
on field conditions to predict vertical spreads and hence concentration 
distribution. Such investigations have boon done in a very limited 
way. The effects of varying aerodynamic rou,^ess of underlying 
surfaces and different stability conditions of atmosphere can also 
bo oxaminod, whioh have been reported in the following ch^ters. 



Chapter Ihree 


DEVBIDSMENfir OP A MATHEMATICAL MODEL IDH PEBDICTION 
OB’ POLEUTAilT COHCMTRAIIONS 

3-1 Loaportance of Mathematical Models 

Tb safeguard against air pollution hazaa^ds to health and 
property^ air pollution control programmes have to be undertaken* 

The concentration of a pollutant at the points of significance 
from communily point of view depends on source strength and the 
dispersal capacity of the intervening atmosphere. Hence, for 
any air-quality -management programme to be rational, simulation 
models of diffusion of air pollutants have to be developed. Such 
diffusion models are the yardsticks for furnishing quantitative 
evaluation of concentrations of pollutants in the ambient air. 

Mathematical diffxision models can be used as an aid to decision 
making in air pollution control, and city and regional planning. 

These models can be used in long-range air-resource-management 
programmes for the evaluation of future air quality wilia regard to 
increased industrialization or increased population. They can be 
used to evaluate the effectiveness of individual control or abatement 
procedures for a given source or combination of sources. Models 
may also be useful for detezmining most suitable locations for 
pollution-causing industrial units as also of sampling stations 
used to monitor air pollution levels. It may be possible to use 
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short-range predictions of meteorological variables with appropriate 
source information to produce, in a few minutes with the help of 
a digital computer, maps of expected air pollution concentrations. 

Biese would be an aid In initiating immediate measures necessary 
to achieve desired air quality levels (Turner, 1964; Middleton, 

1970). 

lfe.|or elaaents of the mathematical diffusion model for urban 
air pollution problem can best be shown through a flow-diagram 
as shown in Fig. 3*1 (ibrtak, 1970). It shows that the mathematical 
model to describe the dispersion of pollutants in the real atmosphere 
plays a very Important role in diffusion modelling. 

3.2 Earlier Models 

To estimate concentrations of pollutants downwind of a source 
in the surface layer of 14ie atmosphere, as discussed earlier, three 
diffusion theories are in use. !Ehe gradient -transfer theoiy and 
lagrangian similarity theory of diffusion have the merit of 
incorporating the turbulence inhomogeneity. However, the 
similarity theory has a doubtful Implication about crosswind spread 
and its use requires evaluating the function ^ , Eq, (2.45),vdiich 
is not always available. 

3*2.1 Bosanquet-Iearson Model . 

Ihe gradient-tisuasfer theory can give good results if; proper 
form of K-profile is available, but it involves complicated math amatical 
analysis. Bosanquet and Iearson‘s equation, Eq. (2.3.4), based 
on olassioal diffusion equation, have been popular in diffusion 
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modelling (Maglll, et. al., 1956) but selection of the horizontal 
and vertical diffusion indices, p and q of Eq. (2.3.4), has been 
empirical, and it is applicable to neutral conditions of stability only. 

3*2*2 Sutton *s Model 

Another equally popular formula used in 
making practical estimates of the distribution of concentration 

has boon that developed by Sutton. His approach, as any other 

approuoh baaed on the statistical theory in general, assumes homo- 

goneoua turbulent field, and is not recommended for the case other 

than neutral atmosphere. 

3.2.3 Baequlll»8 Method 

Over the last decade, Pasqulll's method 
has become very popular because of its practioallty and simplicity 

in UBS. HlB OyCi) and have baan derived fmm the actual 

observations, hence they account for inhomogeneity of turbulence. 
However, the variation of spreads with the surface roughness is not 
accounted for in Easquill's approach. Again, there is a discontinuily 
in incorporating stability conditions. If the spreads could be 
©stinated taking into account these two variations also, the field 
of applicability of this approach wuld become almost universal. 

!Ito this end, at least within the surface layer of the atmosphere, 
the lagranglan similarity theory may be useful in predicting vertical 
spread of pollutant. 

Over and above the effects of enhanced surface roughness on 
vertical spreads, two more points of Easquill's method are to be noted. 
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In fistimating a (x) by Bq, (2.42), the standard deviation of 

j y 

azimuth la taken as invariant with height, and hence a (x) is 

y 

taken as ssraie at all heights, which is not true* Further, in the 

estimation of concentration distribution by Eq. (2.44), u is 

also assmed invariant with height, which is consistent with the 

assumption of homogeneous turbulence, but unjustified in the 

actual almospheric conditions. 

3*2.4 Model based on lagreuagian Similarity 

Bq8.(2.53) & (2»53a), based on similarity theoiy, relate the vertical 

spread of pollutant with the downwind distance x, roughness length 
z , and stability length l. In the use of these equations in 
predicting vertical spread downwind, there is a controversy in as 
much as Moniri maintains, that Bq. (2.53) reckons the spread upto the 
plume boundary and there is no concentration beyond this boundary, 
whereas Paequill argues thht the z in Eq. (2.55) is the ensemble 
mean 5, and Eq. (2.53) becomes Eq. (2.53a). 

The best-developed simulation model of this catego:iy to-date 
is by Monin (l959b) briefly described in Section 2.4. As referred 
to in the section quoted gives equations to compute the upper 
boundary of the plume and also developed curves to directly read off 
z/l if values of bfo^Ii are known. Ihe curves given by Monin 
only cover a range of upto 5 for the unstable and stable conditions 
of the atnosphere. • Since Monin’s approach appears to be one of the 
most promising for application on field, the approach 'has been 
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extended to cover a larger range of ~ and neutral conditions of 

Xi 

atfflospherio atabilily in the later sections. 

Neither of the above approaches have been adequately tested 
on field data for stratified field. 

3*3 Extension of Monin's Approach to Gover Practical Banges of 
z/lt for All Oonditiona , 

lonin ( 1959b), in his wcrk as briefly, referred to in Section 
2.4, assumed that the Integral (Eq, (2.53)) 



Z-Ii 

/ 

VL 


Wis/I) - f(a /l)l , , 

EOT 


(3.1) 


showed the relationship between the vertical spread and downwind 
distance travelled by pollutant particles. He further assumed 
that Bq. (3.1 ) represented the equation of motion of particles at 
the upper boundaiy of plume which was of a finite width. To . 
differentiate Monin’s reference to plume boundary from mean vertical 
spreads, 2 in Bq. (3»l) for Monin’s approach is denoted here as 

«y* ' ' . ' 

Ionia has given a set of graphs of bks^I versus z/l (in 
t)!r0S6iit» 00*80 f for sfo-bl© ss W61X ^0'© misliabl© strEtifioEtioKii 

in his Pig. 3* Die set of z^/I values consists of » 

(i) 0.0010, (ii) 0.010 and (iii) 0.10. Ibr the comparison with 



48 


tho apprufuch presented later in this chapter with mean vertical 
spread, these graphs given by Ifonin need extension to values of 
beyond 5 for the unstable atmospheric conditions. Itor 
both tia® unstable as well as stable conditions, graphs at z^/lf=0.0001 
are required for the comparison purpose. Ibr neutral condilsLons 
of th© aiaaosphere, Monin has not given any graph. 

In the following three paragraphs, the equations and graphs 
are developed using Ifonin 's approach, for unstable, stable and 
neutral atmospheric coalitions. 

3»3*1 Unstable Atmosphere . Sbr the unstable atmosphere, as 
ehovm in Section 2.2.3«3 and later sections, free convection 
occurs at 2 /(“L) » 0.08, and the wind velocity profile can be 
given as 

Zq 

Ibrced convection, i.e,, - -jp- £ - 2/1 £0.08, 

- % 
n B "«»> In 

Kree convection, i.e., 0.08 ^ ~ z/1 

' ' \ ■ / \ 

u » ^ [0.48-2 . 3 log (Zp/-l)-1.3(-2/Ii) i 

Irom eq* (2.16) 

- f(0 A) ® 1*1 2/Zn <,0.08 (3.2) 

& f( 2 /I-)-f(zQA) ^ 0.48-2.31og(::|;)-1*3(-z/l.)“*'^^ for 0.08 < -zA 
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!Qi© standard deviation of vertical .wind velocity fluctuations 
oan be given by 

*= for -z^L <-^ 110.08 

‘’w “ bu^h(z/l) for o .08 < -z/l 

Ibr free oonvoot ion layer, as shown in Sec. 2. 2. 4.1 many expressions 
for hCa/l) (and henoa for are available. Ibr z/(-Ii) > 5, 
th® on3y esqpreasion available is by Hanna (l968) (Bq. (2.22)). 
Adopting his relationship 

- bu^ (1 + 7 2/(-l))^/^ 

h(2/l) » (1+7 ^(-l))’’^® • (3.5) 

Substituting the expressions of Eqs. (5.2) and (3.5) in the 
Eq. (3.1) 


bta 

(- 1 .)“ 


0.08 

2q(-Ii) 


In z/z^ 

"'T' 


d(z/l) + 


/ 

0.08 


0.48-2«31og z^/(-Ii)*'1.3(z/~Ii) 
(l +7 zi/( *.I») j'^^ " 


“1/3 

■ d z/(-Ij) 
(3.4) 


integral of iq, (3.4) was evaluated numerically by using 
Gauss-Iiegendre Quadrature method, for which number of points jused 
were foycct in an interval of ' z/l « 1. computer programae for 

the evaluation of this equation is shown in the appendix A-1. (Ehree 
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casea of different values of Zq/(-l) were studied. 

In the above derivations, it is assumed that free convection 
occurs at « 0.08. When zy(-l) ^0.08, free convection 

conditions can be assumed from itself, In that case 

u»p [1.3 1.3 

f(z/l) - 1.3 1.3 (3.5) 

Bien Eq. (3»4) becomes 


bkx 

(-1) 


^rna/^ 1 . 3 ( z 1 . 3 ( z/-l 

/ . ^ dz/(-l) 

4/(-^) (1+7 2/(»I.))^^® 

(3.6) 


Bq* (3.6) was also evaluated numerically as was Eq. (3.4). 

From Eqs. (3*4) and (3.6), the fallowing liable 3.1 is prepared, 

where bk^(-I() values for different values of 

tabulated, for the oases of z^A-Ii) » (i) 0.0001 , (ii) 0.001 , 

(til) 0,01 and (iv) 0.1. Figure 3.2 graphically presents these 

values. Ihus for unstable atmospheric conditions, a .relationship 

has been obtained between the z„„^ and the mean downwind distance x 

loax 

from the source, if z„ and 1 are known, 

■ 0 . ■ , 
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TABLE 3.1 

TUB VAHJBS OP g^^y (~L) AID COERESIOiroiBG VAIIJES OP 
bk^C-L), IDE DIPFERENT zyC-L), BY USHO EQS. (3.4) 
AND (3.6), IDE UNSTABLE ATMOSPHERIC CONDITIONS 



0.00010 

0.0010 

0.010 

0.10 


blt£/L 



0.03 

0,141 

0.075 

0.013 

^ 

0.10 

0.531 

0.306 

0.131 

- 

0.20 

1.181 

0.244 

0.315 

, 

0.50 

3.144 

2.134 

1.133 

0.204 

0.70 

4*417 

3.053 

1.698 

0.390 

1,00 

6.333 

4.497 

2.544 

0.706 

2 .00 

12.157 

8.785 

5.3O8 

1.789 

5.00 

27.983 

20.566 

13.098 

5.178 

7.00 

37.771 

27.894 

18.018 

7.418 

10,00 

51.762 

38.399 

25.088 

10.706 

20.00 

94.839 

70.866 

47.128 

21.228 

50.00 

208.993 

157.281 

106.428 

50.383 

70.00 

278.729 

210.198 

142.818 

68,558 

100.00 

377.921 

285.487 

194.578 

94.573 



iMMABir ATMOSf’HERE 


0-00010 

0-0010 

^010 


... i* 1x10 ixifr- S; 

tbkx/(-L) if 

%‘bKx/l-L) FOR DIFFERENT VALUES OF 
'.ByLEaS-(3-4)and (3-6) FOI? UNSTABLE 


/ ■' ... ■; o: ^ 
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3 » 5 t 2 i gtablci Atmos-pSaew 

foa? tho atsiiLQ »a.tnio spheric ooMli-tions, 
onn ba giTan by 

. % s 

itL « ^ { lai 4 « — ] 


the Tsinml irofile 


(5.7) 


where ^ is fttmiversall loxaistant. Bia w.alioLe of « , aa flismssed 
in 8601, 2»2.3.3> can, h® “talan as 6*0 I'o:e "tdae stable ai*aiC3®ib.erio 
oondltio'ins# lha seconS, "ierm on the Eiigtot “hand side can fee taken 
na * b/Ii ^ox all prmDtiia.a,l purposes^ iao«ll bhuS| the equs^ftion. can 
ha wxS itwft m 


u. 


u 


1141 f- + « 

lO 


from (2*16) and (3 a) 



(3.7®-) 


f(a/L) - « In f- 4- «:S (5,8) 

O'’'.' 


1?h© Btgendaxd deriation of vQ'Ctjiieal wind valooitliy f liotuations, 
aocoxiing* to Monin iS' given Idy ilXqs, (2*193 snS 1(2.20)) 


o* - I ■’ -4^ 1* 

wh«r®^i(z/L) is as defiwd by Bq.s. (Xill)' and (2.12). 
i{z/L) « 1 

Sahatl tu 1 tng Eq. • ( 3 •« 1 '0} dll E 9 • C 3 ) 


(5.9) 

4l!j3»oxdingly 

(5.10) 


o - l», [t - 

' j ■4* « 


r 


II fc ' 




(5.11) 
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From (5.11) m& (2.20) 

Using Eqs, (3t.8) ani ( to Eq. (5.1) 


hk 



,jh 


JC(s5S^tllnj) f(2; /l) 

— ° 




(3.12) 


(5*13) 


El* (3*13)osnb« ift te gretesdl ntmerio ally. Fox a^/^l'W^ltues 
M (i) 0,0001, (ii) 0,001, (lid) l*h.Ol and (iv) 0,1, iii,€ e^jjui^lifeion 
was evaluatsd numerioaXly u elini 'l•la^ss-Legendze ftuaiaiiatti.ze wtood 
(Caosnahm et al.f 1969 )* Sd^ecsjit' the important values oi btt ^/L 
oor responding to have Tse^s-d. tabulated in Table 5 , 2a^ . 

Fig, 5,5 shows ttie relations ngi®##ilcally. 
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i 

TJffllB 5.2 

TO YAixi»s OP MD oomsioiroiiiG Tjara op 

ak i/L, M DIPPEHENT a^L, BY USIIG BQ. C3*13)f 
FOR miM AIIOSifERIC CONDITIONS 


zjh' 

o.ocoi 

0 

0 

0 

«.0l 

0.1 


bk j^L 



0.01 

0,C46 

0,015 

- 

jm 

0.02 

0.065 

O.Q4O 

0,00^ 


0.05 

0.261 

0 i 14 5 

0.05a 

- 

0,07 

0.590 

0,228 

0.078 

- 

0,10 

0.625 

0.393 

0.161 

m 

0,20 

1.452 

0,993 

0.525 

O.O69 

0,50 

4.590 

5.410 

2.240 

0.900 

0,70 

7.100 

5.450 

3.8IO 

4.875 

1.00 

11.400 

9.100 

£.650 

5.911 

2.00 

50,307 

25.710 

21.000 

15.210 

5.00 

1 !27.100 

115.120 

W5.3IO 

as. 598 






Eq. (3.16) gives the relationship between the z and the 
mean downwind distano® x f»m the source, if the roughness 
length is known, for the neutrally stable conditions of the 
alTO sphere. 

(Dable 3 *3 and RLg. 3 ♦4 show the ®nd corresponding 

values evaluated from Bq, (3.1.6), taking b »= 1.2. 

lABIB 3.3 

THE VAI0K OB z /z^ AMD COHRESK)MDIN(} VAHIES OF Vz , 

max 0 -o’ 

M!0M Bq. (3.16), TOE OTUmAI AIMOSIHERIO COMDITIOITS, 

mms b » 1*2. 


W®0 


10 

29.2 

20 

85.5 

30 

152.0 

50 

304.0 

70 

475.0 

100 

752.0 

150 

1250.0 

200 

1795.0 

500 

5410.0 

700 

8080.0 

1000 

12300.0 ■ 

2000 

27500.0 

3000 

43750.0 

4000 

60900.0 

5000 

78100.0 
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5 *5 *4 Ooiaoentratioa Dlatilbutton by Moala's Approach 


Sbr diteimiaation of the diatribution of pollutant with 
helfihff Mhain ( 1959b) proposed the use of a hyperbo'llc equation, 
in place of common use of parabolic equation. He assumed that 
the plum© had a finite spread, and accordingly use of "telegraph 
equation'’ of hyperbolic nature was consistent with the principle 
of limited propogation velocity. !Ifae eqioations are 


is. 

at 


+ 


M 


- 0 , 


12 

at 


+ 2 a P «» 



(3.17) 


where o is the concentration of pollutant, P is the vertical 
flux, , and a is described as a ’typical frequency’ 

mJ 

of turbulaoce, which is given in the form 



Ibr scxae function i{) of z/I. Bq. (3.17) can be solved numerically 
for any stratifioatioa. lor adiabatic conditions, closed-form 
solution has been given by Monin. 

Basquill ( 1962a) gave a comparison between the concentration 
profile given by Monln’s hyperbolic equation and parabolic equation 
commonly used jtor diffusion studies, and showed almost identical 
results, eacoept near lije top of the spread, vstoere Monin's profile, 
in accordance with the concept of finite propagation velocity, 
suddenly falls to zero concentration. 
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OyltlcaZ Aspraleal of Monin^s Approaoh 

Itoalja assumed that Eq* (3.I) represented the equation ot 
motion of partiolee at the upper bouhdaiy of plume. But looking 
to the Bqs. (2.50) and (2.51 ) by which h(2/I.) ani [ f(2/l)-f(zyii) ] 
are determined, it can be seen that z in Eq. (3.1) shp-uld represent 
mean hei^t of ensemble of partiolee (Basquill, 1966). Although 
the assumption of finite diffusion is more realistic physically, 

Eq. (3.1) should mean to represent mean height, Hie determination 
of oonoentration distribution consistent with finite-width plume 
becomes otanplioated auad estimated' concentrations are not very mtch 
different from those given by simple exponential distributions, 
usually adopted in the diffusion studies. ajhus, it can not be 
denied that fLnlte diffusion is phrsioally more realistic although the 
praotioal difference is probably not great. 

If exponential distribution of concentration is adopted fbr 
simplicity in use, z in Eq, (3.I) i^uld represent some near- 
extreme spread of particles, laaqiiiH (l966) arbitrarily assumed 
that z in Bq. (3*1 ) 'wa® representing the hei^t at vhich concen- 
tration falls down to 10^ of the axial concentration* 

M3ain*s approach was first put fctrth in 1959; Easquill 

(1966) used it to predict vertical spredd* However, very little 

verification has been done on ftejd datSf especially fbr tfaemally 

I , 

stratified field and with large surface ibughness. 
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1.4 

I0 % fidtid of d®v®Wj?lng s ooia|!3?e&6{ti0iy© of 

©t-tooijphtrio iiffuiloia, wblob ©an take into aecomt the effect# 
of boiai the tohoaogeneous turhulehce field of diffejfent sthtiliiy 
conditions of the atmosphere as well as of differeit aet^jdyhiiiic 
rou^atisee of the uMerlSrlng surfaoej and would lend itself to 
sSaplo and reliable lapedlotioni of pollutant eoaeentration at 
various points* 

iEh# present study was taken up to Investl^te the app|i®abllity 
of the lagraaglan similarily theory in predicting the vertical spread 
of pollutants ft>r different stability conditions and different 
aerodyaaalo rouglaneeses. Equations} using similaxily thsoryi have 
been developed to predict mean vertical spread} unlike nea»-iixtre®e 
apreeds given Hbnini Both these approaches are tested os flaid 
data# 

0 ott 8 ist«at with the actual alano spheric observations » the 

varlatioa ©f u and «„(*) with the helgat were also inoorpoiated 

y 

in the foraailatloa of e<Maoeatra 1 d.c» distribution in the downwiM 
direotioa froa a source# ^ 

Based on "he abovif a proeMure mb 'evolved to eon|}Ute 
poUutaat spread trtcing Into aeeouaat surface roughness of Ito© 
asrM as also the ateosphedo ttnrbulehoe mad stability as IniAeaM . 
^ various' aio»o««iteor®l©i^al parsttetwe# Besides wlhsnatleal 
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models^ tablso and obarts have been developed through numerical 
solution of the equation with the help of a digital computer for 
ready use on field* 

3*5 Bvolutlon of a Ctomprehensive Diffusion Model 
3*5*1 Batlmation of Vertical Spread 


In the derivations that follow, the following assumptions 
are madei 

(l) as soon a* the particles are released from a source, they 
instantaneously attain the velocity field of the surrounding 
environment I and 

(ii) the distance x travelled by the particles in the downwind 
direction in time t is the mean distance, x of ensemble of 
particles, each of which has travelled a specifio time. 


It is further assumed that the vertical displacoaent of Itie 
particles as in tlr^ t is the mean displacanent given by 
S (following lasquill, 1966 '). the above assumptions, 

i’asquill’s approach is applicable and the rate of the mean vertical 
displacement is given by Eq« (2.50a) as 


«> bu^ h G/li) 


Al 80 ,as shown Jasquill. 


d1 


/ 


dz / ; c dz 


0 .— ^ 


(3.18) 


(3.19) 


0 
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Unliki Idafi app3?08cli of ]tesquill (1966) vshers he assumes 
b « k » 0.40 for all the conditions of stability, it is proposed 
by the author -feat parameters of Eq. (3. 18) be .obtained from JBq. (3.19) 
separately for ©aoh case of stability, and employing these forms 
of dj^/dt, relationships between mean downwind distance x and mean 
vertical displacement z be established. 

In Eq. (3«19)» 0 is the mean concentration of pollutant 
at any height z in the vertical plane at a distance x from 
the source. lb got 0, one requires the shape of tbe distribution 
of pollutant in the vertical direction. As discussed earlier In 
Section 2.3. 2. 3 f the shape is assumed as an exponential form, given 
by an aquation exp (- bz®) where s varies from 1 to 1.5. Assuming 
a ■ 1, following Bosaaquet and learson (l936), the form becomes 
simply exp (- bz). Writing the equation in terms of for 
an exponential distribution, by definition of 

c »s exp (-1,414 (3.20) 

where c^ is the axial or the ground-level concentration, and 
0 is standard deviation of vertical spread. 

JZ 

lb establish the. relation between the standard deviation of 
a and mean dlsplaceKtsnt z of vertical spread, the first moment 
is taken (lapur and Saxena, 1972) as 
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z * 


Jo, 2 dz 
0 

m 

/ c dz 


/ c exp (-1,414 2 dz 


/ c exp( -1.414 — ) dz 
0 ° ""z 


1 


^ / 


1 


(l,414/o^r (1.414/0 ) 

z z 




1.414 


1.414 


(3.21) 


and o « 0 exp (- ) 

^ % 


(3.22) 


In the Eq, (3.19)> is Ihe eddy viscoalty and within the 
surface layer of the atmosphere, It is given by Manin-Obukhov (l954) 
similazlty theory as (Eq, (2.11 )) 




dw/dz 


(3.23) 


Eor ®Mh case of stability condition di^dz is known and can 
be derived and used in the Iq. (3.19)# 

While di^dt gives the mean vertical velocity of the ensemble 
of particles, the mean horizontal velocity of the ensemble of particles 
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has to be established. As it is assumed that the pollutant 
particles have the same velocity field as that of air around, 
the mean horizontal velocity at any height for any stability 
conditions can be derived from Eqs, (2.4a), (2.14) and (2,15). 

(Do get the mean horizontal velocity for the ensemble of particles 
in any plane, one can write 


m 

/ o.u dz 



UA W 

• y 

/ 0 dz 

0 

(3.24) 

Once 

di/dt and d^dt are established, 

the relation between 

X and z 

can be easily achieved by 



dx dx/dt 
dz d^dt 

(3.25) 


d^dt 

(3.26) 

within the appropriate boundary oonditions. 



Ihe oases of thermally neutral, unstable and stable atmospheres 
are studied separately in the fbUowing paragraphs. 

3. 5.1.1 Neutral Atmosphere ■ 

Itor the neutral atmosphere, the wind velocity profile is given 
by Eq. (2.4a) as 
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- 

u » ^ la ( %/z ^ ) 


du * 1 

**»*-»**• a mrnmm ^ mm, 

dz k z 

(3.27) 

Substituting Eq. (3.2?) in the Eq. (3.23) 


2 

u 

r 


M (u^z) 


Kjl “ 

(3.28) 

dKjj 

Honoo * ku^ 

(3.29) 


Using Bqs. (3.22) and (3.29) in the Eq. (3.19) 

m 

. / o^exp (- |) K (ku^) dz 

dz 0 

dt “S ^ 

/ exp(- |) dz 
0 


SimpliQring 


If " toi# (3.30) 


lb get the mean horizontal velocity of eneemble of particles, 


. / Q exp(- |) K ( ~ In VZq ) dz 

dx 0 

dt " “ 

/ c exp (- 1 ) dz 
0 z 
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/ oxp(-~)(ln z- In z )az 

IT* 

/ exp (- ~) az 
0 z 


u^, [- z (in y/ z) ~ z In Zq 3 
k 

z 

using Laplae Tranefonis (Etraelyi et al., 1954 )? Y being 

0 V 

luler*-Ms .80 heron t constant (« e , ^iiere c = 0.5772 ) numerically 

equal to 1.781. 


dx 
dt 


^ [-m Y+ in z/%U ^ i -0.5772 + m Vzq 3 




^ tin 0.56 + In z/z^] 


- u 
dx * 


^ t In (0.56 a/z^) 1 


(3.51) 


Fjom Eqs. (3-30) and (3.31), 


% * ~ [ln(0.56 z/z^)} 

dz 


•4 .ln(0.56 ^Zq) 
k 


X » 


/ In (0,56 VZq) ' 3-z 


ffi jfc 
"0 


\ [z In (0.56 z/z ) - 


k 


• $ 
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*g 15 In (0»56 i/z^) - z “ In 0,56 + z^ ] 


= -2 Cz (in 0.56 - l) + 1.58 z 1 

0 '0 

JRm 


kx m E i:in(0.56 i/z ) - T] + 1.58 


or 


O**# mm- ■ 

- f- [ 2.3 log (0.56 ^zj - 1 3 + 1.58 

0 "0 . ° 


( 3 . 32 ) 


ibr tho neutrally stable conditions of the atmo sphere » the 
Eg. (3.32) givoe the rehition between the mean vertical spread 
z and the mean downwind diatanoe x from the source, if the 
roughneae length is knovaa. 

feble 3*4 and I*ig» 3.5 show the and corresponding 

x/z^ values evaluated from Eq, ( 3 . 32 ). 

IIABIE 3.4 


!I?HE VmJBS OF z/z AND OOHRESIONDlia VAHJES OF x/z , 
0 0 ^' 

PROM EQ. ( 3 . 32 ), FOR NEUmAI. ADBOSPHERIC CONDITIONS 



x/z^ 



^ Zq 

i/^o 

10 

55 

70 

1180 

400 

11100 

15 

114 

80 

1405 

500 

14500 

20 

187 

90 

1650 

600 

T8000 

30 

350 

100 

1890 

700 

21800 

40 

536 

150 

3220 

800 

25550 

50 

740 

200 

4650 

.900 

29400 

60 


300 

7825 

1000 

33300 




'■■rrr 

fihUTRAL ATMOSPHERE 


i f. 1 4-(- 


... ... , ♦ 


t-1 < 

1-' f ■" 


+ ■ t t 1 t-' t-t. 


t" ■■■:■ -‘""t' T*' 


'I 


"I ♦ I 

* ^ ^ ♦ i ♦ 


4 . I-. . * I j 


T' -;t"--r:.T:1 I'l ■■ n 

f- f 4 I- I 4-,l i 4 , . 






• f ! t ‘ f 


4 ' r - 4 '■ -t f 


. 4 . . • 4 ■* 


1x1 0‘ 


t ' - f *■ • 


5x10 


icajuf- 3.5 PLOT OF z/Zq Vs x/zq.QS given by 

j£q,(3.32) for neutral atmospheric condition: 
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3*5*1*2 Unstable Atnoephere 


Ibr the ynstablu atnospher®, the wind vslooity profile is 

given by Kq« (2«15) when free convection conditions prevail. Such 

free convection conditions occur in the range of z/L = -O.03 to 

-0.08, as discussed in Section 2.2*3. 3* Denoting this value of 

z/L where free convection occurs as (z/d)^, and assuming that 

wind velocity within the rang© of z^/h to (Vd)* is given by 

IcjgarithRiie l«w (Eq« 2.4a), the wind velocity profile for the 

entire rang© of unstable atooaphere can be written for the two 

layers as follows, 
a, 

ibr - i»e. , the forced convection layer, 

u m In (3»33a) 

k 

ili ^ ^ 1 

. . da k a 


Prom Eq. (3«6) 2 

( «# 

Si 


(uyka) 


■ ku^s 


dK^ ' 
_.ku,. 


(3.34a) 


Bbr (- f < - # , i.®M the free convection layer, 



U - ~ E-1.3(-* + oottstant ] 


(3.33b) 


V lil f 2 \-4/3 / \“l 

di * F " 3 i ) 


Eroa Eq. (3-23) 


u 






- 2.31 kU|,(-«/l)*'^^(-L) 


d2 


2.31 eI- 


»• 3.08 fcu^ (-z/L)’’^^ (3.34b) 

Using Bq®. (3.22), (3.34a) and (3.34b) In -the Eq. (3.19) and 
writing * (-a/l»),, (-L), 


dg , 0, 
dt 


/ o^exp(- I ). ( ^).,dg 0 ^ exp(- f ) ( ^) 2 dz 


/* 0 ^ exp (- z/z) dz 


■^^jjere ( •g^)^ is tbe rate of change pf Ejj in the range of 
“ 1." z/li l.(- 2 /Ii)^, and 
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da 0 
dt " 


/* C^0xp(- '§)(ku^)dz + f C_exp(- •§) x (3.08 ku^(- f ))*'/^az 


m ~ ■ 

I 0 exp (- ■? ) dz 

0 z 


lcu*[--z exp(- "l)]^ + 3.08 [/ exp (- ■§ )dz] 


% 


1/3 


2 .-«* 


[-Z ©xp (- ~ )j 


dz 

dt 




ku^J1-tixp(- ~)3 + 3*08 exp (~ ■§ )<!! 

2i Z (“Ii) z# z 


(3.35) 


Uaints Ijiilmu TrtinBforiiis (Brd©lyi et. al., 1954), the integral 
on the right-hand aide oan be evaluated as • 

/ z^'^^exp (- a/z)dz « [z]^ ^ r (4/3, z^z) (3.36a) 
z* 

where r (4/3? z^z) is the inDomplete Gamma function, which is 
given by (J^rdolyi, et. al. , 1954) 


“ l/3 

r(4/3,z^/z) » /_ 03cp C-z) z^ dz 

z^z 

z/i 

B r(4/3) - / exp (-z) z'' ^ dz 
0 


tl- 


W3l 


/ exp (— z) z^^^ dz 3 r(4/3) 


[ 1 - E( 4/3 i z^z) I r ( 4/3 ) 


(3.36 b) 
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where r(4/3) is the complete Qemma function, and l’(4/3, z^/z) 
is the incomplete CSomma function ratio, as tabulated by 
(196S). Oubotitutifig Eqs. ( 3 . 36 a) and ( 3 . 36 b) in the Eq, ( 3 . 35 ) 

^ 1 B3cp(- - 3^)3 + 3.08 r(4/5 ) [l-l(4,~)] 

z z(-L)'^^ . z 

« ku^[l-exp i-^) + 3.08 X 0,8934 (- i/l)”''^^(l-i(4/3,Z4(/2))3 
z 

oinou value of r ( 4 /^) » 0.8934 

|| - ku^ Il-exp (- ~) + 2.75 (- Vl-)"’'^^(l-l(4/3, z^z))] 

25 

( 3 . 37 ) 

Bq, ( 3 . 37 ) shows the mean vertical velocity of the ensemble 

of particloa in the unstable atmosphere. It can be seen that 

z^/z s (- z/h)^ X (-e)/z ■ (2/l)^/{z/l')» and for small values 

% 4 % ' 

(o.g, 0,1 or less) of (z/li)y(z/li)) [l~exp (- ~)3 and P (j, ~) 

tend to become veiy small, and hence the Eq. (3*37) can approximately be 

written as 

If - ku^t2.73 (- 

or If « 2.75 (- for (- i^l) » (- z/l)* (3.37a) 

3b get the mean horizontal velocity of ensanble of particles 
in the unstable atmospheric conditions, the wind velocily profile 
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Eqs. (3/52a) nnd (3.32b) has to be fully establi^ed. At (- z/L)*, 
both u anri K|^| profiles have to be continuous# Bbr K^-profile 
to bo continuous at a » (- x («l) « 


Kjj » ku.^ z^ for the lower layer 


- 2.31 « talj, (-1) I6r the 

■ - 2.31 (- 

OT ( — i'/’. _!_ 

• '• -L^ 2.31 


upper layer 


or (- 7 “) « O.oa 

L 


(3.38) 


i#e*(~^If)^» 0.08 

Ihia shows that upto the heigSat " 0,08 x (-l) from the surface, 
the logarithnio vdnd profile will b© applicable and beyond Ihe 
l/3l?rofile will dominate. At z - 0.08 x (-l), the wind velocity 
given by each of the formulae must be the same, l.e., 

; . ^ in 

f#C 


or 


^t-1.3 ( # constant ] 

■ S* <. S§ .,?f Bs - 1 #3(0*08)*"'*^^+ constant 

O' 

oottstant «■ 1.3(0.08)"’^^^ + 2*3 3ng 

» 3. 08 + 2.3 log 0.08 - 2.3 log (zq/-I.) 
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“ 3.00 - 2.52 ~ 2.3 log (z^Al) 

» 0.48 - 2.3 log (zy-Ii) 

Rubotituting thu wlue of the oonatant, u, beyond z = 0.08 x (-i), 
bocouoe 

S = ~ [0.48 - 2.3 log (z^Ai,)-1.3(- ] 

BiuS| tho overall wind profile beoomeB 

for - A -f 

Eo. 48 “ 2.3 log (Zg/-li)- 1.3{“a/L) ^“^^Ifor (3.39) 

- z/l > (- z/L)* 


Using B(j. ( 3 . 39 ), the mean horizontal wind velocity of 
ensemble of particleB oan be now established 

c. ~ In dz + /o. -[0.48-2.3 log(z .3(- f )“'''^^]dz 
too A z* ^ ^ 

cit “ ^ 

/ 0 dz 
0 . 

(» 

E / ©xp(- '“).ln *1 dz + /exp(- -5) 
o z 0 z* i 

(0.48 - 2.3 log (zqAi<)-1 .3(z/-^r'‘^)dzl 

m ■ ' 

/ / exp (- •§) dz 
0 z 




77 


u„ . # 

* 1 I / 


I / «xp(- )ln z dz - / exp (•« - ) In 


dz 

0 


+ / (0.4a - 2,3 log (z /-I,)) Qxp(-z/i)dj 


^ , yu 

/ (i.3C2/-Lr'’/^) exp (-z/z)dz] 


^ 1 i J axp (- ~) In zdz - / exp (- - ) In z dz 


z 0 


% 2 


In z (“i oxp(-~)) l % (0.48 - 2.3 log (- ~))(-.z exp(- -r))!’ 


sxp(- z/z)dz 3 


Using Laplao© Tranaftorae (Brdelyl, et. al. , 1954) to evalmte 
th© thr©« integrals in above, 


(i) / exp (- •§ )ln z dz a ~ 2 In (Vz), where y is 

0 z 

ajlar-Maseheroni constant (» e°, #16x6 c = 0 . 5772 ), 

r “ % - 

(ii) I exp (- 3) In z dz » z [0xp(~~ )ln z^- E. (- z^z) 3 


iW> f ■•w’l " ' 

where - E.(- z^z) ■ / „ exp (~z ) z dz 


which is the exponential integral, tabulated values being given by 
PagurovB (1961). . 
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(iii) / z oxp (-z/z)dz *= (2/3, zJ'z) 

whory a^nla r ( 2 / 3 > z^z) Is th® iacomplete Gejama ftuactioGj ani^ 
is given by 

CO 

r (2/^, z^/z) « / exp dz 

z/z 


z^i 

«r( 2 / 3 )- / exp (-z) 2”^/^ dz 


Zy^z 

» r(2/3) [1 1 / exp(-z) z“"'^^dz ] 

r( 2 / 3 ) 0 


- r(V3) [1 - I (2/3 , z^z) ] 

where "^{2/3$ z^/'z) is the incomplete Gfamma function ratio, as 
tabulated by Khamis (1965). Substituting these three integrals, 


die 

dt 


— [-Z In (Y/z) - 2 (exp(- z^/z) In z* 
k#z 


- - % 

“ E^(“ z^z)) - In z^iz - z exp (- — )) 

z 

% ' 

+ (0,48 - 2.3 log ( •—)) (0 + 2 exp (- ~ )) 

■Z: 

a)^ r(V3) (1- 1(8/3, VS)) 1 


Simplij^axig the equation, it becomes 
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-0.577a - (->•). to 

0 Z' 

+ + (0*48 - 2.3 log (“i^)) exp (- z^z) 

- 1.76 (i-I(2/3,Vz))] (3*40) 

Bie Eq. (3*40) gives the mem horisDErt^l velocity of -tbe 
enssQbl® of particles. On the evaluation of .dS^at, it was 
observed that fer saall values of z^/i (i.e. (-2/li)^/(i/li))| the 
8UJ.iinatlx3n of first four terms on the right-hand aids becones almost 
negligible for different values of z^. Similarly exp(-. z^z) '-^1 
and l’(8/3fZ^/z) *►0 for this condition, hence 

II « ^ £0.48 - 2.3 log (- i) - 1.76 (- ] 

for (- -I ) » (- ) (3.40a) 

Onoo the mean vertical add horizontal velocities have been 
arrived at, the relation between the downwind distance x and the 
vertical spread z can be established* From Eqs. (3.37) and (3.40), 


dx dx y dz 
“ dt ^ dt 


% « £-0.5772 - In ^ exp (- ) + In z/z, 


dz k 


. E,(-V^) . (0-48 - 2.3 1»8 (- r» 



Here 


- 1.76 [i^xp (--^ ) 

■ ' z '■■■ 

+ 2»7§ (S/~li)*''^(l«I‘(4/3 > 2^/21 ))3 (3»4l) 

% “ (aZ-L)^ x (-1) e 0.08 X (-1), hence 


2 - 
k OX 

dS 


C-0.5772 - 2.3 log— 2:^ (, .0^^) ^ 

(2 q /-I<) (2/-1) 


2.3 log ( ) - (- E,(- 0.0V(i/-l))) 

Wq / 4« 


(0.48 - 2.3 log (-z/l)) exp (- O.OV(V-l)) - 
1.76 (V-Iir’''^(l-i(2/3, 0.08/(i/-I,)))]/ 

[1 - exp (- ~2lS§» . ) 4, 2.75 (i/-I')^^(l“P( 4 » ))] 


(iZ-l) 


(z/-I.) 


Integratlngy 


k^ X / 2 dx , 

•pjj *; ,, , lie ~ 3d(z/-i,) 


Zq /C- Ii ) dz 


/ ,[ -0,5772-2.3 log 

z /(-Ii) (zV-l) 


0.08 . / 0.08 


exp (- 


(z/- 2 j) 


) 


“/ T 

2.3 log ( f^) " (-Ej^C-O.OS/C^-I,))) + (0.48 - 2.3 log(— )) 

exp (-0,08/(^-I/))-1.76(i^-I.)"’^(l-P(2/3,0.08/(^-L))) 1 / 

. . , 93.0 0 , )^ 2. 75( ( 1 -P( 4 / 3 y - - 2 ! Og . , ))] d( 2 /-l) 

(^-1) (z/-I.) 


(3.42) 
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The Integral of Bq. (3.42) was evaluated numerically by using 
Gauao-Logendro Quadrature method (Oarnahan et. al., 1969), for 
which number of points used was four, in an interval of VC-l) = 0»01» 
The ooraputor p.rograi.Mu for Ihe evaluation of this equation numerically 
is shown in the Appendix A*1. Biree cases of different values of 
z^/(-L) were studied. It was observed that for z/(-l))_> i.o, the 
integral of Eq. (3*42) could be almplifiad, as many terms in the 
nujfiurator as well as the denominator become negligible. Using 
Eqs, (3.37a) and (3.40a), a relation can b® established as 


ilE « 1- [0.48 ~ 2.3 log (- z /I) - 1.76 (- / 

dz k*' ° ‘ 

C2.75 3 

»• W £(0.48 - 2.3 log ( -7)) ( -2-)“''/^ 

2.73k‘^ 

- 1*76 

• k^3^(-Ii) » 1(0.48 - 2.3 log ( t^)) C ■^) 


-1.76 ( |=.)“2/^3 d(;/(-i)) + ( ^)i 


wl»ro C '^^)^ is the value of upto (z/-l) = 1, 
arrived at by numerical integration 

. 2 - 


[(0.48 - 2.5 log ( ^))( 


2/-I. 


, 2 - 

. / 

"T ' \ y .j 



Eq. (3.43) givus the relationship between the mean downwind distanc 
X and the mean vertical spread z when i/(-rl) >1. 5Chus solving 
Eq* (3.42) nuraoricnlly the values of can be derived Ibr 

S/(-I*) 11. Sbr ^(-I) >.1, the Eq. (3.43) can be used. 

In thu above derivations, it is assumed that free convection 
oecura at ®/(-L) » 0.08, When the ratio of roughness length to 
stability length, is equal to or greater than 0,08, the 

free convection oonlitions can be aseuned from 2 q/(-e) value of 
a/(-I/) itself. In such oircuaetanoes, Eq, (3.39) giving u has 
to be revised. In the free convection conditions, 

u » ^ {-1.3 (V'C-l))*'^'*^^ + constant ] 

Whm si/C-Ii) ■» z /(-h) ,u«siO 

W‘ ^ 

0" M*5 + constant] 

t # K y 

mmi 

or, constant « 1.3 (zq/(’‘I')) 

Hence, u « [1.3 (z^A-E)) 1»3(z/(-If)) 3 (3.44) 
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/Old, K. 


.31 k u^ (‘~5r)'*^‘^^(-L)''^ 


«’ dK- 
/ 0 J ds5 

dia 0 

» * w* » — > ‘ 1 -HI-- 

at . 

fo dss 


/ f) »< (3.08 ku^ K ( 


/ 0 .xp(~l)dz 
0 i 


f V5 / JZ5 \ 

3.08 ku. / ® (- —) <lz 

o . ^ 


/ exp (- •^) dz 


3.08 ku 


z(-L)^ 


i irUA) . z 


3.08 { 0.8934 ku^ ( 


If r, 2,73 ku. 


(3.45) 




- /«•' « «lss 

dx 0 

dt * 7 

/ c dz 

0 
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1 % 

/ exp(- •5) d .2 




.1...*.?.^^ r ^«£n“1/3 “ .G) 

;|^ \ Y ^ "*" ■WWiWWiPPiPBP 


^3 


k ig 






K- 


4)' 


•1/3 


1.35 (■ 


£-)• 

-Ii'^ 


•1/3. 


(3.46) 


Eqa. (3.4‘i) and (’3.46), the relation between x and z can 
be achieved, an 

to ^ dx/dt 
da da/dt 


li2^ j(^)“1/3^J)-1/3^ 


2.75k 


• . ily - 1;% t(ir^^( J)-'/5- 1 .35 {!5 )-^^,.(-^-l ) 


0.473 I - 1.35 X 3 ( 

a/(-L) 


0.71 I(:4^’**^^^('!t)^- 2.7 ( 1.7 ( 


k X 

T~£l' ^ •■'-L' ^ -I' ■ ■” ' 


(3.47) 
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Kq. (-5*47) rolutuB % and i for the case (z/-!.) > 0.08. 

hqji, (5«4lO •" (3»4!5) and (5.47) the following EDable 
3.5 ia pnujtttriHl, wixu'o A/(-l) values fbr different values of 
z/-I» are tabulated for the oases of z /(- l) = (i) 0.0001 
(ii) 0.001, (ili) 0.01 aM (iv) 0*1. Idgure 3*6 graphio ally 
presents these values# 

Itm yqiwititnio (3.4«), (3.43) and (3.47) or the lable 3*5 and 
Hg. 3.6 c»m he uned to predict the vertical spread fbr different 
condition of unstable atmosphere. 

3.5*1»3 Otable A tro sphere 

As disouoaed in Seotion 2. 2. 3*3 i the wind velooily p3x>file 
under atable atmosphtrio ooMltiona can be given by Bq,. (2.14), viz., 



{In 




0 


(3.48) 


where « io the universal constant, the numerical value of which 
can be taken m 6.0. h'or all practical purposes, the last term 
of the equation (3.48) am be approximated as and the whole 

equation becaoee 

0 ■ ■ ; 

The first derivative w.r.t. z will be 



TABLE 3.5 




™ CORRESPONDING VllDES OP 

T.OR DIBmMT«7(-I,), BY USING Eqs, CJ.aT) ! (j.Jjf 

(5.47)* WH nWfJTAJm*: A'MJSHIERIC CONDITIONS 



0,02 

o.osoe 

0,03 

O.C^78 

0,C4 

0.1334 

0.03 

0. 1674 

0,06 

0,20l0 

o,or 

0.2349 

0,08 

0.2695 

0.09 

0.5038 

0,10 

0,3382 

0,20 

O.6B62 

0,50 

1.C4Q2 

0,40 

1.3822 

0,50 

1,7142 

0.60 

2.0542 

0.70 

2,5442 

O.AO 

2.6442 

0.90 

2.9392 

1,00 

3.2267 

2.00 

5.8570 

5.00 

12.0570 

10,00 

20,3570 

20.00 

55.8570 

50.00 

64.8570 

100.00 

105.6570 


0,0222 

« ' , 


0.0397 


m 

0.057 1 

mm- 

m- 

0.0742 

0.0002 

m 

0,0914 

0,0010 

wm 

0,1094 

0.0031 

m 

0,1281 

0,0062 

, : 

0,1475 

0,0103 

.'m 

0,1770 

O.O151 

pm 

0.3690 

0,0714 

m 

0,5960 

0.1669 

- 

0,8230 

0.2779 

- 

I.Q470 

0.3949 

0.0040 

1,2660 

0.5129 

P.P 42 P 

1,4800 

0.6319 

0,0650 

1.6900 

O.75C4 

0.1000 

1,8970 

0,8609 

0.1580 

2,1000 

0.9866 

0.1710 

5.9800 

2.1400 

0.576 0 

8. 5100 

5.000 

1.7550 

14.6000 

0.9400 

3.5300 

24.5000 

154600 

6,6600 

47.6000 

30.8000 

14.2100 

7S,(XX)0 

51.1500 

24.6600 
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(3.49) 


dU r 1 * 1 

da k ^ a Ii"^ 


% 

k» 




Bi® ®ddy visooslty Kj^ is given by 


and 


u„ 




du/&z 

2 

u 


'li 

la 


( 1 +• <« f/L ) 


K 


ku^ a 


M 


( 1 +-S/L) 
dK^ ku ^ ( 1 +« j'/L ) - a{«* /l ) 


ku. 


da 


(H«a/ri) 


Uoittg ttiio wqutttlon in Bq. (3*19) 

dKj^ 


da ^ o 
dt 


/o * 7 -* da 


da 


/ 0 da 


w» . jj ^ ku^ 


/ “q »«P (-*“) 


a (l •»« 2 /L) 


r 


da 


lid 9'' \ 

/ 0 «xp (“ •£ } da 


(3.50) 


(3.51 ) 



1 


/Wp (-^) 


ku. 


» 


dz 


I oxp (~ 4) dz 


1 


lcu.1. 


/ o 

2 S (1-^g/L) 

t-S exp C- f )3“’ 


dz 


«■ »** / Z V 1 

/ exp (- -) 

0 z (l+“r*; 


dz 


(3.52) 


Subotltutlnij « T, dz *» •* dT in Sq. (3«52) 


ku^ “ 


d t 


(.IS) 1 . 


Z 0 


«z 


(1+f)‘ 


dl 


letting 


P» 


e*?. 


dz 

dt 


## i 

m ku^ / exp (-pT) ^ dT 

0 (1+f) 


Using teplaoe Tranaforms (Srdelyi et. al., 1954) 


|| « ku*p (1 - p exp (p) (- E^(-p))3 
where [(-K,(-p)l «» / exp (-l) T dl, vsfaioh is the exponential 

1 P 

integml. Ihe values of this integral have been tabulated by 
iagurova (1961 )• 

■“ -i « ^ -=■ 


(3.53) 



PJq. (3.W) gives th® mean vertical velocity of the ensemble 
of particles In the stable atmospheric conditions, to get the mean 
horis^tontttl velocity, Bq. (3.24) gives 


/c u dz 
^ 0 

f c da 

0 

(“ «) rKltt i + ■®|' }] dz 


/ 0 exp (- ■£) dz 

o ® z 

^ I Axp (- •!) la ^ dz . /“exp (- f ) ^dz] 

^ /“«xpC-a<iz 

0 z 


£ « — r / 0 xp (- •*) In z dz - La z / exp (- •§) dz 
kz o z 0 z 

W 

+ ~ / z exp (- •§) dz3 (3»54) 
■^0 z 

llMing Laplace Transforms (Rpdelyl et. al., 1954) 


/ mp (- ■§) La z dz « - z In (v/z) (3«55) 


where v is ELler-Maeoheroni constant, numerically equal to 1.781 



Substituting Bqs. (3.b5) and (3.56) in Eq. ( 3 . 54 ), 


•mm 




m 

kg 


I- i In ( v/«) - In z xS 5 2. 

0 Jt J 


“1^ lnv+ In z - In 2 +72) 

0 Ji4' 


% 

'1^"” { 1® { 0.56 ) + in z “• In 2 ^ + z ) 


. “* u 
dx * 

dt " k 


{In ( 


0.56 i v 

m * 



From Bqii. ( 3 33 ) tmd 3.57)i the relation between 
z is derived ma 


dx di^dt 

•MM im m mm Jhm m mm 

dS d^dt 

^{2.3 log ( ) +<*|] 

lu j "(«. exp ( "^) (- ■^))] 

«i «2 «« a 2! 

^2- f/l {2.3 log ( +«|)] 

^ % { 1 --^ exp (i) ("V“ 

Wg «Z ■; «Z «Z 

Again %, (‘j.bS) was solved nnaerioally, talcing ce 1 
different valuefi of Si® iategiml was evaluate^ by 

Ijegindrt Quadrature method (Oamahan ©t.al.f 1969). Bbr 


(3.57) 

X and 


(3.58) 

I 6, and 
Gauss- 



z/l “ (i) 0.0001, (li) 0.001, (iii) 0.01 and (iv) 0.1, Eg. (3.58) 
wmo evaluated and valuta were tabulated (!I!abl® 3.6) and plotted 
(Pig. 3.7). 

'The aquation (3.58) or the toble 3.6 and Pig. 3.? can be used 
to estimte the vertical spread different conditions of stable 
ateoephere. 

•Ehe vertinal spread of pollutants can be predicted by these 
©quttttonw (l*:qe. (3.32), (3.42) - (3.43) - (3.47) and (3.58)) or 
Triblee 3.4 through 3.6 imd Figures 3,5 through 3.7. for neutral, 
unatabl# or stable atoioopherio conditions. Mae validity of the 
prediction by this af^rowh is tested in the following chapter for 
different aerodyniimlo roughnesses of surface and different stabllitsr 
oonditione of atmosphert. 

It sitioudl be eaphasised that in the above approach (Bqs. 
(3.32), (3.42) - (3.43) - (3.4?) and (3.58)) the vertical spread 
eotimtod is the aean of tiie eneeoble of particles. 

3.5.2. Estimntion of lateral Spread 

M stated in Section 2,3. 2. 2, In the Sasquill-Gifford method 
of eotiaating coTOtntMtIon distribution, a set of ouives is used 

in which a wid o are given fer downwind distance frcm a source. 

■ ■ y , I- , ■ 

% estimate rasquill (1961 ), in his original paper, gave 

y 

(reported earlier m Etj. (2.42)), 




(3.59) 


TABia 5,6 

m. VA11JF.S OF z/h m CORRESOTDINO TlXtJlS OF FOE 

I);LFt'T?Ri®T K A, USTOO Eq, (5,58), FOE STiBHl IMOSHEEIG 

coffiinoNs 
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or Oyiz) « xR (3.64b) 

Xft Eq[# (3*64^1) I ^Q(z,j) find 11(0^.) sr 6 usually Di 8 Esux* 6 d Et 8j3y 
height * If <70 is aot availablei but temperatures at two 
heights and wiM velocity at one of these heights are available, 
using Jaasfsky and Iraaada (1965) s» method, Sec. 2. 2.4, 2, 0 ^ 
can be estimated, ishioh can be used in Eq. (3.64b). 

Oie effects of roughness length and stability (given by 
li) are taken care of in determination of o. (or a ), and hence 
lq» C3*64a) or ( 3 *646) can be used to predict standard deviation of 
latexwil spread of pollutant at any distance x from a source. 

3*5»3 Estimation of Gbnoentration Distribution ^ 

Once a (x) and 0 (x) (or z(x)) are determined, the 

y ^ 

concentration distribution of pollutant can be estimated. Unlike 
earlier efforts, it should be noted that in the present development 
of estimation of 0 (x) and 0 (x), the concept of inhomogeneous 
field of turbulence is not saorified. 

Also unlike Easquill-Gifford method of estimating concentration 
dlsttibution using o^Cx) and cr^(x) , where they assume u to he 
invariant with hil^t, Sec. 2.3*2.3, above formulation allows for 
variations of u with height. As discussed in that section, 

oonoentration distribution follows Gaussian shape in the lateral 
direction. Assuming simple exponential variation of concentration 
in the vertical direction, as done in development of a^(x), w 
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the mean concentration of Bonu+o«+ =4. ^ . / x 

pollutant at aoy pojjit (x,z) ±n the vertical 

plane passing through axis can be «riven hv b , , I 

Q- vwi oe given by (see Sec. 3.4.1,Ilq. ( 3 . 4 ))) 

0 (x,z) » 0 (x,0) exp (.1.414 ' (3.65) 

Stauarljr, meai oono.ntration at any potot (x,y) la any hortaontal 
plane can be given by Gausian distribution as 


^ (x,y) » c (x,0) exp (. ) 

2 2 


(3.66) 


Combining Eqs. (3.65) and (3.66), the mean concentration of 
pollutant at any point (xjyiig) can be given by 

o(x,y,z) » c (x,0,0) exp (--I 2-) exp(. 1.414 -^) 

& rffc CJ 


or 


y 

M 2 

c(x,y,a) . 0 (x,0,0) exp C- | 1.414.2-] (3.67) 

c 

y 


Sbr a continuous point source at ground level, the continuity 
equation gives the Q, quantity of pollutant released per unit time, by 


IB 

Q » / / u c ( x,y,2) dy dz 

0 -£» 


(3.68) 


Ihxim Eq, (3i67) 


m m 


/ / U c (x,0,0 exp[-^ ^ - 1.414 ~ ]dy dz 


0 ■ 


m ■ . ' » '2 

c(x,0,0) / u exp(-1.414 — ) [/ exp (- ~ ^)dy] dz 

0 02 ^ CT 


o(x,0,0) / u exp (-1.414 )[/5V ] dz 
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o(x,0,0) / u a exp (-1.414 ~) 


dz 


From Bq. (3.64b), cr ■ xR Oy/u 


Q * o(x^0,0) /& / u ( -^) exp (-1.414 ~) dz 

o u z 

As 0 is liivariant with hei^t, 


W 

Q « c(x,0,0) /2it xR a„ / exp (-1.414 ~) dz 

T a 




Q - 0.707 0 (x,0,0. ) xR 0^ 


or c(x,0,0) 


0.707 xHc^ Oj 


(3.69) 


Substitutlog Eq. (3,69) ia Eq. (3.67) 


o(x,y,z) 


0.707 *^ir xR a„ a„ 

V SH 


exp I - 4 ^ “ 1 .414 * i 


2 ^2 


W ' 

KnowUag that 0.707 «. - z (Eq. 3.4), and xRa^ * u, 


m. . 

o(x,y,z) 


/2tr uo^ z 

¥ 


r 1 jL- 

exp I- 2 **x - r j 

* 0 Z 


(3.70) 


In Eq, (3.70), a and u are to be taken at liie same height. 
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3*6* Oloauie 

Thus, consistent tdth the real atnioaih.eiic inhomogeneous 
turhulenoe field, the spreads of pollutant and hence its 
oonoentration at any point doTOwind can he estimated from the 
formulalion offered in the preceding paragraphs. For any 
diffusion study, iaae surface rou^ess length of the site 
should ho determined first. The stability length L at the time 
of diffusion ©xperinent should be estimated. Once z^ and L are 
known, the spreads can be decided and then concentrations can 
bo calculated* 

This formulation and the extension of Monin's approach 
h«v© to bo verified on the field data, Suoh a verification on 
tho scsrailahl© field data for two widely different surface 
iTO^ess lengths is roprted in the following chapter. 



Chapter Four 


VERIFICATION OF IHS MBMAIPIGAL MDEBl ON AVAILABIE FHID DATA 

The fflathentatic al model developed in Chapter Three Twould he 
meaningful only after it has been validated hy testing on field 
data* The estimation of vertical spreeds, and hence concentration 
distxthution at any downwind distance from the theoretical consider- 
ations disouseed earlier has to he verified with actually observed 
spreads and eonoentrations. This veiiftoation must he done under 
different conditions of surface roughness and atmospherio stability. 
Collection of woh oomprehensive data would be very expensive 
and time consuming job. Fortunately there is a significant amount 
of data oolleotad by various research workers that can be used for 
validating the model developed. The diffusion programme conducted 
at the open, unobstructed terarain of O'Neill, Nebraska, TJ.S.A,, was 
one of the most relevant studies of this type. For urban areas, the 
dispersion study conducted at St, Louis, Missouri, U.S.A. is 
similarly a oomprehensive stidy and has been selected for use to 
verify the proposed mathematical model for areas with large surface 
roughness, 

4,1, Disperaion Study at Project Prairie Grass 

During the summer of 1956, in United States, the Geophysics 
Reseaojoh Direotorate of the Air Force Cambridge Research Center 
sponsored and directed an experimental pirogramme in miorometeorology. 
The field alt© was carefully selected to provide level terrain. 
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uniform vegetation cover (maintained at an average height of about 
6 offi hy mowing)} and prao tic ally complete freedom from large roughness 
©laments. Ihe site selected was flat, prairie country near the town 
of O’Neill in north central Nebraska, U.S.A. Oixe programme was 
given the name, Project Prairie Grass (or, in short, PBJ). 

®10 primary objective of the programme was to determine the 
rate of diffiision of a gas, emitted continuously at a point, as a 
function of meteorological parameters. Ilhe experiments were 
conducted in a wide variety of weather conditions,, viz, wind speed, 
olcud cover, therms^ stratifioatton, etc. 

Experimental Techniques 

The tracer used in Project Prairie Grass was sulphur-dioxide. 
The diffusion experiments involved determinations of average 
concentration at selected points dcwnwind from a point source of 
sulihur-dioxide gas operated continuously for periods of 10 min. 

The tracer gas was released horizontally at a height of 46 cm above 
the grcnind during 64 experiments and at a hei^t of 1,5 m for six 
experiments. 

Midget! implngers wei© used to sample the gas. The implngers 
wore mounted at a height of 1.5 m along 5 oonoentric semicircular 
a»8 located at travel distances of 50» 100* 200, 4 00 and 800 m from 
the source. An angular separation of 2® was used between Impingers at 
the 4 inner arcs and a separation of t® was used at 800 m. 
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Wind apo®d and temperature measurements were made at heigjats 
of 25 t 5^1 100 » 200 , 400, 800 and 16OO cm. Wind speeds were 
measured wi 14 i modified Rikoken anemometers whicli had heen calibrated 
and matohed. air temperature measurements were made with 
radiation-shielded oopper-oonstantan thermocouples. 

Of the 70 gee releases during July >56 and August » 56, there 
were 34 during lapse or unstable conditions and 36 during 
invorsions or stable oonditions. 

4.1,2. Surf SCO Boughnese of PPQ Site 

Analysing the vortioal wind profiles at the PPG site during 
noar-neutrel atoosphorio oonditions, the roughness length 2^ for 
the site have been suggested in the range of 0.6 to 1.0 cm. 

(Hffligen at al., 19^1 1 Germak, 1963} lumley and Panofsky, 1964). and 
tti gwerag© of 0,8 om is adopted here. 

4 , 1 , 5 * Stability lengths for PPG Experimaats 

'Phe stability lengths L for ihe PPG diffusion experimaits 
have b®aa mggosted by two different workers. Takeuchi (1961) 
analysed the wind data and reported L values for all the 70 
experimental runs* Haugen et al. (l 96 l» reported by Pasquill, 

1966) oaloulttted Eiohari son number Bi (Eg., (2,7)) and thence 
suggeoted L’ values, whej® 1’ » L is the eddy oonductivi-ter 

and the eddy visoosity. The ratio of I^/i^ varies wilh stability 
(lumley and Panofsky, 1964), ^ is not unigueiy determined. Hence, 
instead of using L« and Kg/i^ to derive h, the values of L 
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oaloulated directly from wind-temperature data by Tafceuchi have 
been adopted and reported in Table 4.1. 

4.14* Vertioal Spreads at lOO m from Source 

ina^ysing the diffusion data, Haugen et al# (1961) have 

deduced the values of the standaid deviations of vertical spreads 

at 100 m from the source for all the useful runs of the piogramme, 

Oaussian distribution of concentration in the vertical was 

assumed in deriving the standaid deviations 0 . These cr„ values 

z z 

are reported in Table 4,1. Out of the total 70 runs, 48 runs were 
found to givu statistically reliable results (Hswgen et al., 1961), 
Hence 0^^ valuee of those 48 useful runs only have been reported in 
the Table 4«1* 

Out of 48 useful runs 19 were during unstable atmospherio 
conditions, 11 during near-neutral and 18 during stable atmospheric 
oonditiona. Data of 0^, said corresponding L are reported for these 
■three classes of atmospherio stability. As already discussed in 
Ohapter 2 the stability length I* is negative in the unstable 
atmosEhore, infinite at neutral a'boosphere and po si live in the 
atable atmoaihero. In part i of Table 4*1, entiles are made for 
unstable atrooophorio conditions in the seijueiKse of decreasing 
instability. Part B reoo ids near-neutral cases, followed by 
Part C in which entries aie in order of increasing stability. 
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TiBIE 4.1 

STmm DMIATIONS OP VERTICAL SPREiD AT lOO m m 

STJMan LiaTGms l for ppg bxperimhits 

(A^ UHSTABLB ATMOSffiERIO OOimiTIORS ^ 


SR 

No. 

Ran 

No. 

0 , 

m 

li'l 

m 

0 

1 

SR 

No. 

Ran 

No. 

m 

L, 

m 

0 

X 

1 

16 

10,1 

-5.9 

-16.95 

11 

49 

4.5 

-35.0 

-2.86 

2 

15 

7.7 


-12.82 

12 

50 

5.9 

-55.0 

-2.86 

3 

52 

10,0 

-10,0 

-10,00 

13 

61 

4 .8 

42.0 

- 2.58 

4 

25 

9.0 

-11.0 

- 9,09 

14 

62 

6 .9 

“44*0 

-2.27 

5 

43 

7.0 

■-21.0 

- 4.76 

15 

30 

5.6 

-51.0 

-1,96 

6 

19 

6.0 

-27.0 

“ 3.70 

16 

20 

5.8 

-65.0 

-1.54 

7 

51 

4.8 

-50,0 

- 3.33 

17 

31 

6.5 

-74.0 

-1.55 

8 

27 

5.8 

-31.0 

- 3.23 

18 

35 

5.1 

-89.0 

-1.12 

9 

44 

6,5 

-51.0 

- 5.23 

19 

48 

5.8 

“69 .0 

-1.12 

10 

26 

5.7 

-53.0 

- 3.03 






(B) 

NE/sR-NlUmEi 

/jMDSHIBRIC CONDITIONS 





20 

54 

5.3 

-100,0 

— 1,00 

26 

55 

4.4 

200,0 

0.50 

21 

45 

4.5 

-1C4,0 

- 0,96 

27 

42 

4.4 

160,0 

0.63 

22 

57 

4.9 

-250,0 

— 0,44 

28 

56 

4.4 

150.0 

0.77 

25 

24 

4*6 

580.0 

0,26 

29 

46 

3.5 

116,0 

0.86 

24 

21 

4.2 

250.0 

0,40 

30 

38 

5.8 

114.0 

0,88 

25 

22 

4 ■•^6 

240,0 

042 








TiBiB 4.1 (CONTD.) 
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SR 

Ho. 

Run 

No. 

Ogjf 

SI 

L , 

m 

m-'^xlO"^ 

SR 

No. 

Run 

Ho. 

02 , 

m 

L , 

m 

1 

L ’ 

X 10"^ 

51 

67 

3.5 

88.0 

1.14 

41 

68 

1.9 

27.0 

3.70 

52 

17 

5.7 

88,0 

1,14 

42 

66 

2.1 

21.0 

4.76 

55 

35-S 

4 .0 

85.0 

1.18 

43 

59 

2.0 

16.0 

6,25 

54 

37 

4.5 

83.0 

1.20 

44 

58 

1.5 

12,0 

8.33, 

35 

41 

4.0 

80,0 

1.25 

45 

36 

1.9 

12.0 

8.53 

36 

65 

3.4 

60.0 

1.67 

46 

32 

1.6 

12.0 

8.33 

37 

54 

3.9 

58.0 

1,72 

47 

14 

1.4 

11.0 

9.09 

38 

60 

3.9 

55.0 

1,82 

48 

40 

54 

5.2 

19.23 

59 . 

29 

3.4 

29.0 

,345 



- 



40 

18 

2.7 

29.0 

345 
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4.1.5* IfStiogj; Spread from Proposed Poinmlation 

Knowing and L for the runs, the vertical spreeds at 100 m 
from the oouroo wore estimated from the proposed mathematical 
formulation given in Chapter Three. 

For the unstable atmospheric conditions, firstly values of 
zyC-L) wore calculated, followed by computation for the value of 
k ^(-L) where k is the Von Kaman’s constant and is equal to 
0.4, X ■ lOO m. Prom Fig« 3.6, z/(“L) values were read for 
each pjiir of z^/(-L) and k^3^(-L). There values are given in 
Table 4*2* Multiplying z/(-L) values with corresponding (-l) values, 
tho meiBT. vortical spreads a at 100 m were evaluated. To compare 
tiiese z values with the standard deviations of vertical spread 
observed, z should bo converted to o„* As the standard deviations 
of vertioal spread given in Table 4.1 are based on the Gaussian 
distribution, tho standard relation between z and 0^ of Gaussian 
distrlTmtion can be used as (Kapur and Sarena, I972) 

■S— m 0,80 (4.1) 

0 _ 

z ■ ■ 

Thus, for ©aoh run, values wore computed. These computed values 
of 0JJ, should compare with tho observed values of as reported 

in Table 4.1, Tho ratio of computed *0 observed for each 
run has been shown in the last column of Table 4*2* 

The case of near-neutral atnospherio conditions was taken 
ae ooaparhble to neutral oonditions given in Sec, 3.5.1.1. 
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X m 100 m, and « 0.006 m, was calculated and from Pig, 5,5, 

, I 

corroaponding value of z/z^ was read. Prom z/z^, z and th^' 
using Eq. (4.I) wore oalculated. The ratio of computed 0^ to 
observed wan calculated for each run. Part B of Table 4.2 
reports the reailts. ' ■ 

The oaso of stable atmospheric conditions was handled similarly 

and roportod as Part C of Table 4,2. Prom calculated values z /l 

0 * 

and k for each run, corresponding z/l values were established 
using Pig. 5«7* Jfeiltiplying z/h with b and using Bq. (4.1) the 
standard deviations of vertical spread under stable atmospheric 
conditioaa were ooaputed. Hero again, the last column of Table 4.2 
gives the ratio of computed a» and observed 

M 21 . 

Beaidea comparing the agreement between observed data with the 
approach proposed by the author, the validity of the extended Monin's 
approach to estimate the vertical spread as given in Chapter Three 
has also been checked. For this the standard deviations of vertical 
spread were computed using that approach, and compared with observed 
values of Oj,* In this oaso, for unstable atmospheric conditions values 
of bkj^(-L) were read from 

Pig. 5.2. The numerical value of b was taken as 1.2 (as suggested 
by Pasquill, 1966), Here again, k is Yon Karmaa's constant and is 
equal to 04, Prom evaluated. The value 

corresponds to the height at which concentration is 1C^ of axial 
(or ground -level, for ground -level sources) ooncentration* for 
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Gaussian distribution, the \ax relates to by (Pasquill, 
1962a) 

Vj/"® ■ ^•■'5 (4.2) 

Thus, for each value of corresponding values of were 

evaluated and compared wite observed values of a^. Table 4.3, 
part A, reports these results. 

The near-^ieutral ateo spheric conditions were treated 

separately in Part B of Table 4*3» E*or known value of ^21^, 

IS ,, ®*iS» 3*4* ®is ratio z „ /z was used 

max o ” max 0 

to evaluate and hence o^* Comparison of computed and 
observed Ojj is again made in the last column. 

Part G, Table 4*3* reports results for stable atoiospheric 
oonditions. For this cose again, off for each 

pair of z /L ond bk^L from Pig. 3*5* aud values cr were computed 
and compared wite observed values of o^,. 
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Ti*BLE4.2 

COMPimSON OP B^mm devi aeons of verecal spread as calculaed 

FROM PROPOSER FOHMEfMEON 


(a) 

ummm AWOsmERic 

CCHDIEONS 




SR 

No* 

Run 

No. 

“o 

Fi7 

k^x 

Pl) 

Z 

Computed 

m 

computed a 

observed a 

1 

l6 

0.00149 

2.710 

1.35 

9.95 

0.985 

2 

15 

O.OO1O3 

2.050 

1.00 

9.75 

1.266 

5 

52 

O.OOOflO 

1.600 

0,70 

8.75 

0.875 

4 

25 

0,00073 

1.455 

0.62 

8,52 

0.947 

5 

43 

0,0CX)38 

0.765 

0.26 

7.35 

1.050 

6 

19 

0.00030 

0.593 

0.22 

7.45 

1.076 

7 

51 

0,00027 

0.553 

0.19 

7.12 

1485 

8 

27 

0,00026 

0.516 

0,185 

7.16 

1.258 

9 

44 

0,00026 

0.516 

0.185 

7.16 

1.138 

10 

26 

0,00024 

0.465 

0,17 

7.02 

1.250 

11 

49 

0,00C^5 

0.457 

0. 16 

7.00 

1.558 

12 

50 

0,00023 

0.457 

0.16 

7.00 

1,188 

13 

61 

0.00019 

0.381 

0,121 

6.35 

1.327 

14 

62 

0,00018 

0.364 

0.12 

6.60 

0.955 

15 

50 

0.00016 

0.314 

0.098 

6.25 

1.110 

16 

20 

0.00012 

0,246 

O.C76 

6.17 

1.067 

1? 

51 

0.00011 

0.216 

0,065 

6.00 

0,924 

Ifl 

53 

o.oooc^ 

0.180 

0.054 

6.00 

1.178 

n 

48 

0,00009 

0.180 

0.054 

6.00 

1.053 
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T;iBLE 4.2 (CONED.) 

(b) NEija-MJERiL AmoSEEEmC GONBIwnWfl 


SB 

No. 

Ian 

No. 

MU' 

2L 

% 

0 

% 

0 

oominited 

' ^25 

m 

Computed 0^ 
Observed 0 

z 

20 

54 

12500 

440 

4.4 

0.830 

21 

45 

12500 

440 

4 .4 

0.978 

22 

57 

12500 

440 

4 *4 

0.898 

25 

24 

12500 

440 

4 14 

0.956 

24 

21 

12500 

440 

4,4 

1.C48 

25 

22 

12500 

440 

4.4 

0.956 

26 

55 

12500 

440 

4.4 

1,000 

27 

42 

12500 

440 

4.4 

1.000 

28 

56 

12500 

440 

4.4 

1,000 

29 

46 

12500 

440 

4.4 

1.555 

50 

58 

12500 

440 

4.4 

1.158 
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Tmm 4.2 (CONTO.) 

(c) STmjM iimOSHffilIC ookditioi® 


SE 

No. 

Hun 

No, 

% 

0 

1" 

, 2“ 
k X 

L 

z 

1 

Computed 

Computed 

Observed 

31 

67 

0.00009 

0.182 

O.C4O 

4.40 

1.535 

52 

17 

0 , OOOC9 

0.182 

O.C40 

4 .40 

1.190 

33 


0,00009 

0.188 

O.C41 

4.36 

1.082 . 

34 

57 

0.00010 

0.193 

O.C42 

4.36 

0.970 

35 

41 

0.00010 

0.200 

0 ,C 45 

4.30 

1.067 

36 

65 

0,00015 

0.267 

0.050 

3.75 

1.100 

37 

54 

0,00014 

0.276 

0.052 

3.77 

0.968 

38 

60 

0.00015 

0.291 

0.056 

3.86 

0.990 

39 

29 

0,00028 

0.553 

0,085 

3.06 

O.9O6 

40 

18 

0,00028 

0,553 

0.085 

3.00 

1.140 

41 

6b 

0.00030 

0,595 

0,094 

5 . 18 

1,760 

42 

66 

0.00058 

0.765 

0,120 

3.14 

1.500 

45 

59 

0.00050 

1.000 

0.140 

2.80 

1400 

44 

58 

0,00067 

1.355 

0.0175 

2.62 

1.750 

45 

36 

0.00067 

1.333 

0.175 

2.62 

1.580 

46 

32 

0,00067 

1.553 

0.175 

2.62 

1.640 

47 

14 

0.00075 

1 .455 

0.180 

2.48 

1.770 

48 

40 

0.00154 

3.08 

0.520 

2.04 

0.612 
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TiiBlE 4.3 

coMPiaisoi OF sTMum miumm of imtigm. spbe/j) is galculmed 

FBOM EXTENSION 0FMDN1N«S APHOi'CH 


U) 

UNSTABLE ATMDSEHERIC CONDITIONS 



SR 

No. 

Run 

No. 

0' 

C^J" 

"bkic 

z 

max 

W 

Computed 

Computed 

Observed a 

z 

1 

16 

0*00149 

8.140 

1.900 

5.21 

0.516 

2 

15 

0.00105 

6,150 

1.350 

4.90 

0.655 

3 

52 

0,00060 

4,800 

1,000 

4.65 

0.465 

4 

25 

0,00075 

4,360 

0.920 

4.71 

0.525 

3 

43 

0,00058 

2,290 

0.440 

4.50 

0.615 

6 

19 

0,00030 

1.780 

0.350 

4.40 

0.658 

7 

51 

0,00027 

1.600 

0.510 

4.52 

0.900 

8 

27 

0,00026 

1,550 

0.300 

4.47 

0.771 

9 

44 

0.00026 

1.550 

0.500 

4.47 

0.710 

10 

26 

0,00^4 

1.455 

0.290 

4 »4 5 

0.780 

11 

49 

0,00023 

1.371 

0.260 

4.24 

0.940 

12 

50 

0,00025 

1.371 

0.260 

4.24 

0.719 

15 

6l 

0.00019 

1.142 

0.200 

5.91 

0.815 

14 

62 

0.00018 

1.090 

0.195 

4.00 

0,580 

15 

30 

0,00016 

0#94 1 

0.170 

4.05 

0.720 

16 

20 

0,00012 

0.740 

0.155 

4.08 

0.7 C4 

17 

31 

0,00011 

0,650 

0.120 

4.12 

0.655 

18 

33 

0.00009 

0.540 

0,100 

4.14 

0.810 

19 

' 48 ;. 

0.00009 

0.540 

0,100 

4.14 

0,714 
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WJIiE 4.5 (CONTD.) 


(b) AraOSHDSBIG COroiTIONS 

SR 

No* 

Bun 

No* 

2L 

z 

0 

mast 

z ' 

0 

Computed 

Computed 

Observed a 

' z 

20 

54 

12500 

1000 

5*82 

0.72 

21 

45 

12500 

1000 

5.82 

0.85 

22 

57 

12500 

1000 

5.82 

0.78 

23 

24 

12500 

1000 

3.82 

0.83 

24 

21 

12500 

1000 

5.82 • 

0.91 

25 

22 

12500 

1000 

3.82 

0.83 

26 

53 

12500 

1000 

3.82 

0.87 

27 

42 

12500 

1000 

3.82 

0,87 

28 

56 

12500 

1000 

5.82 

0.87, 

29 

46 

12500 

1000 

3.82 

1.158 

50 

38 

12500 

1000 

3.82 

1.00? 


T/J?LE 4.5 (CONTD.) 

(C) STABLE ATMOSHfflEIC COlOirmNS 


SE 

No* 

Eun 

No. 

0 

bkx 

L 

’ 2 

max 

L 

Computed 
c i m 

Computed 

Observed a 

a 

31 

67 

0.00009 

0.545 

0.090 

3.68 

1.117 

52 

17 

0.00009 

0.545 

0.090 

3.68 

0.995 

55 

35 -S: 

0.00009 

0,565 

0.091 

3.60 

0.900 

54 

57 

0,00010 

0,580 

0.092 

5.55 

0,790 

55 

41 

O.OOOlO 

0.600 

0.095 

3.53 

0.885 

56 

65 

0.00013 

0.800 

0.123 

3.45 

1.010 

57 

54 

0,00014 

0.827 

0.131 

5.54 

0.905 

30 

60 

0.00015 

0.875 

0.140 

3.58 

0.919 

59 

29 

0.00028 

1.653 

0.250 

3.10 

0.912 

40 

10 

0,00028 

1.653 

0.230 

3.10 

1.15 

41 

68 

0,00030 

1 .780 

0,250 

3.14 

1.655 

42 

66 

0.00038 

2,280 

0.510 

5.05 

1.442 

45 

59 

0.00050 

5.000 

0.380 

2.83 

1.415 

44 

50 

0.00067 

4.000 

0.520 

2,90 

1.932 

45 

36 

0.00067 

4.000 

0.520 

2.90 

1.526 

46 

32 

0.00067 

4.000 

0.520 

2.90 

1.812 

47 

14 

0.000? 3 

4.360 

0,560 

2.86 

2.C4O 

4 B 

40 

0.00154 

9.240 

1.02 

2,46 

0,725 
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Bosidofl ooroparinfc ths ooinpj.tcd and observed values of o in 

z 

"th© individuftl onsosj ths ©"bKjnved values of can "be compaxed 
graphic iilly with 0^, values computed by the foDanu la tion given earlier. 
Pig. 4*1 ©hows values versus i/L values. The vertical spread 
depends on tho stability conditions of the atmosphere, measured in 
terras of stability length L. The numerical value of stability 
length L approaches very large value as the atmospheric conditions 
approach neutral stability, end becomes infinite at the neutral 
oondltions of atraosphorio stability. Hence instead of plotting L 
graphically, the inverse values, i.e, 1/L values, are used in 
Pig. 4*1* Hi® observed values of a are dbown by simple dots in 
that figure. For the sake of oomparlson, tho values are 
calculated from the proposed formulation for different values of 
stability length L, 

Table 4.4 shows the values of standard deviation of vertical 
spread nt 100 ra from the source as calculated from the proposed 
forraulation. Calculations were done for three different cases 
of Table 4.2, For six different values of stability length L, in 
the oiwo of unstable ahnosphexlo conditions, a^, values wea^e 
eraoulnted from Pig, 3,6. Calculations were carried out exactly 
in ti’io oiiinu way as wore done for Table 4»2. Similar calculations 
wore ojirilod out for neutral atmosphere (using Pig, 3.5) and six 
oauuo of otabl© atmo^iere (using Pig. 3*7)* The computed 0^^. 
values at® shown by enoirolad points in Pig, 4*1 • 
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Similarly, values of standard deviation of vertical spread 
at 100 m from tiio source were also computed using extension of 
Monin’n appro, 'wh diom in the previous chapter. Table 4.5 shows 
tho viiUion for rdx oases of unstable atmosphere (calculated 
uoinr P’iff. 3 * 2 ), a case of neutral stability (calculated using 
l'’ii":. 5.4) and six cases of stable atmospheric conditions (derived 
from Pig. 3 » 3 )* The computed values are inserted in Pig. 4»1 
art triangular ixjints. 


TJBIE 4.4 

STANDiED DEVIATIONS OP VEHTICAL SPBEAD AT 100 m, AS CilCtlLAIED 
BY HiORlSED POHMULADION 

(a) UNSTBLS ATMOSHIEEIO CONDITIONS 


L , 

m 

% 

0 

1 1 

■ z 

HI 4 Ii ^ 

-1 

m X 1 0 

-10,0 

0.00060 

1 *6 

0,700 

8.75 -10.00 

-13.3 

0,00050 

1.2 

0,480 

8.00 -7.50 

-20,0 

O.OOC40 

0,8 

0.290 

7.25 -5.00 

-26.7 

0,00050 

0,6 

0,210 

7.00 -3.75 

-50,0 

0.00016 

0.52 

0.095 

5.95 -2.00 

-lon.o 

O.OOOOd 

0.16 

0,040 

5.00 - 1,00 

(B) NWJTl?AD AmOSFHlRIC CONDITIONS 

0 , ^ 

m , ■ 

2 

0 

m 

g, „ ' 

g" ' ' 

0 , ' 

m 

1 

L ' 

m ' xlO^ : 


4 *4 


0.008 


12500 


440 


G.OG 
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T.-IBIE 4.4 (CON TO.) 

(C) r, TAB IK ATMOS riffiRIC CONDITIONS 


1 , 

ffi 

% 

0 

2» 
k X 

"IT 

z 

L 

Oz f 

m 

1 

L » 

m X 10 

10,0 

0.00080 

1.6 

0.195 

2.44 

10.00 

15.5 

0,00060 

1,2 

0.160 

2.67 

7. 50 

20*0 

0.00(40 

0,8 

0.115 

2.88 

5.00 

26*7 

0,00050 

0,6 

0.094 

5.15 

3.75 

50*0 

0.00016 

0,52 

0.058 

5.59 

2.00 

100,0 

o.oooon 

0,16 

0.034 

4.25 

1.00 


TABUS 4 .5 

STANDARD DEVIATIONS OP VERTICil. SPREADS AT 100 m, AS 
C/iCUMTED PK3M EXTENDED MONIN’S AEEROACH 

(a) UNSTABIE ATMOSffiERIC CONDITIONS 


L , 

m 

% 

0 

C-L) 

<m 

bkx 

C-L) 


CTz » 

m 

1 

L ’ 

-1 

m X 1 0 

-1O.O 

0.00080 

4 f 80 

1.00 

4.65 

1 

0 

» 

0 

0 

-15.3 

0 . 0 (XS 60 

3.60 

0.70 

4.55 

-7,50 

-20,0 

0.00(40 

240 

045 

4.19 

- 5.00 

-26,7 

0.00050 

1,80 

0,53 

4.10 

- 3 *7 5 

-50,0 

0,00016 

0*96 

0,17 

5.95 

- 2.00: 

-100,0 

0,00008 

0,48 

0,08 

5,82 

-1.00 
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T/iBLE 4.5 (CONTD.) 

(B) NTSUTR/il. .VlMGSraKRIC CONDITIONS 


K t 

0 

m 



0 

z 

max 

z 

0 

Og f 

m 


^0 

X 

1 

0.008 

12500 

1000 

3.72 


0.00 ^ 

Tcl 

STABUB SMOSHERIC 

CONDITIO IB 




L , 
m 

58 

0 

L 

Blcx 

z 

max 

L 

m 

1 

x10“^ 

10*0 

O.OOOOO 

4.80 

0,62 

2.88 

10.00 

15.3 

0.00060 

3.60 

0.47 

2.92 

7.50 

20,0 

0,00040 

2,40 

0.32 

2.98 

5.00 • 

26.7 

0,00050 

1.80 

0.25 

5.10 

3.75 

50.0 

0.0(X)16 

0*96 

0.14 

3.26 

2.00 

100,0 

0,00008 

0#48 

0,072 

3.55 

1,00 
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^ * "^ *^ * of Obaerved a nd Comwted Values of Vertical Sp3::ead 

Fron Tabluo 4,2 and 4.3, and Figure 4.I, the foUowing 

ohaewations ofin bo malo 1 

(i) 'Pho iifjrooment in the values of vertical spread oompated 

by tho Extereled Monin's Approach with those observed in the unstable 
atmospheric conditions ia rather poor. The Extended Monin's 
Approach prodiots mch lower values of vertical spread, ratio of 
ooffliiutod 0^ to obiiorvod Cj, falling in the roige of O465 to O.940, 
and the disfigreoioout increases with the increase in the instability 
of tho utmouiiJieru, Tima, lower estimation of vertical spread by 
thl u method mriy leal to higher estimation of concentrations. 

(ii) In the neutral atao sphere, the vertical spread calculated 
from Monin*® AiiiTOaoh is loss then the observed one, however, 
thtJ ontimato io leas by about 10-15?^ for most of the near-meutral 

oases* 

(iii) ^ri.ng stable a'baoapherio conditions, the spreads estimated 
from the Extended Monin's Approach agree well with the observed 
ont)B nnd altw with tho spreads oalenlated from the formulation 
profiotjod by tho rfutiior for z/L £, 0.1* Again, bpyond ** ^#1, 
tho j^roemunt ia poor find could be regarded as stemming out mainly 
from the Inadequacy of formulation of wind profile for very stable 
atmofsihoH) ao ibirthor diooussod in point (vi) below, 

(iv) Kw formulation proposed by the author for -the unstable 
ataao sphorlo oonditions offers the estimation of vertical spread 



very oloso to the observed values. Prom Table 4.2, it cam be noted 
that tho veariation in the ratio of computed vertical spread to 
obtsorvod vertical spread is in the range of 0.875 to 1.558. Although 
thorii is aomo variation oi the observed vahies, keeping in mini 
tho atatiatioal nature of the obseirved values of vertical spread, 
tho agr(5omont shouH bo considered rather close. Prom Pig, 4.1, 
it can bo very clearly seen that the values of vertical spread, 
comiutod from the proposed formulation for the different degrees 
of itiotr^i 11 ty oi the atmosphere, matoh very closed to the observed 
vortical spreads of pollutant. This agreemeit is observed for 
tho entire tfvige of instability oonsidered, i.e., from near- 
noiitral ntmooiihorio conditions to highly unstable atmosphere luith 
stability length as low as (-)5.9. 

(v) At neutral atmosphorio oonditions, the vertical spread 
oaloulatod from the proposed formulation agrees excellently well 
with the obBorved vertical spread# 

(vi) It is observed that, during stable atmospheric conditions, 
Uiu vortical sproalo calculated from the proposed foirmulation 
?igrt'iO!! vory closely with the observed spreads till the ratio of 
'I/l. < U,1, Beyond this, the agxsement is rather poor. It should 
bo natud horc! that the formulation for estimation of vertical ' 
nproat i;i bnJJod on a wind velocity profile given by Eq.. (2.14 ) 

( re-^'ritUtn m) (3#4 b)), which does not hold good beyond z/h w' 0,3, 
and no other appropriate formulation is available (lumley and 
I'miof ulry, 1964 ) . aean value of vertical spresl z (in the ratio 
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of k/l) 'tho moan of particles found in the vertical section 
at a partioulnr downvdnd distance x, For the concentration 
di iibrlbu lion tiJtod in the derivation of the vertical spread, i.e., 
{?A* 

C « 0^ oxp ( - ~) 

z 

“ ®xp ( - , 

z/X. 

v./h « 0 , 3 f k/L » 0.1, o/Cq - exp C-3) - 0.05, 
xihiah raufunn tiiat conciinbrntion at z/L « O.3 is O.O5 tine axial 
aoncentration, Htrnou some particles are beyond 1h.e hei^t 
z « 0 , 3 L, which do not obey Eq,. (2.14) and derivations therefrom. 
;.n z/li inortiiu.!un laoro /ind more particles fall under this undefined 
ro^on. Till n io probably the reason for the discrepancy in the 
obsowad iffid ooraputed values. Obviously such discrepancy is only 
Xmomtm of a proper velocity profile not having been available 
and not duo to any flaw in tho approach proposed. 
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4*2 Dlopernlon Study at St* Loula 

rtilfi dloporsiott study wsis OAXiisd out at St, louisj 
StiB 30 UJ?i| from Msy 1965 to Maxoh I 965 with a priiuaiy ohjscti've 

to describe dispersion over an urban area and to relate the disper- 
sion to moteorolo^cal parameters. 

St. Louis was 00 looted ae the experimental area because 
it io locat«*d In s roaaonably flat terrain, removed from signifioant 
topofTraflvlc fontureB that would influenoe the large-scale airflow. 

MoElroy rmd Pooler (l96fia, 196Sb) reported the data and 
primary wmlysos purrortiied on the data collected during the 
diaperoion abidy over m'*.tropolitan St* Louis, ' 

4*2*1, ExiiCtrtmontal Teohnibues 

The fluorescent particle system utilizing zinc-cadmium 
ailfide was used as the tracer system, ©le tracer was dispersed 
frtm an nOTOOol diBseminator, and the dissemination of -ttie tracer 
from either of the two pre-selected sites near ground level was 
uaially 1 hour long. 

For each ©xporinent, meaairemmts of total dosage at the 
mrfflcc) wore obtfvinod on three nearly oircular.arcs at distances 
between ^ wd 10 mllei from the dissemination site. On each arc, 
sampling oiteo wore at nomtoal 6 -degree intervals, 

A roe teoro logical network was oompsed of three stations 
on tlw porlphory of ttie urban area and an instrumented television 
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towor (niimod KMOX) in the downtown area. They provided continuous 
rocojxin nf wind, tanporature and relative humidity. . The TT tower 
waa inatrumnntod at thro© levels to provide informatipn on the 
vurtioali grcriiants of wind and temperature, lerovane wind 
transmitters recorded wind speeds and directions, while aspirated 
’thormohae* operating on th® resistance principle with oopper- 
wixtj flensing eiemmts transmitted temperature data eleotrioally 
to IX multi-point recorder. Valid turhulenoe data weire not h 
availdsl© during ax^' of the tracer experiments. ■ 


Over the period of 2 years, 26 daytime and 16 evening 
oxporimunts waro conducted in varying moteoirolQgical oonditions, 

4 •2.2, Vertloal Spreads 


From the omcentration data collected for the dispersion 
experiments, MoBlroy and Pooler (l968h) calculated the standard 
deviations o„ of vortical spread for each experimental run. They 
aaaimed Oaussiin distribution of concentration in -ftie vertical 
to dorivo values. 

For each experimental run, they analysed the meteorological 
dntft find OBtimatod the standard deviation of horizontal wind 


direction fluctuations, a. , and the 'bulk^ Eiohaidson number, 

0 . ■ ' 

Th© Rio wao based on temperature and wind measurements at 
127- and 4 59 -foot levels of the TV tower, and was defined as 




B 


£ 

f 


[at/az + r^3 

■ ■ . 


(4.3) 
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whoro K 3 3 ..ooolemtion duo to gravityj aT , temperature differenoe 
totwoon top ami bottom of the layer} Aa, hei^t difference between 
top wd bottom of the layer; I, mem absolute temperature throu# 
the Inyor; dry adiabatic temperature lapse rate; z, height 
of upi^or anemomoter ; and u , the mean speed for anemometer at 
hoiKht z. 

Those Ojj Values at different distaMoes from the source 
and Rijj valuee nro reported in Table 4 . 7 . Out of the total 42 
experiincntal runs during different atmospheno oonditions, results 
of ton oxitorl wonts wore considered doubtful (loElory and Pooler, 
196f?b) nnd honoo 0 ^^ values for 52 experiments were reported, Prom 
theao 52 seta of the roailts, 4 sets were oonsidared unreli^le, 
MoElroy nnd Pooler grouped these useful 28 runs into four different 
stability oategoiios, as shown in Table 4 , 6 , 

TABIE 4.6 


STAUILITT CAffiCSORIES OP MoEIEOY iKD POOIEE 


00 f dogrooo 

h-b 

Experimoat los. 

24 ^ 

< **• Of 01 

5,7»11,19»20,25,28 

IB - 22 

< OaOl 

2,5,4,6,8,9.14,22,25,51»40 

15 " 20 

Of Ol 

16 , 24 , 55 , 42,45 

B *» 13 


'' 2 , 18 , 53 , 57 , 41 " 


Thoy plotted values of versus x for these four categories and 
drev^ bo 0 t«^it lines for each category by statistical analysis, 

Tfible 4,7 reports Og values for each of the useful experiments in four 
ontogorlOB of tiio atmo spheric stability. This table is compiled 
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TJfflIE 4.7 

iniMflARY UK DUIIKHDIUN DATA AND METEORDLOGICAl INDICES FOR 
m\ LOUIS DISIERSION STUDY 

(A) i * Oq I + ? Rijj < - O.Oi 


Exi^oriniGnt 

No. 


Oq 

at 127 ft 
degrees 

X , 

’ m 

m 

Diurnal 

period 

5 

••0*25 

53 

1937 

1943 





4267 

1672^ 

Day 




8065 

1485^ 


7 

u,03^ 

29 

1994 

997 





4258 

4237 

Day 




7999 

2100®* 


11 

m ■ 

25 

1998 

323 





4197 

814 

Day 




7964 

1511 


19 

m 

25 

712 

149 





5198 

999 

Day 




6440 

318O 


20 

0,(X)^ 

26 

825 

480 





3284 

11930 

Day 




7421' 

3445®" 


25 

"^•02 

31 

698 

5 16 

a 





3245 

1204 

Day , 




6384 

1106^ 


28 

-040 

24 

2019 

2148 





4258 

3466 

Day 




1596O 

1241^ 
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TAB IF, 4.7 (CDNTD.) 


(B) Oniofwry Tl 

* 0Q t 18® 

- 22® j HIb < 

-0.01 



Export fflcnt 

KLb 

a§ 

at 127 ft. 

4 m» 

X , 


Biumal 

No. 


V f 1 V * 

degrees 

m 

m 

period 

2 

- 

21 

732 

152 





3152 

559 

Day 




6445 

4154 


3 

-0,05 

17 

761 

56 





3156 

220 

Day 




6770 

890 


4 

—0, 05 

19 

882 

77 





3525 

529 

Day 




7590 

1429 


6 

-0,13 

20 

2022 

2901 





4286 

'St 

1752 

Day 




7958 

1154 


0 

-O.CTf 

20 

1994 

263 





4171 

765 

Day 




• 7988 

952 


9 

■"0,06 

21 

715 

101 





3296 

6 25 

Day 




6607 

1238 


14 

-0,(4 

19 

2050 

123 





4157 

459 

Day 




7932 

657 


22 

-0.51 

17 

945 

1(4 





3378 

1771 

Day 
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ivjni'i -1 .7 (odNTO.) 


Export tnont 

No. 

Hij 

at 127 ft 
degrees 

, X 1 

’’ m 

az» 

. m 

Diurnal 

peilod 

25 

-0,10 

10 

825 

99 





5274 

984 

Day- 




7601 

4112 


31 

-0,11 

19 

940 

56 





5571 

2535 

Day 




1357 5 

1982 


40 

-0,05 

21 

755 

147 





3207 

872 

Day 




6445 

1211 


(C) CatofJX^ry III * Oq » 

15® - 20® 1 Big + 0,01 



16 

0,00 

1? 

710 

64 





3241 

277 

Evening 




64C4 

470 


24 

0,01 

20 

949 

44 





3584 

345 

Day 




7434 

24I8 


35 

-0,01 

15 

754 

250 





3245 

662 

Day 




6395 

597®' 


42 

0,00 

15 

1989 

97 ■ 





4205 

154 

Evening 




7501 

204 


43 

—0,01 

15 

757 

65 





5186 

3C4 

Evening 




6607 

420 
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TiJIiS 4*7 (CONID.) 

(D) Catogttry lY j Og , 8* - 15« } Bij > O.Ol 


Experiment 

No. 

“3 

Cg 

at 127 ft, 
degrees 

X , 

m 

m 

Diurnal 

period 

12 

m 

13 

2052 

98 





4294 

86 

Evening 




8058 

224 


18 

0,14 

0 

914 

24 ^ 





3577 

65 

Evening 




7200 

173 


35 

0.05 

13 

959 

65 





3398 

144 

Evening 




7C47 

109 


37 

0,01 

13 

753 

68 





5162 

171 

Evening . 




6514 

345 


41 

0,12 

12 

627 

66 





3223 

102 

Evening 




7459 

98 



* Values considered to bo aignifioemtly affected by a 
limited mixing layer in td® vertical. 


^ Values oonsidered to be doubtful in view of otbor independent 
moaaares of stability during laie ex^riment. 
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4.2.3. Evaluation of Rou^eas Length 

For tOie approach proposed in Chapter Three by the gpthor, 

the estahliehment of roughness length is very iu^rtant. Since 

no reliable estimation of z was available for the urban St. liouis 

0 

area, this was determined as described in the following paragraihs 
from the wind profile data at near-neutral conditions. 


Out of total 42 experiments, iiiose experiments were 
considered near-neutral ones which had bullc Eichardson number Mg 
in the range of + 0,01, Such experiments were Nos. 16 , 24 , 35, 42 
and 45* Wind velocity records at two heights of KMOX - TV tower 
v/ore analysed for the five selected experiments and z^ was estimated. 

For the neutral or veiy'^ear-neutral stability, the mean 


wind speed can be given as 

(4.4) 

k 

0 

d is *th6 displsDsnisnt leugthf which depsn-ds on *tihs 

average height of buildings in the area, as discussed in Section 
2. 2. 3, 2 and k is Von Karman’s constant » 0.4. For St, Louis urban 
area, d can be estimated as 5 ® l^t, (Hanna, 1969 ). Writing Eq., (4.4) 


for two heights 






(z2-d) 

^0 
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or, 


or, 


In (z^-d) - In 2 



«*. In (ZgMi) - In In (z^-d) - In 

Ug In (z^-d) - In (z^-d) 


Ug 


log z. 


Ug log (z^-d) «u^ log (Zg"^) 


(4.5) 


U« •• Um 
2 1 


Hero and Ug can be taken as half- 4 ionrly averages at two heights. 

The wind spood data at KEOX - TV tower, as reported by 
MoElroy and Pooler (1968a), axe lO-min averages. These data were 
used to obtain half -hourly averages of wind speeds at two heights. 
These average wind speeds are in miles per hour and are reported 
in Table 4.8. It was observed from comparison of wind speed data 
found from KMOX » TV tower with that from pilot balloon measarements 
that the wind speed data during experiment No. 43 were inconsistent, 
and henoe were not included in the analysis. Out of the four 
os^oriroonte used for this analysis, wind data of anemometers fitted 
at 127 ft. and 255 ft were used for experiment nos, 24, 55 and 42, 
whoroao during experiment no, I6, anemometer at 255 ft was not 
fitted (Mollroy and Pooler, 1968a), and henoe wind speed rooords 
of 127 ft and 459 ft were used. Using d « 50 ft in all oases, • 
value of roughness lengidi was oaloulated for each half <^ourly 
interval. Those values of z^ are reported in Table 4.8. 
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TmiE 4.8 

of as^ FBOM HAIF-4I0UHLY wind data,, using EQ,(4,5) 


Experil® nt 

No, 

Time 

at 

127 ft, 
mph 

Ug at 

255 ft, 

mph 

z , 

0 ’ 

ft 

16 

20,20 

m 

20.50 

14.50 

21.00®* 

346 


20.50 

m 

21,20 

14.25 

' 21.75^ 

5.75 


21.20 

- 

21.50 

13.50 

21.75®" 

8.12 


21,50 

- 

22,20 

13.50 

21.50^ 

7.50 


22 .20 

- 

22.50 

16.50 

24.00®' 

5.55 


22,50 

mm 

23.20 

17.75 

■ 24.25® 

1.70 


2 3.20 


23.50 

17.25 

24.75® 

3.02 

24 

JO. 00 

- 

10.50 

16.00 

20.75 

5.45 


10,30 

- 

11.00 

15.50 

20.00 

5.40 


11,00 

- 

11,50 

17.00 

22.25 

6.50 


11.50 

m 

12,00 

17.00 

22.75 

740 


12.00 

•m 

12.30 

17.25 

22,25 

5.25 

35 

12,00 

m 

12,30 

14.00 

16.25 

049 


12.30 

- 

15.00 

14.00 

• 16.75 

1.28 


15.00 

m 

15.30 

15. 00 

19.00 

5 .98 


13.50 

- 

14.00 

14.25 

17.00 

1.28 

42 

20,00 

m 

20,30 

11,25 

15. 00 

7.84 


20,30 

■m» 

21,00 

11,75 

16.00 

9.75 


21.00 

m 

21,30 

14,25 

19.25 

8.70 


21,50 

m 

22.00 

14.50 

19.75 

9.35 


22,00 

m 

22.50 

16,50 

19.75 

1.33 


22.30 

m 

23.00 

14.25 

16.50 

0.49 

a t ItioBo war® at hsip^tit of 459 ft, 


Ttim& wer® at hsip^t of 459 ft, 



TiiBlE 4.9 

ME/JJ MD BUmm DEVIATION OF z 


« 

B » 

0 ' 

ft 

z^-^ean z , 

0 o' 

ft 

2 ' 

(z - Mean z ) , 

0 „ 0 

ft 


346 

-1 .42 

2.016 


5.75 

0.87 

0.757 


0,12 

5.24 

10,500 


7 .50 

2.62 

6 .864 


3.55 

-1.33 

1.769 


1.70 

-3.18 

1O.112 


3.02 

-1.86 

3.46 0 


5.45 

0,57 

0.325 


"^40 
6.50 
Y *4 ® 
5»25 

0.49 

1,28 

3.9a 

1,28 

7.84 

9.75 

8.70 

9.35 

1.55 

049 


. 107.^ 


0.52 

142 

2.52 

0.57 

*4.59 

-3.60 

-0.90 

-3,60 

2.96 

4.87 

3.82 

447 

-3.55 

-4.59 


0.270 
2. 016 
6.350 
o;i57 
19.272 
12.960 
0.810 
12.960 
8.762 
25.720 
14.592 
20,000 
12.603 
19.272 

' ' ■ ' ' r 

E (z^'^ean 


.= 189.5?/ 


Meai a. 


12US. 4.88 ft * 1,488 a 
22 
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• 955^ confidence interval ■ Mean z + 2.Ce c 

X 

- 4*88 + 2,06 X 0.64 

- 4»88 + 1.55 ft 

« 6,21 ft or 5.55 ft 
i»©. 1,895 ™ on 1.082 in 

This shows that there ia 95^ probability of the value of falling 
in the rang© of 1,082 to 1,895 m. 

The tttrfaoo roughneao length, z^, is sometimes oalculated 
vdth the assumption that displacement length d » 0, !lhe Eq., (4,5) 
booomoa 

Up log z, - Ui log z 

gp ' """S" Z" " (4.8) 

2 1 

Using and Ug, values from Table 4.8, was oalculated for each 
half -hourly interval. The aXithmetio averege value of z^ worked 
out to be 10,5 ft (« 3.14 ®). ®3.is shows that the estimation 

of z^ from wind data assuming d ■ 0 leads to over-estimation of 
value, 

O' 

It is interesting to compare the mean value z^ computed from 
the wind data during near-neutral atmospheric conditions, and z^ 
value given directly by approxlaation foimlae given by iSitzbaoh 
(1961, quoted by Hanna, I969). According to 'Kitzbach, as discussed 
in Section 2, 2. 3. 2, Eqs# (2,6a and 2,6b), 

d 7 .a‘" 0»3 

^ « .(A) 


(4.9) 
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F* * U) (4.10) 

wh02?o d. iB displsoQosii'b 161181311} h. is ilis aivsiag© iisighi; 

of the iMi Mings, and i. is the ratio of the total area of a given 
region to the area of the bulMings in that region, TaMe of I 
Uoo u at ally in the range of 2 to 10 for urban areas# 

Mollroy and Pooler (1968a) have showi (their Hg, 15) the 
land UBago in vicinity of M3X - TV tower. Prom this, it can be 
assumed that tide average height of the buildings in the area is 
around, 60 ft. With d * 3O ft, from Eq,, (4»9)> natio A can be 


worked out, as 


M (A) 


A « (2) ** « 
Using A - 10 in Eq. (4.IO), 


s t 

O 0 

F "IPS 


■ 60 w » 4.75 ft “ 145 m 

It is observed that roughness Imgth calculated from 
Eq. (4«10) in of the same order of the value of found from wind 
data. The former vaM© li©s within the 95?^ confidence interval of 
deMvod from wind data, and is very close to mean 
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For thi present study of diffusion, average value of rou^ess 
length given in Table 4.9, as 1.488 im (« 4.88 ft) om be approxim^d 
to 1*50 HI and used, 

4*24* Stability Lengths 

For the St, Louis dispersion study, the stability length L 
for nny of the oxperimint has not been derived. It was necessary 
to establi* L for the present study. 

As doBoribed in Section 4.2*2, MoElroy and Pooler (l9$8b) 
used and as the stability parametors and derived relation-* 

ship between and downwind distanoe x for four different categories 

of at»s]^@rio stability. Studying these graphs (Fig, 2 of McElroy 
and Pooler, 1968b), it oan. be observed that there is a considerable 
scatter of data points about best-fit lines, Henoe, instead of 
0 stab lilting stability lenglh L and then evaluating Cj, from 
theoretioal considerations to compare with Cj. derived by McElroy 
and Pooler for each individual experiment, more meaningful 
oompiirison oouM be done with the beet-fit lines of plot 

for oaoh of the four oategories of atmospheric stability* Bms, 
stability lengths L representing each of Ihe four stability 
oategortes axe ©stadjlltiied as shown below, : ; ^ . : 

As stated in Section 2,24.2', 1^® standard deviation 0^ of 
wind dlreotion fbaotuations oan be takoi as a predictor of 
turbulaioe. In the pre^nt case, it is used to establish stability 
length L, Following Inone (1959), 2<1. (2.24®), 
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" - (4.11) 

Wt 

whoro e ia a constmt end 0^ ia in radians. The average values 
of Og for the four categories of stability were taken as shown 
in Table 4.10. 


TilBIil4.lO 

msum TiUJES OI Og IOR tour STABIIITI CAaSOOEIES 


Category 

Hg 

stability 

conditions 

Average 0., 
degree 

I 

< «0,0l 

Unstable 

27 

II 

< ••OtOl 

Unstable 

20 

III 

0* 01 

Neutral 

18 

IV 

> 0.01 ■ 

Stable 

12 


'Rno Eq., (4.11) modifies to 


tr 



(4.12) 


for in degrees. 

For tho ntjutral stability condition, the mem wild speed u in 


Eq,. (4.12) can be subatitutid from Eq. (4.4)i so 1^®* 




"e * «5 - u, 

ta (— ) 


k - 


where d ia the *©ro-pl«»o displsoeBant lengtih (•> 50 ft) and k is 
Von Kamm’s oonstmt (■ 04)* 
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• - 18^ 


2.3 log (~) 
0 


♦ .# 




or 


0 - 0,1 0. log 

V 


(4.13) 


1?ho flt«!dft£d deviation of wind direction finctuation was 
jnoawrod at ttte hoight of 12? ft. SalDstitating mmeiioal values 
of »,d and for Og- 18“ 

0 • 0.1 X 18 X log 

•« 1 «8 log 17 •55 ■ 

m 2 f25 

It oan feo recalled 1diat, for any stability condition, mean wind 
at any hei^it » is given by lq.» (2.16) as 


u - ^ IfC^^A) “ f( 
Substitating I<1. (4.14) in E<1. (4.12) 


(4.14) 




mmmmm « — 

0 100 u 


$k 

f(z/L) - f(2QA) 


(4.15) 


For tom ml... of o, (.t Wi#* .)- 8 
longth . om ^ ^oruostoOI. oase, toriog tof at 

height 127 ft .h.™ “j ooiweotion 
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TOUU be provnlllne (V(-l) > O.O 0 , see. 2 . 2 . 3 . 3 ), u ie giva, by 

Ecl* ( 5*39) aa 


u 


^ [04S - 2.5 leg (zj{^)) . 1.3 (e/(-L))-'/5 ] ( 4 . 16 ) 


Waiting 2 in th@ Iq.# (4*16) as above the zero^-plane, its effective 

vnhiii ** (z-cl) 




* .■» # 

, . u - ^ 


[ 04 fl - 2.5 log (zy(-i)) - 1.5 ((z^)/(-x))-V 3 3 

... (4.17) 

or , f(|) - f(-^) - 048 « 2.5 log (Zq/(-1) - 1.5 

• •• (4.185 

gubstituting Eq.. (4.18) in Bq. (4.15)» for unstable atmospheric 


conditions 


ek 




® 0.48 - 2.3 log (z /(rL)) - 1.3 ((z“<i)/(-€.))'''/5 


Insorting vahios of g ■ 2.25, ^ » 127 ft, d = '^0 ft, (z-d) » 97 ft 
•< 29.6 m, « 4.88 ft « 1.50 m, and k » 04 

^ 2.25 >« 04 

^ *185 '“■idS'-iiTTog 

-1.3 (29.6/(-l))'^^^) 

0,9 


0*48 ** 0#405 + 2*3 (^L) 


1.3("L)“'/^ 


(29.6) 


W 
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0.9 


0.075 + 2.3 log (-1.) « 0.42 (-1)V3 
B’or sterility oatogory I, CTq «• a?* , 

0.9 


(4.19) 


27 X 


180 


0.075 + 2.5 log (-1) - 0.42 


/. 0,075 + 2,3 log (-1) - 0.42 ^ 

or, 2,5 log (-L) - 0.42 1.855 (4. 20) 

Solving Eq . ( 4,20) by trial -and -«rxox, the stability length 1 works 
out to bo -20 ro. 

Similarly for stability category II, » 20® | substituting 
this in Eq . C4,19)f ^ simplifying 

2.3 log (-L) - 0.42 . 2.495 (4*21) 

wMoh on trial solutions gave L » - 65 a* 

For stable atmospheric conditions, u is given by Eq. (2,14) 


as 




k z ’ 1 

0 


S 0^1 


(4.22) 


whore w is a constant • 6. Hitting El, (4*22) in Eq. (4*l5)j for 
stable atroosphores, with adjustment . . 

0k 


“e “ TsS ■ 

2.5 log-^ +« i — ) 

0 


<4.23) 
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As - 12° , 0 « 2.25, z » 127 ft, d = 50, (z-d) = 97 ft = 29.6 m, 
« » 6, » 4.88 ft « 1.50 m aM k = 0.4, the stability length L 

works out from Eq. (4.23), as 128 m. 

Table 4.11 summaries the stability lengths for the four 
stability categories. 


TABIE 4.11 

STABILITI lENGTSS FOR POUR STABIIiTTY 
CATEGOHIBS 


Category 

Og , deg 3 ?ees 

L,m 

I 

27 

-20 

II 

20 

-65 

in 

18 

00 

IV 

12 

.128 


4*2.5. Vortloal Spreads from ProTX)sed Formulation 

Having established the values of ard h, the vertical 
aproads for any distence downwind can be derived from the 
thoory dosoribod in Chapter Three, These values of vertical 
spread can be compared with the observed spreads, figures 4*2 
through 4.5 show the comparison betweai these values. In each 
oatofjc^ry of stability, the 0 ^ values obtained by McElroy and 
Pooler (1968b) arc daown corresponding to their downwind distances, 
X. MoElroy-Poolor also showed best-fit lines for each category. 
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which aro shown hy solid lines. 

For known values of and Jj, the standard deviations of 
vortical spxesd, were calculated from the proposed formulation 
as well as from the extension of Monin* s approach. The procedure 
of finding out ihese value was the same as descrihed in Sec, 4*1 • 5 * 
These vaLuos are shown in the gra]^s . 
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of Obaerved and Computed Values of Vert ical Spread 

Itoni i'igures 4«2 through 4*5 f the comparisoEi. of observed 
vertical spreads aiad computed vertical spreads can be made. The 
following observations should be noted: 

(l) The Extended Mbnin*s Appxsoaoh gives smaller values of 
virtioal spreads for all atmospheric conditions. As shown in 
Fig. 4.2, for very unstable atmosphere, the values of vertical 
spread oomputed by lonin's approach are much smaller than given by 
lo Elroy and looler. As the downwind distance x from the source 
increases, the divergence increases. 

(il) Ibr tite moderately unstable atmospheric conditions also, 
the values of vertical spread predicted from the Extended Monin's 
Approach continue to fall below McBlroy-loolers ‘ values. The 
discrepancy, however, is less than that observed in very unstable 
atanosphare. 'Hais observation is in confirmation of PPG values 
where also, Monin's approach showed that higher the instability, 
wider was disagreement. 

(iii) Fig, 4.4 shows vertical spreads for neutral atmosphere, 
liionin's approach, he3?e again, gives smaller values of vertical 
spreads, which was also in line with the observation made for PPG 
valiws. The discrepancy increases with the downwind distance x. 

(iv) Monin's approach applied to stable atmospheric conditions 
at St. loulB, gives results which are closer to McElroy-Poolers ' 
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valueo, ao shown in Pig. 4.5. However / as discussed earlier, 
applio ability o£ mathematical formulation for very stable 
atmoopheres is doubtful. Although Monin’s values are better in 
comparison in this case, his values for very stable cases of PPG 
Bite were poor in agreement, 

(y) fbr very unstable atmospheric conditions shown in Pig. 4.2, 
the vertical spreads computed from the proposed formulation agrees 
closely with the ones given by McElroy and Pooler (l968b). At a 
distance of around 600m from the source, the discrepancy is maximum, 
via#, observed spread is 12Qm whereas computed value is nearly l6Qm. 
As the distanoe downwind increases the agreement be-fcween computed 
and observed spreads improves, and beyond l50Qm, the agreement is 
exoellenti 

(vl) It is observed from Pig. 4.3 that for moderately unstable 
atmosphere the formulation proposed by the author gives vertical 
spreads which are very closely in agreement with Me Elroy -Poolers' 

VttlUGS, 

(vil) 'JSi® proposed formulation gives values of vertical spread 
for neutral conditions of the atmosphere which are in excellent 
agreement with MoElroy-Poolers ' values, as shown in Pig. 4.4# 

(vili) Bbr the stable atmospheric conditions, the estimation of 
vertical spreads by the proposed formulation is rather poor. Oliis 
is, again, due to the fact that the wind profiles beyond V]j=0. 3 has 
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not baen properly expresaed mathaaatically. looking to the nature 
of slope of the times drawn in !Pig. 4.5j it is observed that for 
X < 60Chij 0 ^ values given by the foimulation would be comparable, 
but as distance x increases, more and more particles of pollutants 
reach the undefined region beyond z/lf=0.3 and the poorer is the 
agreement. 

(ix) Itor the case of stable atmosphere, the implications of urban 
heating art not property investigated as yet, and hence it is not 
very surprising that the vertical spread predicted by the proposed 
formulation is not in good agreement of McElroy-Boolers’ values. 



Olmpter Five 


IROCEKJRB TO IREDIOTION OP CONCMIRAillOIS OP lOLIUTAUTS 

IN AIR UNDER DIPFERENT METBOROIDGIGAL CONDITIONS 

Oi6 HiEliheDia1;ica,l equations (and tables and charts) developed 
in Chapter Three can be used to predict vertical spreads and hence 
air pollutant concentrations at various points. As seen in Chapter 
Bbury euoh predicted spreads and concentrations closely agree with 
the field observations. It would hence be useful to enumerate the 
steps involved in the prediction of concentration of air pollutants 
at any point downwind from a source, following the proposed mathema- 
tioal models. 

5.1 Miorometeorlogloal Data 

According to the proposed mathematical model, determination 
of surface roughness length and lonin-Obvikhov stability length 
D i© very essential. The procedures for their determination 
depend upon the type of available micrometeoro logical data. 

IBiree sets of miorometeorological instrumentation are 
considered* 

(i) Case one* instrumentation consists of wind speed records 

at two helghtej temperature records at two heights and continuous 
wind direction records at one height. 
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(ii) Case .fwos Temperature recording instruments are available 

lit two along with continuous wind speed and wind direction 

ricorcis at one hei^t* 

(iii) Case Kn^ee* Ihe miniiaum of micrometeoro logical instrumentation 
is required in this oaa6| viz., temperature records at two heights, 

and and vdnd velocity records at one intermidiate height Zg. 
height Zg is so selected that it is a geometric mean of 
and itS# Zg « In a practical case 1m, 4m , 

and hence 2m* 

d 

Ualttg data from any one of these three cases, surface roughness 
length and stability length L can be evaluated as described 

below. 



Check the typ© of the surface rou^ness found at the site 
of difi'usion study. If the site is open unobstrix ted level 
terrain, the zero-plane displacement length d = 0. If the area 
oontaioi moderate to large obstructions, the surface roughness would 
be large, and the zero -plane displacement should be taken into 
considerations and can be approximated by Eq. (4.9), i.e,, 


d 


ix 



(5.1) 


where h is the average height of ttie obstijuctions (including 
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build lng«, trees, etc.) found in the area, A is a constant 
which depcjndc on the amount of area occupied by the obstructions 
in tteit region and for urban areas lies between 2-lO.When the 
area occupied by buildings in the area is large, value of A 
would is smaller. localities like very thickly populated residential 
or coimiicrclai complexes have A nearing 2A, while relatively less- 
thickly populated areas have higher values of A. 

It is im.perative that the heights at which wind and 
tumpextiture records are obtained should be more than d. This 
ia very important in mlcrometeoro logical observations within the 
urban area. (VUnd velocity profiles and all the formulation 
derived in the previous chapters apply only above the zero-plane. 
Below toe zero-plane, i.e., at height z < d, the wind velocity- 
profiles treated earlier and foimulation for spread do not apply. 
Bollutont distribution in the local geometry of streets, lanes and 
by -lanes cannot be eathaated by the procedure described here . ) 

5. '5 Determination of Surface Roughness length . 

Surface roughness length z^ for the site can be 
detormined as shown below, 

5«3»1 Instrumentation as Case one . 

(i) Select periods duriog which neutral to veiy-near- 
neutral conditions of atmosphere prevailed. (Hais is done as 


follows 
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At iiBUtral or adiabatic conditions of atmosphere, the 
ttffliporaturo Ifipee rate ia 1.0°c/100itt = 5.5°P/lOOO ft. (Magill 
tit ul, 1956)* If the heights at which two temperatures are 

available are and Zg ( Zg > z^ ), for this interval of 
height I adiabatic lapse rate of temperature correponds to 

(zo- z. ) (m) 

I.Ox — = i (5.3a) 

100 (m) 


in dwgrtus Qentigrades, or 


(z„- zj (ft) 

5.5x— ^ J — 

1000 (ft) 


(5.3b) 


in degreei IhJarenheit. 

When the fall in actually observed temperature between 
these two heights is veiy-n early this value, neutral to very- 
near-nwtral conditions of atmospheric stability are assumed. 

(ii) Obtain the vdnd velocity records of such periods of 
neutral to veiy—near^neutral conditions. Calculate the half- 
hourly averages of wind velocity for these periods, and denote 
thaa as u^ at height mO. Zg. Ihen, the 

Burfaot roughnoas length is given by Eq. (4.5 )> rewritten here 

as : 
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u log (2 - u log (z.-d) 

log 2 ^ —U- ] _2 (3^2) 

Ug - 

5»3»2 IttS'tii*uni6Q,tia1;loii as Pass Two and Ga.se Three 

In those cases data to compute wind profiles are not 
available* 3Qie surface jxiughness length has to be approxi- 
mately QstimJited from Eq. ( 4 . 10 ) 



where h and A are appropriately aesxmied as discxassed in 
Section 5.2* 

5.4 Detarmination of Stability Length 

At the time of any diffusion study, stability conditions 
of the ateiosphere are to be quantitatively specified. This can 
be done in terms of Monin-Obukhov stability length I. Ihe 
method of determination of 1, again, depends on the type of the 
mieitjiaeteorological instrumentation available. 

5.4.1 Instrumentation as Case One 

(i) Ibr the periods of neutral to very -near-neutral 
conditions of atmospheric stability (as decided in Sec. 5.3.1 ), 
taka the wind direction fluctuation records. Calculate the 
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aiittanetic laeatt dirsctlon of th© wind, th© deviations of individual 
■ruf-idinKfi troia tho mean direction. Square the deviations and sum up. 
Thu otundard. deviation, o„ , is given by 

r ^ n 2 -ti/2 

I (0 “Mean e ) 

Cg = 1 (5.5) 

whuro 0 is tho Individual value of wind direction, and N 
number of obaorvatione for the period. 

(ii) Svniuttta the constant g given by Eq. (4,13) under neutral 
stability conditions, as 

0 » 0.1 X 0 X log ( ) ( 5 . 6 ) 

0Q being in degrees , and z being the height at which wind -direction 
is recorded. 

(ill) Similar to step (i), determine (in degrees) for the 
period at which stability length 1 is required. 

(iv) If it is a day time experiment, unstable conditions may 
bo tjQcpoctod, and the following equation can be used to determine 
L by trial solution 

ir B lE 

pf V ms ■' " ' ' ' ■ . /z (5.T ) 

® 0.48-2^^1og(zQ/(-I>))-1*3((z“d)/(-L)) 

k bflog a universal constant of value 0 # 4 * . 
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tor the evening or night time experiment, inversion 
QoMitiom may be prevailing, the equation can be used as 



B xk 

2.3 log ( f^) + cc( 


where « a 6 und k = 0.4. 


(5.8) 


Substitute appropriate values of B, a » z and d, 

6 

and obtu.in value of stability length I. 

S»4*2 Inutrumontation as Case !l?wo 

Hero also wind-direction fluctuation data are available 
and stepo to be followed to evaluate stability length L are the 
same as onimerated in Sec. 5.4.1. 

5.4.3 Inatr mentation as Case Ohree . 

Here wind direction records being not available, the 
following prooedur© is adopted. 

(i) Obtain temperature records, t(z^ ) and tCz^), at heights 
S 5 ^ « lia und 4ffl and wind records, uCzg)? at height Z 2 = 2m. 

.Mod tile average values of each of these, Find absolute mean 
tuiipe nature 5* between height z^ and z^. 

(li) Evaluate stability parameter B as given by 

B 

I U (Zg) 

where g is the acceleration due to gravity. 


t(z^) 


(5.9) 
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(ill) ^ trial solutionSj obtain, stability length L from the 
equation (iteiuohi, I 96 l) 


B as 


1 

I 


3 -j + 1.368 


(2 ^ + 2.3 log j-f 

% 


( 5 . 10 ) 


whfer® <* la a constant) which can be taken as 0.6 for unstable 
alaaoophur© and aa 6,0 for stable atmosphere (See. 2 . 2 . 3 . 3 ). 


ffiia above procedure is only for open country with d = 0. 

If d 0, z,^) Zg, can be taken as (z^-z^-d), (z^-z^-d) and 

(zj-z^-d) rospactiVBly. Hence temperature records should be 
at z,j» (Z|^+ d + 1 ) ffl) and Zj(z^+ d + 4) m, whereas wind anemometer 
should be plaosd at Z2*(z^ + d + 2) m. 

5*5 Estimation of Vertical Spread 

Once the surface roughness length z^ and stability length 
L are known, the vertical spread of pollutant" at any downwind 
distance x from the source can be estimated from the formulation 
given in Chapter Oltree. 

2-p 

(l) Calculate the ratios 2 q/|i 1 and k x/jlj where k = 0 . 4 . 

2^' ■ ' 

(ii) Ibr the pair of values of Zq/I i*! x/\ls\ , selecting 

appropriate graphs of Mgs. (3.6) or ( 3 . 7 ), value of z/| L | can 
be read off. ^ 

(ill) Obtain z f 100 the ratio i/ 11 ' 1 * 
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As the surface roughness length fbr a particular 
Bite remains almost constant, the Mgs. (3.6) and (3.7) can ■■ 
ro-drawn by denoting z/z^ as ordinate and x/z as abscissa. 
Pigs. (5*1 ) and (5.2) show the re-evaluated graphs for different 
conditions of z /III • For a known value of z , valiies of x 
corresponding to values of the ratio x/z^ and similarly for z 
can be pemanently recorded on Mgs. (5.1 ) and (5.2). Ihis can 
facilitate the use of the graphs, Men for any ratio of Zq/|Ii|j 
at a particular x, value of z can be directly read. 


5*6 laitlaatlon of lateral Spread 


(i) Prom the known values of u(z^), I and z^, obtain 
the value of u^ from 


% 


u w [ 0.48-2.3 log (zQ/(-Ii))-1«3((z^-d)/(*'l‘)) (5.1l) 


u, 


'# 


or u » {2.3 log 


z^-d 


^ a 


(z^-d) 


(5.12) 


wh«r© Eq. (5.11) is for day time tnostable and Eq. (5.12) is for 
night time stable atmospheric conditions. 

(ii) At acy height z where concaatration of pollutant is 
required to be computed, evaluate u (z) from %. ( 5.11 ) or 
(5.12) depending upon stability conditions. 

(iii) If the micrometeoro logical instrumentation is of case One 
or Case Iwo •mlue at height z^ would be known (Sec. 5.4.1 ). 
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!rht‘ valuo of oTg at height z is 


'^e(z) u(z.) 

Jg) « - ' 


^(z) 


(5.13) 


Ibr instrumwtatioa of case Three, the 0 ^ value can be given by 


$u^ 

Oft(z) “ 1 — “ 

® U(2) 

whoro gi om be aesumed as 2.25. 


(5.14) 


(iv) Rind the averaging time s for reaching upto distance 

X by 


X 

u. @ 


(5.15) 


Here u ia what is calculated in Step (ii) above, and 3 = 4 
(Hlno, 1968b), 

(v) Rind the sampling time t for vshich the sample of 
pollutant would b« collected. 

(vi) Ibr the values of s and t » from lig. 5.3 (after Hino, 
1968b), obtain the value of 

i X [0.1 Ja. 

(vil) Rrom Eq, (3»59) 

Oy - 


(5.16) 
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5*? Mtte$j.on of OonoentratiQBi Distr±!mtlon ■ 

mm Sq. (3.70), for any knovm value of quantity of 
pollutanl; released per unit time, q, the concentration at any 
point (x,y,z) oan be oaloulated 

- Q .2 

0 (x,y, 2).— 3 ^ exp [-•! ( 5 . 17 ) 

U , CJ . z 

y y 

«q, (5.17) is for a source located at ground-level. 

® 3 r an elevated source, the concentration of pollutant 
at any point (x,y,z) can be computed from (Turner, 1970 ) 

e (jc,y,a,H) • 1 - 2 “ 

{expt-~ (— )^I+ exp [-{ ( 7 ^)^ 3 > (5.18) 

i^trt 1 is the effective height of the source, 

Tbtjs fbr any given meteorological conditions and pollutant 

source data, the pollutant concentrations at various points in the 

1 ' 

air can be predicted. 2 he effect of various source control 
measures or of changed m4tebrologloal conditions can also be 
predicted to enable planning of a successful pollution control 


programme * 



Oljapter Six 

mm STcmms done u kanpde 

6*1 Introduotion 

Bss foittuls.tlon developed by the author is for a souioe at 
or vexy near grouiad level* fflae field data used to validate the 
formulitioa were also for eouroee at or near ground level* If 
there ie a point eouroe at moderate height above the ground level, 
this method jsmy still give praotioally comparable results, 
fflaeoretioally, however, some modifications are warranted, but without 
going into much of theoretical modifications, applicability of 
this formulation on the spreads at relatively short distance away 
from moderately elevated sources was checked. Ibr this purpose, 
eooe diffusion studies were carried out in the city of Kanpirr. 

Kanpur city is situated on the right bank of the river 
Gtoga# Its geographical location is aromd the latitude of 
26® 30'N, and the logltude of 80° 20'E. It is situated in the 
gangatio plains, having flat terrain of very moderate slope. 

!aai population of Kanpur is nearly 13 lacs, which makes 
it the eighth. most populous city of idae oountry. It is one of Ihe 
important induslwlal centres, <if India* , dlie . campus, of Indlatt-Jkxstitute 
of !Bi®hnalpgy at Kanpur is located at a distance of about 5 km from 
the outskirte of the oily. 
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^•2 MftMlPg of milutlon Studies 

Ibr tils prssQiit study , two sites with different surface 
roui^osees were selected. Site One was near a power-house 
situated near liie bank of river (hence, named Riverside lower— 
House), !Eh@ prevailing winds are fiom the river side, and 
smoke emitted by the power-house spreads on to the built-up area 
of the oily aroiaid Ijeathar Research Institute. Second site, named 
Sits 3?wo» was within the campps of Indian Institute of Technology 
near the boiler house in the Southern laboratories. 

It was decided to use sulphur-dioxide as the tracer, as 
the coal bufflat by the power-house has 0.495 sulphur content, and 
smoke contains easily measurable concentrations of siilphur-dioxide. 

Bie proposed method of estimating vertical spread is 
valid only in the surface layer of the atmosphere, which has 
as stated in Sec, 2,2.2, a thickness in the range of 50-200 m, 

Ih® plvaae of smoke ooiLtted from the moderately high chimineys in 
the surfao® layer would spread vertically as well as laterally 
and after travelling some distance doTOwind may spread into the 
re gie n above the surface layer of the atmosphere. Hence downwind 
dlstancfs selected Ibr the study were relatively ^ort. 

It was decided to use hydrogen gas-filled tethered balloon 
to measure the vertical temperature profile and collecting samples 
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of ftix“ at diffsrsnt hsights* As no othsr oicroiiietcoiological- 
infoxTsatloa was availabls for ths urban arsaj vane anemoiaeter 
was used to record wind speeds# Thus, instrumentation of case 
thre®! as described in Sec# 5#1> i.e# temperature records at two 
heights and wind speed record at one inteimediate height were 
obtained, 

Bia observed concentrations of sulphusr-dioxide at different 
points were compared with those computed using procedure outlined 
In the preceding chapter# 

6*3 Materials and Melhode 
6.3.1 at Souroe 

(a) attte One 

The cluster of five chimneys of the Riverside lower- 
House wm taken as the source of power-house 

has a power generating capaoi-ty of 1.1. million units per day. 

It has 15 boilers, fad on coal from Ranlgang series of coal mines. 
Ibble 6.1 gives the details of the chimneys. 6.2 gives the 

I . ■ 

partioularlB of the boilers. 



!EAB1E 6.1 




DEIAIIS OP CHIMNKrS OP RIVERSIDE POWER-HOUSE, KANPUR 


Gr* 

No. 

Height, 

ft 

Dianueter, 

ft 

oype 

Boilers attached 

1 

186.33 

10 

Brick 

1»2,3,4 

2 

204.5 

12 

Brick 

5, 6, 7, 8 

3 

204.5 

12 

Brick 

9,10,11,12 

4 

213.25 

8.86 

Steel 

12A,14 

5 

200.0 

8.86 

Steel 

15 


32ABIE 6 .2 

lAMlOUlAlS OP BOimS 01 RIVEESIDB POWIE-HOUSE, KANPOR 

Sr. 

No. 

a^pi' 

Steam-fbw 

rate,1bAr 

Ooal-bumlng 

rate,ton8/hr 

1 

liE 

20,000 

1.5 

a 


20,000 

1.5 

3 

w 

20,000 

1.5 

4 

I® 

20,000 

1.5 

5 

w 

40,000 

2.5 

6 

I® 

40,000 

2.5 


w ■ '■ 

40,000 

2.5 

& 

I® 

60,000 

2.5 

9 

HP 

100,000 

6.0 

40 

HP' 

100,000 

6.0 

11 

HP 

120,000 

6.0 

12 

' tan 

, ■ Ws ' 

120,000 

6.0 

t2A 

V HP 

120,000 

6.0 

14 

HP • 

120,000 

6.0 

15 

IP 

120,000 

6.0 
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sulphuz" content of the coal used was 0.4^, as given 
from thu analysis of coal supplied by the power-house authorities. 

On the days of experiments, the boilers nos* 1,5, 6, 9 and 14 
wejTQ not operating, but remaining 10 were in operation. 

(b) Site Two 

!]Sie ohimney of the boiler-house at the end of the Southern 
Laboratories, I.I.T. Campus, was used to emit known amount of 
BUl^imr-tlioxide. Sulphur-dioxide was produced by burning sulphur 
ut the bottom of the chimney. 

Eht chimney height is 15.2 m and diameter is 75 cm. The 
rate of bimiing of sulphur was maintained at 5.0g/min. 

6«5«2 Vertical OMperature Measuarements 

As described in Sec. 5*1 , different sets of micrometeoro- 
logioal Instnmients can be used for estimation of spreads of the 
pollutant, but it was difficult to go for continuous wind speed 
and dlx«otion records. Hence it was decided to go for the Case 
Three of Instrumentation. Tbls requires temperature measurements 
at two hei#its* Ibr this purpose, a hydrogen gas-filled balloon 
along with a temperature sensor was used. 

6. 3. 2*1* Materials ■ ^ ’ 

(a) Balloons 

Hubber balloons, specification No. TWiite S.R. 875, of 
Swastlk Rubber Iroducts ltd. Iioona-5, were used for this 'study. 
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Th^ae balloons when fulljr inflated with hydrogen gas assume 
n dirifflotur of about 150 cm. Ihe price (l973) of the balloon 
waa Ra* 30/~ ouch. Ibr each fillingj it requires about 1.5 
cubic motors (at atmospheric pressure) of hydrogen gas. 

Ilydrogun gas obtained in cylinders from Indian Oxygen Co. at 
the coat of Re. 350/- per 100 cubic meters was used. 

(b) Hwnnistors 

A thorraiator was used as a temperature sensor. Different 
tyjjes of thermistors are available for the different types of the 
use®. Thermistors were purchased from l//s. Tempo Semiconductors 
(18, Auronjape B Scheme, Subhash Boad, Vile Parle (East), 
Bombuy-57AS)» Two different types of the thermistors were 
triad* 

(i) Rxigged beads (Price; Rs. 12/- each) 

(ii) Glass mini-probes (Price; Rs. 20/- each) 

6*3. Methods 

(a) Billing the Balloon 

Commercial hydrogen gas was used to fill in the rubber 
balloon. A pressure— laibber tubing was attached to the hydrogen 
gas cylinder and the other end of the 2-m long tube was put 
Will inside l&e balloon, and tightened wftia a piece of cord. 
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oolilurod tt) thy thynaiator, were about 100 m long. Ihe 
mittrt5-iUiU.iutor wao used to measure the unbalanced current, and 
wttO protuctud by putting 1 kilo-ohms resistance in series. A 
9-vult batteiy powered the circuit. 

®ie oalibratlon of the meter was carried out in the 
lioouTfttyly-iiifilntained hot water-bath of a Warburg apparatus. 
Oic thujiaiutor probe and a mercury thermometer readang upto 
0.1® 0 wore put in a test-tube kept in the hot water-bath. 
Unbfilnnoed current passing through the ammeter was correlated 
with the temperature. 

Out of the two different types of the thermistors, 
it W'io found that the rugged-bead type was more suitable for 
the out-of-door work as it was in the present case. 

(l) Field Test Layout 

Once the themistor was calibrated, the thermistor 
aeneur and the lead wires were connected to the nylon cord 
of the balloon. As it ves also proposed to collect samples 
of air at different heights, a plastic tube of 3 mm diameter 
was fileo attached to the cord. ’Sam the combined cord was 
made up of the plastic tube^ two insulated copper wires as the 
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lund t{» thursaiator circuit, and a nylon cord. All these were 
ctsnnuctud by light transparent adhesire tape at every 30 cm or 
luiUJ. 'Bigu carrying distance (or height) mark from the probe 
(and thu mouth of the plastic tube) were also attached to the 
combined cord at every one meter. 

file nylon cord attached to the inflated balloon was 
OKtonUod upto about 2m and there it was tied to the combined 
cord# ®iw balloon was slowly allowed to rise lifting the 
oumbinud cord along with it* If the atmosphere is not very 
turbulent and the balloon is properly filled in with the hydrogen 
gas, the away of the balloon is almost negligible. Ihe balloon 
aaoande straight up# She box containing the complete circuit of 
ths thaniiistor was kept on the ground (and so was the instruments 
for oollcoting seoplee of the air through the plastic tube). 

'Cho balloon was slowly sent upwards. As the thermistor 
probu (and also the mouth of the plastic tube) reached specified 
heights, the balloon was kept floating at that height and the sample 
uf air wiw CO Hoc ted and the reading in the Wheatstone bridge 
circuit was recorded. She same procedure was repeated at different 
heights. . ■ — : 

6.3.3# Meaeurment of Wind Velocity 

A Vi.uiu anamoneter , was used to record the average wind 
vwlooily at a specified hei^t. Hie recording was done in a cumu— 
latire wacr for every 15 miautes. 
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^’*54. MgaaureriBnt o f Suliiiur-Dio^de Goncentratinn.. 

(a) Collection of Sample 

An It wao decided to collect samples of air at different 
lu'lf'htu, a plastic tuts of 3 mm diameter was sent up along with 
tlio to the .red balloon, as described in the preceding section, 

A hand-operated midget vacuum pump was used to suck 
nrtfflploo of f.iir. The absorbing solution was put in a midget 
liiipingor iind ttie rubber tubing provided at the end of the 
pio.Htic 'tube Web connected to the inpinger. The vacuum pump 
connoction was trade at the other outlet of the impinger. The 
vwtmm pjttip wan re-oalibrated to know the rate of flow of air 
through the entire length of the plastic tube (as the head loss 
iiiorRaaOB with the length). The pump was operated for 10 
ffiiuuten at each sempling. As the rate of flow and period of 
sampling wore known, the total volume of air sampled could 
bo calculated. 

(b) Analysis 

e 55 ) 0 sed absorbing sample was analysed by the standard 
V/oMt and Gmlm method which is described in Appendix A-2. 

6.4* Results and Discussions 

6.4.1. Zero-plane Displacemt Length ■ 

(a) tut© Ono 

The average height of buildings axound that area 
in abtfit 40 ft « 12 m, mgking h » 12 m in Eq. (5.l)v ISie value 



175 


u{ A it> L(i. ( (. 1 ) ia the ratio of total axea of a given region to 
Uu- niVH of tho iTui Wings in that region, lesser the density of 
hoi liiingn, highur the vaWe of A. It was ohserved that the build nags 
aroiunl UiO Mto Ono weie relatively sparsely situated, the ratio of 
bu:' It uj' area being around i2^ , hence value A was assumed as 9* 

'Hies value of d is given by 

■| *» (a) 

wlioi't) h * 12 m 

.% d « 12 X (9)“°*^ 

or, d « 6 in 
(b) Site Two • 

'Hio obatruotions were mainly trfees and single -storeyed 
IjulWingt) with average height around 6 m. The obstructions are 
not very thickly located the area occupied by obstructions being 
around 5-75^, and A value, consequently, was assumed as 15 . 

d « h(A)"°’^ 

-ecu)'®-’ 

a 3 m 

6,4.2. flu r face Roufhnosa Length 
(n) Site One 

From Bq, (9.4)> the surface roughness Imgth z^ cm be 


oalculatod* 
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» 12 X 

« 1.0 m 

(b) Sit© fwo 

Th® surface roug^ess. length is given by 

f 

* 0.4 m 

Tiioao values fall well within the independently derived 
vnluoD of for different roughnesses, shown in Table 2.1. 

6.4*% MioroBB teoro logical MesBuraaeaiits 

(ft) Site One 

As discussed in Sec, 54. 5> temperatures were measured 
at huit>iitfl z^ ■ d + Zq + 1 * 6 + 1 + 1 = 8 m and = d + + 4 

m $ + 1 4 4 ■ 11 m. Wind velocity record was obtained at hei^t 

zg » d + z^ 4 2 - 6 + 1 + 2 - 9 a. 

For the tomperaturo moasurements, the themdstor was 
oalibrfttod. The calibration graph of the current (in p A.) versus 
ton!|»<jfa1«r® (degrees, O) showed that for the range of the 
tomiK^ratiiro of interest, i.e. 10»0 to 50 “C, the plot was lineax. 
®i© slop® of the Un® was 0,1 ® G per PA. 
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IiVtry time when the thermistor was at 1 m above the surface 
Uu' ..■nrih, Uiu mnmeter was adjusted to the central reading of 
hO A. Iho nnrronponding actual temperature at that height was 
lii't.fd down \/i th tho help of a mercury thermometer (reading upto 
0* r’C), When the Ihermistor sensor was sent up in the air, the 
In tho ammotor reeding was noted down and was converted 
iiiiu dlfftirt^nau in temperature in ®G, This could avoid the use 
of tinUbrntion onrvo for every time converting the ammeter 
rc.Twlitig to oc'rroHixmding temperature reading. 

'IVinporntu ro readings for 8 m and 11 m were interpolated 
from tomj)t>rafcura profile obtained upto hei^t of 40 m. 

(h) Site IVo 

Itio temperatures wore measured at heists = d + + 1 

t» 3 0.4 1 w 4 *4 m, and » d + + 4 5 i* ^*4 4 “ 7 *4 

whurofio wind velocity was measured at a height of z^ = d + z^ + 2 = 

3 h 0,4 + 2 - 54 ». 

At each Bite, three different runs were carri.ed out at 
throe: dlfforemt dlstanoos from the source, as shown in Table 6.3. 
lArring o/joh run, oonoentrations and temperatures at three different 
henghto worci moasurod. Stability conditions were different for 
thono throe runs of study. Table 6.3 reports the temperatures and 
wir^ speed m0a»J'r®mehts , 
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TjfflLE 6 .5 

\/INi) ..IW'iD /iND VERTICAL LEMEERATOBES FOR KAWHJE STODY 


Run 
Hu • 

Downwind 

diatonoe 

X, m 

Temp, 

@ 24 m, 

®C 

Temp, 

@ H'x m, 

=G 

Wind 
Speed 
@ !^m, 
m/ sec . 

Period, 

hrs 

(n) 

Site One 





1 

7OCI 

1940 

18.60 

0.780 

ho on 

2 

400 

16,90 

16.84 

# 

0.34 

Late 

afternoon 


200 

15.10 

15.12 

0.122 

Evening 

(b) 

SI to IVfo 





4 

100 

28.15 

28.025 

0.480 

Noon 

5 

200 

27.60 

27.63 

0.360 

Late 

evening 

6 

100 

32.50 

32.51 

0.265’ 

Early 

evening 


aibstituting those values in Eq.. (5.9)j the values of 
otitWlity parOTOtor B were calculated; from which, using Eg.. (5.10) 
vuhion of stability length L wore derived, as reported in Table 64. 
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TABLE 6.4 

VALUED 01-' STAI^IUTY PAEAMETBRS K)R KABHIR STUDY 


('0 

Run 
No • 

Site One 

B, 

-1 

m 

L , 

m 

(b) 

Bun 

No. 

Site Two 

B , 

-1 

m 

B , 

m 

1 

0,C42 

-70 

4 

0.0178 

-50 

D 

0.022 

A 

o 

5 

-0.0075 

70 

3 

-0.016 

45 

6 

-0. 001.6 

115 


6 «4 4 » Cnloulntlon of Vortical Spreads 

From thB oatimated values of and L, the vertical spreads 
otttn be oaloiilatod from Pigs. 5*1 and 5*2, for appropriate values of 
X. Hurt} , being 1 m, the figures 5*1 aud 5*2 represent x and z 
vuluou^ fur Site One. 

Tli0 results are shown in Table 6.5* 

TABLE 6.5 

VERTICAB SPREADS z POE KANKfR STUDY 


(a) 

Bun 

No, 

Sit® One 

mm 

ZfXa. 

(h) Site Two 

Run 

No. 

■ 

1 

115 

4 

20.0 

2 

62 

5 

11.2 

5 

15 

. .6 

9.2 
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of Latoral Spreads 

/^iurmling to tho stops given ip. Sec. 5.6, the lateral spreads 
Vv'nt calculatu.U Tito remits of inijortant steps are shown in Table 6.6. 

TiiBIE 6.6 

LA'ffiHAL SPEBiDS 0 FOR KfflHJR STUDY 

y* 


^ • 

£ i 

m / Soo . 

Height 

K, 

m 

m/ 800 * 

radians 

T J 

min 

X 

V 

m 

B — 

u(z )8 

(a) 

Site One 







1 

0.461 

20 

2.890 

0.558 

10 

60,5 

650 



30 

3*300 

0.314 

10 

55.0 

595 



40 

3.340 

0.292 

10 

49.5 

575 

2 

0 . 17 R 

20 

1.160 

0.345 

10 

86.2 

550 



50 

1,360 

0.294 

10 

75.5 

294 



4O 

1.475 

0.282 

10 

67 .8 ■ 

288 

3 

0,122 

20 

0.575 

0,200 

10 

87.0 

96 



30 

0.810 

0.142 

10 

52.0 

77 



40 

1.060 

O.1O8 

10 

47.0 

61 

(Ij) 

Site 'Two 







4 

0,123 

10 

O.9I8 

0.300 

10 

27.2 

102 



20 

1.250 

0.220 

10 

20.0 

79 



30 

1,450 

0.191 

10 

17.5 

71 

5 

o.an 

10 

0.700 

0.260 

10 

71.5 . 

150 



20 

1,050 

0.174 

to 

47.6 

97 



30 

1.325 

0.158 

10 

'37.8 

80 

6 

O.C 62 

10 

0.500 

0,279 

10 

50.0 

76 



20 

0.717 

0.194 

10 

54.8 

57 



30 

0.870 

0,160 

10 

28.8 

54 
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^’* 4 *f>* QoMiy ttJ-tlon of Oonoontrationa 

I'irntiy, fcho ^luantity of SOg emitted at the source is computed 
tu onch cnnuf and rop(.irted as Q,, g/seo. As the distances at ■siiioh 
Llic concuntrationa aro measured ars Tery diort, the correction for 
the plumo rioo nood not be applied (Burner, 197O). The actual 
hu-lf]:htn of tho ohimnoys are taken as the effective heights. 

For Si to One, tiie total qnantities of fuel burnt on those days 
VKU'kwi nut to bo 59 #5 tonnes/hr. Taking 0,4^ as the sulphur content, 
fi 4 nr5Uiin.ilHr weight of 302,52 molecular weight of S , the quantity of SO2 
being jflumiod per second works out as 

^ (g/tonnes) x 0,004 (g/g) x 64(mol.wt of SO2) 

Q M ”• J ■" ^ J2 (mo i . wt . 0 f ' S ) ” 

m OM (fj mo of 302 * 

For 3 it<j Two, the 5 g of 'sulphur was burnt per min. Hence 
5(g/g) X 64 (mol. wt. of SOg) 

^ *■ 1 (min} X 60 (seo/min) x 32 (mol. wt. ofs ) 
w 0,1^7 g/o®c of SOg. 

Bofwlth of IS10 computed and observed concentrations axe 

/'Ivuu in Ttihlo 6*1 * ■ ^ ^ ^ 

.‘Sfeuaplljug stations were seleoted as near as possible to 

tho mtloipfttwd mem oeter line? hence, y « 0, for each sun. 
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TiBLB 6.7 

('Olil.ll'f'iJ'IT CONCfclN'HlATION VAIUES FOR KiNFUR STUDY 


Ituri 

No, 


f. , 

m 

Observed 

concentration, Co, 

Ppn 

Computed ' 
concentratio n, GO , 
ppn 

olo 

0 [0 

(-0 

HJi U» 

Ono 




1 


20 

0 ,C 420 

0.0427 

1.020 



50 

0,07!3 

O.C42O 

1.120 



40 

0.0580 

O.C417 

1.100 

2 


20 

0.3170 

0.2780 

0.877 



50 

0.3620 

0.2770 

0.766 



40 

0.2800 

0.2760 

0.986 

3 


20 

0,1510 : 

0.1620 

1.070 



30 

0,6000 

0.7650 

1.280 



40 

1.4210 

2.5900 

1.680 

(b) 

Si isj 

Two 



>' *' 4 

4 


10 

0.0067 

0.0075 

1.120 



20 

0,0055 

0.0065 

1.170 



30 

0.0050 

O.OC47/ 

0.950 

b 


10 

0.0120 

0.0098 

0,816 



20 

0.0105 

0.0093 

0,880 



30 

0,00 00 

O.OC4O 

- . 

6 


10 

0.0350 

0.0293 

0.890 



20 

0,0312 

0.0290 

0.950 



30 

o.c^oo 

0.0790 

0.880 
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ft. 4. 7 PiscuBulun of Rosulta 

thu vciluo 0 of pollutant concentrations given in 
'Inhlu ft .7, thu following points can be observed: 

(i) toaed on the proposed formulation, the vertical and hori- 
aprucida and lienee the pollutant concentrations can be 
CHloulutud which are generally comparable to the observed concen- 
trntifsrin. 

(li) During unatable atmospheric conditions, for both the sites 
of InvutJtignbton ysonoentrations computed agreed very closely with 
nbaurved tinwa, Ih® variation of the ratio of computed to observed 
conoentmtione was within the range of 0.766 to 1.17, which for 
the Mtntiaticul variation of the entity may be considered as 
roaaotuibly narrow. 

(ill) Dtublw at.ioaphuric conditions gave generally close agreement 
botwuun cwaputud (.umi observed concentrations for the area of low 
nurfneu ivughnuss like Site fwo was. Bae ratio of computed to 
obnurvud vnluuu of concentrations was between 0.816 to 0.93. 

(iv) nt tho area having large surface roughness, like Site One, 
t^m ttgru*Atunt botwuon computed and observed concentrations was 
rftthsr poor^ during very stable atmospheric conditions. (Hie ratio 
of uiMuputcd to obsurved concentrations lies in the range of 1*07 to 
1*68# As diaounaod earlier jthis may be regarded mainly due to 
imuleqiwit# wind profile formulation during very stable atmospheric 
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0 tiji'i 1 til tUt! fifui uImo (lull to effscts of urbcm hsating. 

Al'j’li'- ability uf the fo.xo,ulatioa during neutral atmospheric 
wutttllti'UUi not bo verified in this study at Kanpur j as none 

of tiiu uxptiriniuntril run was conducted in such conditions. 

Fttju this study, it appears that the proposed formulation 
lii.’iy hu uuwd to prod lot ooncentrations at any point downwind from a 
iuuui*t!w uitbor at or near the ground-level. !I!he procedure used 
ijuru to ouinputu conoentrations of pollutant requires the micrometeoro- 
d.'itn which cm be relatively easily collected. The procedure 
inouri'Mrntuo, in n quantitative manner, effects of different svirface 
roughnoasoo and of different atmospheric stability conditions. 



Chapter Seven 
CONCIiUSlOIS 


I/igrangiun slinllarity theory has been applied to diffusion 
pliunvii.iunun to predict the vertical spreads of pollutants released 
frtji.i n iiuurco situated in the surface layer of the atmosphere. 

litijjvd on the investigations given in the preceding chapters, 
tho To Hawing oonoiusions can bo drawn j 

(1) hxtondwd Ifttiiin’s Approach to represent plurae boundaries by 
hia uquationa generally yields lower estimation of vertical 
sproad® for all atoosjJierie conditions, and is hence not reliable. 

(il) Matbjuatlnal equations can be developed to represent the 
vurticttl tjprond of tho mean of the ensemble of particles at any 
noon tlovmwind distance from the source, as has been done by the 
Huthiir in Chapter 'fliroo. • 

(lii) Bor both low as well as high surface roughnesses, the 
foTOuliitiurw proposed by the author give very good predictions of 
vurkloul spruods In unstable as well as neutrally stable atmospheric 

condittona. 

(iv) for B table atmospheric conditions, the formulation gives 
couparable results upto i/l Beyond this level, 'th.e 

divergence between observed vertical spread and computed vertical 
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usiruacl buouwur, conoldwable. This is piobably due to inadequacies 
in th«i I rtjuunti-tiny knowlodge to express wind profiles mathematically 
in vury ritnblo utauapheric conditions. 

(v) Unliko other available fonnulations jthe proposed mathematical 
ioOTulntinn (nnd tables and charts prepared therefrom) can be 
nuccuonfully used to predict the quantitative effects of different 
nurl’nou i%<U|?luiuBHyo nnd different stability conditions of the 
ntuonphum. 

(vl) Itiio fbnaulution oan be used with even minimum micrometeoro- 
logifml inutruauntution (viz., temperature records at two pre- 
do tunuiawd heights and wind velocity at one height) without recourse 
tu thu uuu of u nutiber of empirical coefficients, as is necessary in 
variout) oxinting funaulations. 

(vii) Very little micrometeorological information is presently 
available which can be useful for the prediction of concentrations 
of pollutaritB oven in urban and industrial complexes where air 
poimtii<n jxstentlal is veiy high. Itoring planning of air qualily 
mnniigarttfit of a region, along with instrumentation for monitoring 
uxtunt uf air ixjllutlon, instrumentation for collecting micrometeoro- 
loglcnl data nhould also be planned. Instruments, for recording 
turbulenoe fieM would be ideal, but if they are not available, as 
shown by tbw present investigations, at least instruments for 
recording wind velocity, wind, directions and temperature profile 



fihuuid Uu itiotnliwwl, fxnd such records must be routinely collected 
fiir nix Wiu urbmi and industrial areas* 

Lliuitntlono u,t tho Irojosed At.'pxoach and Direction for Future Research, 

Ihu major limitation of the proposed formulation is that it 
a|.ipille8 strio tly in the surface layer of the atmosphere* Althou^ 
many of thu «X)uroos of nir pollution are located within this layer, 
t.nll ohiJtmwyw of induntrtol complexes may emit pollutants in the 
rwMlon beyond thia l/iy^T and diffusion phenomenon for such cases 
douu not lend itaulf to ©nay mathematical formulation, 

Hiu proniawd formulation applies to the case of a single 
point auurew. ^ extend its use to a multiple-source urban 
diffusion problem would be a little unwieldy, because contribution 
from each aource located upwind has to be calculated separately and 
funaulutlon proj-osed here does not relate x and z by a simple 
foai. 

Hencjo, it would be interesting to investigate, in future, 
thu uawfulneoa of the proposed formulation to the study of multiple - 
yuurou diffualun problem, and to suggest, if possible, simple 

M ' .mm , 

mlutim butmmm x and i # . 

probliffli of diffusion in., vsry '-S'tebls a1iiiosph'6Tic. ■ co.n..dit,ions ^ 

alocj nciycii further inwitig&tions.» . Valid lao-thsiaaticai'forisiulatloii 

fiir tl'io vertioal wind profile during such stable conditions of 
atoKJSphere is very much essential fop ary work of the type proposed here 



Mare uxi’l0rtttli)n ia roquirsd for this purpose. looking to the 
fact th.afc Uu< miy stahio atnospheric conditions impede the 
tja|fusitiy ui thu air to disperse pollutants and hence the pollutant 
conewitruttons build up, adequate mathenEtical relations are 
ruquirod to predict concentrations during such conditions of the 
atimaiihurw for th® successful managament of air quality in any 



MIMDIX A-1 


COMiSi fi-at iHOCaUiMME TOR NtMERIGAI. INIEGBATIOI USIHG 
GAIIS3-LBQMDRE QUADRATURE 

C C1M!T4»HITO thbory 

n SODVIHG (DX/DZ) NlJMERIOADEf 

c a;jn;:i-u-;ra'3n)iiK gu/ajHATORE WITH TO® TOiiTS 

(1 IN'P-lNTKHVnL 

G W’-FirinL l.lMiT 

0 OF 10 mu IN GAUSS-LBGBNDRE QDADRA!roRE=4 

C N;I-'NI,I14I»M( of STEI’S 

0 A»INiTlAIi LIMIT 

C Ai* AH FA OF INTERVAL, OR lOINAD 

U i\N«NEW AREA 

C NTO««NUMBHi OF ZO/L 

00Mi40K aI,I 

DIMENSION Z0L(10),AI(10) 

RKiJ) 100,XNT,IJ’,M,N20 

1UO TOH14A'T('.I> ^ 

UBmD 110,(AI(I),I»»1»NZO) 

no tormat(?i7.3) 

R BAD 1 20, ( ZOLC I ) » I»1 »1IZ!0 ) 
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r^O K)HMnT(Vi'7.4 ) 

AO '■■{).() 

NS«I,1?/INT 
DO 2 IftIjNZO 
iklNf 200,E0I,(I) 

200 i-X5Hi4AT (li!,1S>X,5H!f50/lP,2X»F7.4) 

IHINT 201 

201 TOHMAf (1H,10X,3H^I,,22X,5HKBX/L,15X,4HZ/ZO,15X,4HX/Z0) 

IX/ 1 1,100 

H^AtFTXiA'KlNT) 

AHtlAoO/AIODCAjBiM) 

/iH«*AO OUjEA 

IF {j.EQ.lOO) 00 ID 10 

AO--'^AN 

,A*B 

«) TO 11 

10 A«0.0 
A0«0.0 

11 H1*»D/(W)l(l)) 

AI<AN/(2OIi(l)))/0.48 
ifim* 202,D,AN,B1,il1 

202 raMii.TC/1H,5X,f10.3,10X,FI5.6,10X,f15.6,10X,P15.6) 

1 COOTIKUE 
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2 COMTINtli'] 

r.Wi' 

HN!) 

RmCTlON roNCTN(X) 
cxmwjji Ai,i 
lHM}«ai'*'nN Al(lO) 

HIM* : W -( A1 ( 1 ) -1 . 3/( l**i n/3 . ) ) )/( ( 1 . +7 . )** ( 1 ./6 . ) ) 

HK'ilHjll 

P") 

HJMliTluN OAimfl (A,B,m) 

C-(H-4i,)/2.0 

tmM) 6'’>2l4*(KirfC'ni(0.33998»^+D)+mJNCm(-0.33998*0+D)) + 
1 (»,''-47l'.»(«lli0TN(0.B6113*C+D)+KJMCm(-0.86l13*C+D)) 

j.ujsn 

HKTOHN 

END 



iiPFEIDiX A-2 


IfliJ’ilJjJiD V/FiLi'l' ijND GiEKE IffilRQD FOR IffiASUEBMBNT OP SGLHIUR DIOXIDE 
lii AIH (TnULan Standnxd IS i 5182 (Part II) - 1969 ) 

1 • }kuif':titit.n 

(-0 if'Uorbing aolution 0,i M sodium tetrachloxomercurate* 

Dissolve 27.2 g mercuric chloride and 11,7 g sodium 
chloride in 1 litre of distilled water. 

(h) p-Roni«iiline hydrochloride. 

Dianolve 0,20 g of p-rosaniline hydrochloride in IOO ml 
of’ liuitliled water and fitter the solution after 48 hrs. 
!‘l}>ott ,0 20 ml of this into a 100-ml volumetric flask. 

•Vld 6 ml of Qoncentrated hydrochloric acid. Allow to 
i t.fuul for 5 udnutea, then dilute to mark with distilled 
w ntci r* 

(o) |'‘o'rmaldohyde solution. 

Dilute 5 wl of 40?^ formaldehyde solution to 1 litre 
with distilled water. , 

(d) Bodiuffl sul^ito ijolution 

DiBBOlve 640 mg of sodium me tahi sulphite in 1,0 litre 
of water. Standardize. 1 ml = 150 Pl of So^. 

2 • I'rooednco 

(<i) Pitipfirution of standard curve 

(l) Pi, potto exactly 2 ml of standaiii ailphdte solution 
inh) u lOO-wl volumetric flaic and dilute to mark with absorbing 
'Bil 0 final solution oontained 3*0 vl of Sog per ml. 
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( ii) Md ticcuratoly 0,5, 1,0, 1*5 and 2.0 ml portions of the 
di lut... nilphite solution to a series of l0-ml volumetric 

i (uul di.luto to the marks with absorbing reagent. Continue 

with tho iBinlyais procedure given in 2(b). 

(iii) Plot the absorbance as the ordinate against the ^ of 
per 10 ml of absorbing solution. 

(b) /jruilyoln 

(1) Ooi up a uoinpling train consisting of absorber, metering 
tiiiviru, ond air pump in this order. 

(il) I’lFU'ttd uKaotly 10 ml of absorbing reagent into the 
abttorbor. /,a{d.rato tho air sample at a rate of 0.2 to 2.5 litres 
poi* min. id the end of the sampling time, stop the pump and 
tho ttbaorber for titration. 

(iJi) !M 1,0 ml of p-rosaniline solution and 1.0 ml of 
rtumaWohydo solution and mix well. 

(iv) /i’tor 20 min read the absorbance at 560 my in a 
npuiiltn photomotor with a lO-ml unexposed absorbing solution as 
Uu- bl.'iJik for roferenco, 

(jfllcu lation 

(i) Cancan tration of Sog in ain is s 

yl of SOg 

Sog, pp by volume « volume of air sampled 
(I I.) Sag, Pinn by volumex- “ SOg, yg/m . 
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A9&ENDA 


Replies to thtt oommts of Thesis Examiners 

Bximiners* Report No. I 

CD The point observed is corxect» as z in Bq,f (2.4a) is from the 
EartA*s surflaco, whereas * in. Eq. C2.4b) is from a basis where u.s 0^ ' 

C2) Although it is true that Pasquill (1966) "has proposed that z in 
equations (2*50), (2,S2) and (2*53) is the mean height of an enseable 
of particles, he could not get validation of his proposal in stratified 
flow. Thus, Pasquill (1966) concluded as : 

♦If the apparent failure of the current application of 
fiailaiity theory in stratified flow does indeed arise as st^gested above, 
it if difficult to see how to represent the vertical spread in terms 
of I alone, A possible alternative is to revert to the type of 
representatim used by Monin (1959)*’ 

In view of the observation made it is suggested that line 10 of . 
the tynopsis b© road as 

♦Besides, unlike Nbnin’s approach, vdiere he assumed that the 
equationi of motion represented boundaries of the plme, the author 
assumes, following Pasquill ’s suggestion, that the equations of 

motion bo token to »p« 80 it tlw mem height of ensemie of perticles.' 
(S) The lnt<igrel In Bq. C3.S5) we evaluated using the Table of 
Integral Tranafoma (Erdelyl et al,, 1954). It is correct that laplace 

TranaftBaa aa auch wre not uaed. 

(4) lh» valuea of J/C-« 9“*' 

b. d.nv.d fr» equatlona C5.«). C5.45) m4 C3.47). But within the 



2 


sngineorlng limitations of accwacy of measurements of various 
partatttrs, the valuos of I/C-D derived from the graphs as shown 
in Fig. 3.6 will b® nomlly good enough. 

(5) All these misprints ax® noted and corrected, 

(6) Within tho mpinn of 0 < i< the velocity profile is undefined 
,nd hmc. it is not theoretically correct to evaluate the integrals 

as macatod in oc, nations just after C3.30). C3.34b), C3.39), (3.51), 
(3,54) from the hounary 0 to hut they should have been Integrated 
fn» Ig to However, the region of 0 to t, being very small, the 
dlfferloe in the numerlca vaue. is not expected to be significant. 


Al 


suggasted by the eamlner. the evaluation of lnt..gral in 
the ren*. of to - i. arrled out below for n.utr.1 stehillty 
conatlooe fcr equations ftllewlng Bq. (3,30) 


/ c S d* 

Ak ^ *0 

It Jcds 

K 

m 

% 


f Cq ®3CP C- r-) 


Writing * ■ *’ + H» 
it •• tE», 


. snssuring height from s. basis and denoting 


/** ©39 C“ 


Z’+Zq 

% cm -— •) 


^ r- 

dt h 
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Table 


1 

0 

m 

E-. 

*0 

(new) 

m 

2L 

*0 

Cold) 

m 

z 

*0 

m 

(new) 

■ ' ■ ■ 

2o 

(old) 

1 

1,94 

M 

50 

764 

740 

S 

MM M 

Ml 

60 

976 

950 

10 

m 

ss 

70 

1205 

1180 

IS 

145 

114 

80 

1145 

1405 

20 

217 

187 

90 

1690 

1650 

30 

579 

3S0 

100 

1935 

1890 

40 

ses 

S36 





Thfii cticaiittlanf ittidacat® that for larger values of z/Zo,(such as 
440 in the PP6 itudiei) tho values of 2 /Zq obtained from the procedure 
followd above will be within II of the values computed as in the 
thesiSf 

Thus, in Table 4*2 of thesis, the values in columns 3,4 and 5 
would ramtiii «nchaii|id. 

In the east of other stability conditions, both those stable 
and aiatabla mm, the influence of this aange would be even smaller 
and nagleeting it would be permissible. 
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Exawiners* Report N0« 2 

( 1 ) Th© integration on page 6S was checked again and it was found 
that the rwsult* reported in the thesis was correct, 

(2) 1h© auggestion is good^ but the assmption of 95% conversion of 

S to SQ^ ii nlao arbitrary as is of 100%, 


Examiners* Report No, 3 

Under th« heading 'Connents’ i 

CD Th, v»lu«, ot a liapand on tha jtablllty conditions. 

Tor miUbU »t*>sph.«. O'Brim C1965) snalysod over 2800 

wind ptoniai »■» 

M .ugj«t«l otlglnauy by H>nln-0b«khov. Tbls value was adopted 

for the preawftt atu#. 

Por ftable atnoaphtre, however » the variation in the values 

,, by diyyermtautbontlesvas aider, m lOSS. S^nbenB 

ft A hnt later Thylor i%9m reanalysed the nunerous 

ttt 0 • 0ii| few ImXMT v. 

alnd rrcflla. raported by S,lab-.y. Mcnln-Cb^bov. etc. md arr^d 

1 ii ftf a • 6 fbr stable atiotphere, which is adopte ere 
at ^e value of « ■ « 

ft>r the present study* 

(2) Eoimtim 2,15 should r«t<?‘ 


u* 


U . P- t-cC- «/«■''** 


S.S.1 la *l»» » W" «-’*• 
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Wn^t thi htttding ’Concentmion Distribution* : 

CD It it ttw that PwKjuilX suggested a general profile of vertical 
diftributioti of conc«itmtlon as exp («b z*)T'however, the values s has 
not been uniquely offbrred. Pasquill (1962, p, 188) himself observecU 
*Vlieit the tourcti is at ground level the vertical distribution 
evidently requitet a dlfiiriiit value of a, which has yet to be generally 
^temiiiid, hut liildli on preliminary evidence is nearer unity, ' 

iotMiquet4'e arson (1936) and Calder (1952) also observed and 
adopted f • i« tm such the simple exponottial distribution with $ - 1 
itas adonl^dL 

Sli^larlyi the asiuarption of (Saussian distribution In the 
hotitontal is also eaplsieal, hut there is more evidence &r s nearing 
2, and henee (Saussita fittstrihution was adcpted, 

(2) In Table 4.2, values of i/C-D, i/^o or i/L were read from the 
appropriate piphi# Prom these values, a was calculated. To cospare 
these fxxsputed t with observed Gaussian 02 values, z was converted to 
Oji by Sq, (4*1) assuming Gaussian distribution. Thus the last coluon 
in Table 4.2 gives '^e ratio of oonputed z/0,8 and observed Gaussian 

H ||t 

Althmijih there art several models existing, ncne of them 
ineei^fetes the effects of both the inhomogeneous turbulence field 
of dif^r^it stibility conditions of the atuosphere as also of 
different aerodynasdc rou^esses of the underlying surface. In me 
model propoted here, attempts have been made to incorporate both these 



and in tUa differs frem the other 

tdcliting i04iltt 

C«) The anparUNSta 4fntid mTSI Kanpur uaraireally meant to he 

CMTily pilot study j to ussoss ttie po8si1^l!|i^df using the p^oposod 
tJOTiulation to t point scurc© at moderate height above the ground 
lUffhct. Thort wtrt many shortcooings in the expefimental work* 
tht vtnlcal ptofllo of teiftperature and measurement of pollution 
concentration at different heli^ti were difficult to get. The 
rtcourae was to use meteorological balloons which are helpful for the 
purpose only during relatively low wind speed conditions. 

the nap of location of powerhouse in Kanpur city is attached. 

The hadtground ctmcintration of SO2 was negligible in both the 

•I tea* 
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